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Abstract 



In recent years entanglement has been recognised as a useful resource in quantum 
information and computation. This applies primarily to pure state entanglement which 
is, due to interaction with the environment, rarely available. Decoherence provides the 
main motivation for the study of entanglement distillation. A remarkable effect in the 
context of distillation is the existence of bound entangled states, states from which no 
pure state entanglement can be distilled. The concept of entanglement distillation also 
relates to a canonical way of theoretically quantifying mixed state entanglement. 

This thesis is, apart from a review chapter on distillation, mainly a theoretical study 
of bound entanglement and the two major open problems in their classification. The 
first of these is the classification of PPT bound entanglement (separability problem). 
After having reviewed known tools we study in detail the multipartite permutation 
criteria, for which we present new results in their classification. We solve an open 
problem on the existence of certain PPT states. The Schmidt number of a quantum 
state is a largely unvalued concept, we analyse it in detail and introduce the Schmidt 
robustness. The notion of Schmidt number is exploited in the study of the second 
major problem, that of the existence of NPT bound entanglement. We reformulate this 
existence problem as particular instance of the separability problem, present a uniform 
treatment of distillability and numerically demonstrate that bound entanglement is 
primarily found in low dimensional systems. 
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Preface 



Quantum error correction allows quantum information to be reliably sent through noisy 
channels. Tremendous effort has been devoted to the study of quantum error correct- 
ing codes and capacities of quantum channels. This thes is d eals with entanglement 
distillation, also known as counterfactual error correction 13^ . The idea is to purify 



the noisy channel before trying to send any information over it, rather than the more 
indirect method of error correction. Personally I feel very strongly about the idea of 
distillation since it also relates directly to a canonical way of quantifying how much 
useful entanglement is present in a state. The distillable entanglement is a physical 
property of quantum states, bridging between a practical and theoretical study of en- 
tanglement. Bound entangled states are states from which no pure entanglement can be 
extracted, although entanglement is needed to construct the state. They are two main 
open problems in the classification of bound entangled states (the separability and the 
NPT bound entanglement problem) and they constituted the bulk of my research. 

During the course of my Ph. D. I have published a number of papers, but in writ- 
ing up, I have tried to make this thesis more than just a compilation or medley of 
these papers. The primary goal has been to convey the context and motivation of the 
different problems. Secondary, I have tried to provide a review of the state-of-the-art 
knowledge in entanglement distillation in the bipartite scenario. A notable gap in this 
presentation is a treatment of experimental entanglement distillation. The structure 
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of the dissertation and the interrelation between the chapters, is depicted in Figure [TJ 
The figure also contains a reference list of my papers, and their position in the thesis. 

Chapter one deals with the separability problem. After having reviewed the theory 
of positive maps and entanglement witnesses, I present new classes of PPT entangled 
states with particular ranks. Next, the permutation criteria for separability are studied 
in detail and a new result on their classification is presented. This is the only place in 
the dissertation which deals with multipartite systems. 

The Schmidt number of mixed states is a generalisation of the separable versus 
entangled division, imposing a finer structure on the set of quantum states. This is the 
topic of Chapter two. The tools for the separability problem are readily generalised to 
deal with Schmidt numbers. I review in some detail known results. In a final section, I 
introduce the Schmidt robustness and derive upper and lower bounds for pure states. 
The motivation for the topics discussed in this chapter becomes clear in Chapter four. 

In Chapter three I review the literature on bipartite entanglement distillation. Af- 
ter the necessary introduction of quantum operations, the main results on pure state 
entanglement transformation are sketched. Then I move over to mixed state entangle- 
ment. I start with a discussion on entanglement measures focusing on the distillable 
entanglement and end with a discussion of irreversibility in mixed state transforma- 
tions. This leads to the introduction of bound entangled states. The chapter ends with 
a review of mixed state entanglement distillation protocols. 

Chapter four is a discourse on distillability and NPT bound entanglement and 
describes my main results. As usual, the first sections introduce known results. In 
particular I review the three fundamental theorems on distillability. Next a uniform 
treatment of distillability theorems is presented. I illustrate the power of this reformu- 
lation by means of four applications. The most important ones are the derivation of 
reduction- like maps and an estimation of the volume of distillable states. In the next 
section the problem of the existence of NPT bound entanglement is reformulated as 
the problem of determining whether a particular set of states is entangled. It follows 
that the two major open problems on the classification of bound entangled states are 
really just one problem (although the problems are by no means equivalent). I end this 
dissertation with a short review on activation of bound entanglement. 
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CHAPTER 1 



Entanglement and separability 



88] . entanglement 



From the 1930s, starting with the EPR thought experiment [27 
was mainly studied in the context of the completeness of quantum mechanics. It was 
not until the 1960s when the problem of com plet eness was put on solid mat hematical 

Hunds by the remarkable works of Gleason [106], Kochen and Specker 175 ] and Bell 
, [l7]. The work of Bell dealt directly with the EPR thought experiment and is of 
enormous importance as it showed that entanglement is incompatible with a certain 
local classical inequality which can be verified experimentally. We will discuss the EPR 



argument and Bel 
found in Ref. 2221 



s ine quality briefly in Section Tl. 11 More in-depth treatments can be 



244] 



In recent years, the problem of completeness seems pretty much settled and consid- 
erable effort has been put into understanding the mathematical structure of entangle- 
ment. A primary problem is, given an arbitrary quantum state, to determine whether 
it is entangled or not. The study of this problem came together with the realisation 
that Bell -like inequalities are pretty weak tests for entanglement. The breakthrough 



132, 



223] came when notions from convex analysis and C*-algebras were applied to the 



problem. In Section 11.21 we review these mathematical results leading to the general 
framework of entanglement witness and its dual formulation in terms of positive maps. 
We also introduce the so-called PPT states and relate them to classical problems in 



1 



1.1 The EPR argument and Bell's inequalities 



2 



linear algebra. The separability problem has been shown to be NP-hard [llOfl , and 
therefore there is little hope of finding a simple analytical method which is able to 
distinguish all entangled states from all separable. Subsequently, interest has grown 
in finding good numerical methods in tackling the problem. In Section 11.2.51 we apply 
tools from semidefinite programmingto the separability problem. In particular we fo- 



751 ]. which besides a complete numerical solution, 



cus on the work by Doherty et al. 
offers appealing analytical features. 

In Section 11.31 we step away from the entanglement witness app roac h and study 



the permutation criteria for detect i ng e ntanglement, initiated in Ref. [2511 ] and further 



developed in Ref. 



63, 



144, 



254, 



3241 ] . Although in general, these criteria cannot de- 



tect every entangled state, they are particularly appealing because of their operational 
simplicity and can be applied in the bipartite as well as the multipartite scenario. 
Reviews of the entanglement versus separabil ity problem, from a vari ety of different 



viewpoints, can be found in Ref. 



37 



39, 



84 



85 



138, 



184 



263, 



283, 



284]. 



1.1 The EPR argument and Bell's inequalities 



The existence of entanglement 



2601 ] follows naturally from the quantum mechanical 



by Einstein, Podolsky 



formalism. This was first made explicit in the famous paper 
and Rosen (EPR), where it was used to argue that quantum mechanics as a physical 
theory is incomplete. Their argument runs as follows. Consider a particle with known 
position decaying into two equal particles. Without measurement, all we know is that 
the particles will drift apart with opposite momenta and that their centre of mass re- 
mains constant throughout. Assuming that both particles are well separated, there is 
no way a measurement on one of the particles can affect the other particle. This is 
the famous local realism assumption, which dictates that well separated systems can be 
completely and independently described. Now, measuring the momentum or position of 
the first particle enables one to predict either the momentum or position of the second 
system, without disturbing it. In quantum mechanics position and momentum are non- 
commuting observables and the theory cannot predict precise values for both, hence 
EPR are led to conclude that quantum mechanics is incomplete. A formidable philo- 
sophical debate followed on the foundations and interpretation of quantum mechanics, 



1.1 The EPR argument and Bell's inequalities 



3 



starting with Bohr's reply [27j]. For a historical account see Ref. 167]. 



In 1951 Bohm 26J] formulated a version of the EPR argument with an entangled 
spin system. Here both well separated parties (^4 and B) share one half of the singlet 
state 

liM = -^(|oi) - |io». 

Measurement of the spin of the first half in any direction reveals the spin of the other 
particle in that direction. By an EPR argument the quantum mechanical description of 
the second half cannot be complete. It is for such a system Bell [r| derived his famous 



66] 



inequality. We will present here a proof of the more general CHSH inequality 

Suppose the two parties share a large number of singlet states and measure them 
one by one. For each measurement, observer A chooses to measure the spin in direction 
a or a' at random (with results a n or a' n ). Analogously observer B measures the spin of 
his half in direction b or U (with results b n or b' n ). Let us denote with Q(a, b) the average 
value of the product of the measurements outcomes ab. This quantity is a measure of the 
correlation between the two measurement outcomes. Now the EPR argument implies 
that for each singlet pair, a n ,a' n ,b n , b' n have definite values (±1) before a measurement 
has taken place. By going through all combinations, it is easy to see that 

a n {K + b' n ) + a' n {b n - b' n ) = ±2. 

Because the left hand side is either -2 or +2, its average over all n is not larger than 2, 
thus we get on average 

\Q(a,b) +Q(a,V) + Q(a',b) -Q(a',V)\ < 2. 

This inequality gives an upper limit to the correlation coefficients. Choosing the unit 
vectors a, b, a' , b' in the same plane but separated by angles of vr/4, one can verify that 
this inequality is violated quantum mechanically for the singlet state (with a value of 
2\/2). Bell type inequalities can be derived for any number of parties, for more tha n tw o 



distinct measurement outcomes or for any number of alternative setups, see Ref. 



2251 ] . 



Violation of the Bell inequalities has been experimentally demonstrated using the 



polarisation of photons 



99J and experimental results are in excellent agreement 



with the predictions of quantum mechanics. Apart from entanglement in polarisation, 
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4 



violation of Bell inequalities wit h ph otons has been demonst rated using entanglement 



based on po sitio n and time 



momentum 



98, 



2881 ]. phase and momentum 243] and orbital angular 



1951 ] . Exp erim ents also hav e be en conducted without the use of photons, 



namely with protons [l80l ] and atoms [2481 ] . It should be mentioned however, that 
to date no experiment has been loophole free, so that technically, local realism has 
not been ruled out experimentally. An example of su ch a loophole is not having proper 



time-like separation of the detectors (see Ref. 



65, 



lOOl ] for a discussion of these matters) . 



1.2 The main road 



For any entangled pure state shared by an arbitrary 
that there exists a Bell inequality which is violated 



number of parties, it was shown 



104 



235]. For mixed states this 



does not seem the case. In this context it has been shown that although a s tate migh t 



not violate any Bell-like inequality, it still may be useful for teleportation 



154 



2331 ] 



Other sta tes only violate some inequality after a certain local generalised measurements 



1051 . 



2341 ] . Yet o ther states only show a violation when multiple copies are measured 



collectively 



2241 ] . It will be clear that Bell inequalities are a poor test of non-classical 



behaviour. It has proven convenient to characterise the set of entangled mixed states 
with a certain reasonable mathematical definition. The advantage of this approach 
is to have a workable definition of what is entangled and what not. Though only 
recently it has been rigorously proven that all states defined in such way behave truly 
non-classically (see Chapter 

We will give the definition of entanglement for bipartite states, generalisation to 
mul tipar tite states is straightforward. Let us first recall the definition of a pure state 



94 



260j. Let H be a Hilbert space such that H = H A <8> U B = C dl ® C d2 (with integers 



di-, d,2 > 2). A pure state \ip) e H shared by two parties A and B is called separable if 
and only if it can be written as 



IV) = Wa) ® IV 



B 



with \ij)A) £ TYa and \iJj r) g H-b, in other words, if \ip) belongs to the Cartesian product 
7~(-A x Hb C Ha ®Hb [2791 ] . Otherwise we call entangled. 



1.2 The main road 
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There is a simple way of determining whether a pure state is entangled based on 
the Schmidt decomposition: 



Theorem 1.2.1 (Schmidt decomposition 94 . \2l4 . 12591 ]). Let £ 7i, then there exist 
orthonormal states la,) and 16,;) such that 



^Aila*)^), (1.1) 



with \ > and Y^i \ = 1- 27ie coefficients A, are known as the Schmidt coefficients. 
The number of non-zero Schmidt coefficients is referred to as the Schmidt rank S(ip) of 
The state \ifj) is separable if and only if it has Schmidt rank one. 

Proof. Let 

with ^ = {ipjj \ a matrix of complex numbers. Recall that by the singular value 



decomposition 



1 261 ] every matrix ^ can be written as ^ = UDV, with U and V 



unitary and D a diagonal matrix with nonnegative elements. So we have 

ijk 

Now putting \a,j) = Uji\j), \bi) = Vik\k) and A, = da we obtain the form (|1.1|) . 
The orthonormality follows from U and V being unitary. To check whether a pure state 
is separable it is therefore sufficient to calculate the rank of □ 

Moving to mixed states, let us denote with j£(TC, TL) = Cy, = C-Ha^Hb — ^-'Ha ®^Hb 
the set of linear operators on TL. Since we are working over a finite dimensional space 
Cn = M.did.2 t ne space of complex matrices of order d\di- The subset of (Hermitian) 
positive semidefinite operators will be denoted The set of density operators on TL 
constitutes then the normalised operators in £i and we will denote them with Dn or 
simply T>. A bipartite mixed state p £ T> is called separable if and only if it can be 
written as a convex sum of pure product states, i.e. 



P = Y,PMt)W?\®\^fM 



Bi 
i I > 



with pi > 0, Yli Pi = 1) l^Pf) ^ and \if)f) G TLb- Otherwise the state is called 
entangled. The set of separable states will be denoted S and can thus be identified 



1.2 The main road 



6 



with the convex hull of the normalised operators in £^ (g) C^ b C £i ^ . If 

243] follows 



146 



such a decomposition exists, then from the Caratheodory theorem 
that we can always replace it with a decomposition of at most m 2 terms where m 
d\d2 = dimW. This definition w as fir st coined independently by Werner and Primas 
in 198sQ, but popularised in Ref. 



321]. The physical idea behind it is that a separable 
state can always be constructed by two independent parties with the aid of classical 
communication only (but again, a full physical justification came first later). 

1.2.1 Hahn-Banach theorem and entanglement witnesses 

To check whether a pure state is entangled it is sufficient to expand it in its Schmidt 
decomposition. For mixed states the problem is much harder as no such canonical 
decomposition can be obtained in a straightforward way. However, by exploiting the 
convexity of the set of separable states a useful characterisation can be obtained, which 
we will use to construct some necessary conditions for separability. 

The following theorem gives a geometrical characterization of the problem of deter- 
mining whether p £ D is contained in a certain convex subset C C D: 

Theorem 1.2.2. Let C C D be a convex set of states in a Hilbert space 7i, then for 
each p £ C there exists a Hermitian operator A S B(TC) such that 

Tr (Ap) < and Tr (Aa) > 

for all a G C . 

This t heorem is an immediate consequence of basic theorems in functional analysis 



181 



2471 ] . Namely, the Hahn-Banach theorem states that a convex set and a point lying 
outside it can be separated by a hyperplane W, and the Riesz-Frechet representation 
theorem t hen c haracterizes such hyperplanes. The hyperplanes W are commonly called 



witnesses 



2841 ] (they witness states outside C). 



Corollary 1.2.1 (Horodecki et al. [1321 ]). A state p S T> is separable if and only if 
Tr(pW) > for all Hermitian operators W such that Tr([\ipA){'^A\ <8> \^b)('4>b\]W) > 0, 
with \ipA) & 7~(-A an d IV'b) 6 Such operators W will be referred to as entanglement 
witnesses. 



1 R. Werner, LMS Workshop on quantum information and computation (York 2005). 



1.2 The main road 



7 



We have encountered an entanglement witness when we discussed the CHSH in- 
equality. From our discussion, it is easy to see that 



Wchsh = 21 - [a ■ a ® (6 + V) ■ a + a' ■ a ® (b - V) ■ a], 



(1.2) 



is an enta nglement 



matrices 



161 



163 



witn ess. Here a = (a x ,a y , a z ) is a vector containing the Pauli spin 



2841 ], Note also that a general entanglement witness W can be 



measured locally, by decomposing it as (non-convex) sum of product states. The prob- 
lem of finding an optimized decomposition both in the optimal number of projectors 
on product vectors and in t he o ptimal number of settings of the detectors has been 



partially solved in Ref. 



1071 . 



1081 ] . 



Let us now give two exa mples of entanglement witnesses, which will be important 



for the sequel. The first one 



1321 ] is the so-called swap or flip operator F, which acts 



on state vectors as F\iPa)\iPb) = \iPb)\iPa) or i n a basis of H: 

ij 

It is not so hard to check that this operator is indeed positive on product states. The 
swap operator is an entanglement witness because F has negative eigenvalues and will 
therefore necessarily act negative on some state (for instance for \ifi—})- The second 
example is given by 

W R = 1- dP + 

with P + the projector onto the maximally entangled state \ip+) = Since 
the maximum overlap of a separable state with P + is 1/d, and since W has negative 
eigenvalues, it is an entanglement witness (for instance for 

Intuitively the value min^ Ti(pW) can be associated with the distance of p to the 
set of separable states and in some way is related to the degree of entanglement of p. 
However, as such the minimum is ill-defined as the set of entanglement witnesses is 
unbounded. This can be overcome by imposing suitable constraints on the entangle- 



ment witnesses W, and it turns out that in t 
entanglement measures can be obtained [33J, 



lis w ay a variety of geometrical intuitive 



2301 ] . We will come back to the topic of 



geometrical entanglement measures briefly in Chapter [2j 



1.2 The main road 



8 



1.2.2 Jamiolkowski isomorphism and positive maps 

In this section we will reformulate the characterisation of separable states (Corol- 
lary [L2TT]) in terms of positive maps. This can be done by making use of the Jamiol- 
kowski isomorphism which establishes a one-to-one correspondence between linear op- 
erators and linear maps between operator spaces. We will outline the isomorphism 
and illustrate it for Hermitian preserving, positive and completely positive maps. We 
then apply it to the examples of last section, and obtain the reduction and the partial 
transpose criteria, respectively. The r esult s outlined in this section date back to the 
1970-80s' 53] but have only recently 1321 ] been used in quantum information. This 



section is meant as an introduction to the isomorphism and we will see more about 
positive maps and the context in which they were studied in Section [1.2.31 

A Hermitian map A : C-h a — * £-h b ^ s a ^ near map that maps Hermitian operators 
onto Hermitian operators. A positive map A : £t-Ca ~* is a linear map that maps 
positive operators p > onto positive operators A(p) > 0. A A;-positive map is a 
positive map such that the induced map 

A fc = t k <g> A : C £k C Ha -> C £k <g> £n B 

is positive for k > 1. A complete positive map (CP map) is a /c-positive map for all 
positive k. 

In the proof of the Schmidt decomposition we found it useful to interpret the coef- 
ficients ipij in the expansion of \t/j) = Y^ij ^ij'Ki) as the matrix {ipij} of an operator ^f. 
With this notation one can equally write 

|^,) = CS> 1)1-0+) = Vd(t ® * T )|^+). (1.3) 

This relation in effect establishes an explicit isomorphism between TCa <8> T~Cb and 



£(Ha,H r) 



295, 



296, 



251 ] and comes back frequently in quantum information (see for example 
340[]). Suppose now that we have a general Hermitian operator W £ C-Ha'&'Hb > 



using the spectral decomposition we can write 

i 

with all coefficients uii real. Using the above isomorphism we associate with every 



1.2 The main road 



9 



vector \(pi) the matrix <£j and write 

W = d^ Wi {l (8) $i)P+(l ® $}). 

i 

Defining a map A^y as acting as 

A w (A) = Y,Wi$iM> (1-4) 

i 

for all A, this can be written as 

W = d{l® A W )P+. (1-5) 

It is clear that when a linear map Aw can be written in the form of A1.4H it is automat- 
ically a Hermitian map. Conversely, using the above correspondence (|1.5|) . it is easy to 
see that any Hermitian map Aw will correspond to a Hermitian operator W and there- 
fore can be written in the form (11.40 . When W is a positive operator, all the coefficients 
Wi will be positive and Aw will become a completely positive operator. Conversely, it 
is again easy to see that any completely positive map can be written in the form (jl.4|) . 
Thes e results were first derived by de Pillis and Choi respectively (see also 



Ref. 



124 



131 



174 



178 



232 



276 



2961 ] ). we summarise them in the following theorem. 



Theorem 1.2.3. A linear map A : Ch a ~* ^H B is Hermitian (completely positive) if 
and only if for all operators A £ C-h a it can be written as 

A{A) =Y,hViAvl 

i 

for some operators Vi : ~^ ^~-'Hb> an< ^ rea ^ (positive) numbers Aj. 

In the proof of this theorem we made use of the relation W = d(l ® A\y)P+ which 
associates with each linear map Aw an operator W . Note that W defines uniquely the 
map Aw since a linear map is determined by its action on the basis Indeed, it 

is not so hard to see that (|1.5|) can be written as 

(k\A w (\i)(j\)\l) = (ik\W\jl). (1.6) 

Thus in this way we have defined an isomorphism between linear maps A : — ► Cy_ B 
and operators W £ Cn A ® Cn B ■ This is the famous Jamiolkowski isomorphism and 
can also be written as 

A W (X) = Ty a [W(X t 
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Linear maps A 
(A:C Ha ^£h b ) 



Hermitian maps 



Positive maps 



fc-Positive maps 



CP maps 



{C HA ®C n J 



(C Ha <g> C Hi 



Linear functionals on 



Linear functionals 
assuming real values 
on {C Ha ® C H]3 ) h 



Linear functionals 
acting positive 
on £+ a ® £+ b (or S) 



Linear functionals 
assuming positive 
values on Sk 



Linear functionals 

acting positive 
on (C Ha ®Ctl b ) + 



Figure 1.1: The Jamiolkowski isomorphism and correspondences for important classes 
of linear maps. The correspondence for /c-positive maps will be outlined in Chapter O 
Except for the addition of these maps this chart is from Ref. [531 ] . 
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Jamiolkowski 166] derived this isomorphism in order to characterise positive maps. 
As a positive map is Hermitian we can always write it in the form of (|1.4j) . However 
it is unclear how to distinguish between Hermitian and positive maps from this. The 
trick is to associate with the operator W the linear functional fw{\^)) = (i>\W\il)), this 
linear functional can be characterised as follows. 

Theorem 1.2.4. A linear map Aw : C-^ A — > £h b * s positive if and only if the corre- 
sponding functional fw on Ha <8> T~Lb is positive on product vectors. 

Proof. We can prove this theorem solely by examining (11.60 for the special case where 
K) = \j) = \ x ) and \ k ) = 10 = \y) : 

{y\A w (\x){x\)\y) = {xy\W\xy). 

If fw is positive on product vectors then Aw is positive on projectors, and therefore 
on arbitrary positive operators also. If Aw is a positive map, then from the above 
equation follows that fw is positive on product vectors. □ 

Thus in the language of the previous section, positive maps correspond to entangle- 
ment witnesses. Figure flTTI summarises the isomorphism and its relations for different 
classes of linear maps. With this we are in the pos ition of stating the main result of 



1321 ] (note that it can be extended to 



this section, which was first published in Ref. 

n 

multipartite systems, for details see Ref. [139|). 

Theorem 1.2.5. A density operator p acting on Ha is separable if and only if 

(1® A)(p) > 0, 
for all positive maps A : £h b ~^ ^Ha • 

Proof. With each linear map A : Cy_ B — ► Cy_ A there is associated an adjoint map 
A" 1 " : Cn A -> C-Hb defined by 

Tr[AA(B)] = Tr[A* (A)B] 

for all operators A, B. It is easy to verify that the adjoint map of a positive map is 
again positive. Suppose now that p is entangled, then from Corollary 11.2.11 follows 
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that there exists an entanglement witness W such that Tr(pW) < 0. By virtue of the 
Jamiolkowski isomorphism and Theorem 11.2.41 this is equivalent to 

Tr[(l ® A w )(P + )p] = Tr[(l ® A\ v )(p)P + ] < 

and since P + is positive it follows that (U ® A| 1/ )(/)) ^ 0. Conversely, let Aj^ be 
a positive map such that (II ® A\ v )(p) has a negative eigenvalue with corresponding 
eigenvector \<f>) = A ® 1| ?/>+)• With this we get 

(<f)\t ® A^(p)|0) = Tr[(A ® l)iy(^ < 0. 

Since W is positive on separable states, the same applies to (A® 1)W(A J< ® 1). This 
shows that p must be entangled. □ 

The proof of Theorem 11.2.51 reveals two important facts about the relation between 
entanglement witnesses and positive maps. Namely, given an entanglement witness W, 
the corresponding positive map A* w will always detect more states than the witness, as 
it detects all states detected by the class of witnesses (A<S> H.)W(A< (8)1) for all operators 
A. Furthermore, if a certain positive map A\y is negative on a state p, then the proof 
shows us how to construct an operator A such that the local transformation^ 

p^ (A* ® l)p(A ® 1), (1.7) 

will yield a state that can be detected by the entanglement witness W. 

In the previous section, we have seen two important examples of entanglement 
witnesses. The map associated with W = 1 — dP+ is given by A(A) = 11 Tr(A) — A. 
From this follows that all separable states p satisfy 

PA &> 1 - p > and 1 <g> ps - p > 0, 

while this is not necessarily the case for entangled states. It is not so hard to prove 



the positivity of this map directly. The associated 
reduction criterion and was first studied in Ref. 431 . 



sepa rability criterion is called the 



131]- 



Another entanglement witness was given by the swap operator F. It is not so hard 
to see that the associated map is given by the transposition. For a separable state 



2 See Section [3]T] for more on quantum operations. 
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p = YliPi^t')^^ ® \' l Pf){i ) f\ we have that (1 ® T)p > 0, while this is not necessarily 
the case for entangled states. We write (1 (g) T)p = p Ts and call this operation the 
partial transposition. When p Te > we say that the state has a positive partial 
transpose (PPT), otherwise the state has a negative parti al t ranspose (NPT). The 



223j | independently of the 



partial transposition criterion was first discovered by Peres 
notion of positive maps. This criterion turns out to be powerful, in particular, positivity 
of th e partial transposition is necessary and sufficient for separability when dimTC < 6 



1321 ] . It w as subsequently shown also t o be necessary and su fficie nt for states with low 



rank 



1511 ] , in particular for pure states 2841 ] , rank two states 1531 ] and rank three states 



1511 ] . In general though, it is not sufficient, and there exist PPT entangled states. We 
will see an example of such a PPT entangled state in the next section, where we also 
analyse the importance of the partial transposition in the structure of positive maps. 
We will denote the set of PPT states as V. The set V is, just like the set of separable 
states S, compact and S C V C V. 

It can be shown th at t he reduction criterion is strictly weaker than the partial 



transposition criterion 



1311 ] . except when dimWg = 2 or dim 7^ = 2 in which case 



both criteria are equivalent. This last statement follows from the identity 1 Tr(A) — A = 
(a y Aa y ) T for two dimensional matrices A. The reduction and the partial transposition 
criterion are important in the context of distillation, they will be studied in detail in 
Chapter [3] 



1.2.3 Some deeper results 

We start this section with the following classical theorem, which gives a complete 
characterisation of positive maps between low dimensional operators (for elementary 



proofs see Ref. 



84, 



183, 



2981]): 



Theorem 1.2.6 (St0rmer and Woronowicz 
map, with n < 3. Then A can be written as 



2771 



3281 ]). Let A : £ C 2 — ► £c n be a positive 



A = Ai + A 2 o T, 



with Ai , A2 completely positive maps. 
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We call a positive map decomposable if it can be written as a combination of 
a completely positive map and the combination of a completely positive map and 
transposition. It is easy to see that entangled states detected by a decomposable map 
will also be detected by the partial transposition alone. Thus from Theorem 1 1 . 2 . 5 1 we 
obtain that for all bipartite states p acting on TC with dimW < 6, positivity of the 
partial transposition is a necessary and sufficient condition for separability. 

The first example of an undecomposable map was given by Choi 
Ac : C c a — ► £,£3 be the linear map defined as 



50, 



51 



5J. Let 



A u + A 33 
-A 2 i 
-A 31 



-A 12 -A 13 
A 22 + A n -A 23 
—A 32 A 33 + A 22 



(1.8) 



for all 3 x 3 matrices A = [A{j]. He originally proved that this map is positive on 
real positive semidefinite (and thus symmetric) matrices , but it easy to see that this 
implies that the map is positive on all positive semidefinite matrices [52|. With a 
counterexa mple, we will show that this map is not decomposable. The exa mple i s du e 



to St0rmer [2781 ] but analysed in the context of entanglement in Ref. 
Define for < a < 5 

1 



75, 



138 



1491 ]. 



a 



7 



[2P + + aa + + (5 - a)a. 



with 



t7 + = i(|01><01| + |12)(12| + |20)<20|) 



and ct_ = Fo + F. This state has a positive partial transpose for 1 < a < 4, for 
other a the state is NPT entangled. Now it is easy to verify that (1L (g> Ac)(c) is not 
positive for 3 < a < 4 and 1 < a < 2 since it has an eigenvalue A = (3 — a)/2. This 
immediately shows that Ac is non-decomposable and that a is PPT entangled in these 
regions. For the region 2 < a < 3, the state can be shown to be separable by virtue of 
an explicit decomposition in terms of product states. Following Choi's example, many 



other non-decom posable maps have 



instance Ref. 



56 



111 



115 



179 



187 



jeen found and studied in the li terature, see for 



200, 



2171 . 



2461 . 



281 



282, 



289, 



3281 ] . 



The notion of decomposability can also be extended to entanglement witnesses via 
the Jamiolkowski isomorphism. Namely it is easy to see that any decomposable map 
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A\y corresponds to a witness 

w = p + q Tb , 

with P and Q positive semidefinite (and vice vers e). We will refer to such witnesses 



as decomposable witnesses. There exists 



1851 . 



18a] a canonical representation of inde- 



composable entanglement witnesses, which we will derive now. Central is the idea of a 
so-called edge state, which is a PPT entangled state 5 such that for all e > and all 
separable \ab), 

5 — e\ab)(ab\ 

is not positive or does not have a positive partial transpose. Edge states can be thought 
of as lying on the boundary between PPT and NPT states, and thus in this sense, they 
are the most entangled PPT states. From the definition it is clear that a state 5 is an 
edge state if and only if there is no product vector \ab) in the range of 5 such that \ab*) 
is in the range of 5 Tb . If such a product vector exists, we could subtract it from the 
state and still obtain a PPT state. Repeating this process will eventually yield an edge 
state (see also Figure [L2j) : 

Theorem 1.2.7. Every PPTES p can be decomposed as 

p = (l - p) Psep + p 8, 

with p sep a separable state and 5 an edge state. 

The following theorem provides a canonical construction of an entanglement wit- 
ness that detects a given edge state. Conversely it states that any non-decomposable 
entanglement witness can be written in a similar form. 

Theorem 1.2.8. (i) Let 8 be an edge state with P and Q the projector onto the kernel 
of 5 and S Tb respectively. Let 

W 5 = P + Q Tb , 



and 



e = mf(ab\W s \ab), 

\ab) 



then W\ = Ws — el is an non- decomposable entanglement witness which detects 5. 
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Figure 1.2: Schematic representation of the nesting S C V C T> and a PPT state as 
the convex combination of an edge and a separable state [39|] . 

(ii) Any non- decomposable witness W can be written as 

W = P + Q Tb - el 

with 

< e < inf(a6|P + Q Ts |a6). 

\ab) 

Here P,Q are positive semidefinite and satisfy Tr P5 = TtQ Tb 5 = for some edge 
state 6. 

Proof, (i) As the range of 5 contains no product vector \ab) such that the range of 
contains \ab*), there must exist such a product vector in their respective kernels, 
so that e > 0. By its very definition, W\ is positive on product states and since 
Tr(Vl^i5) = — e < it is an non-decomposable entanglement witness which detects 8. 

(ii) Let W(\) = W + All with A > 0. We denote with Ao the smallest A such 
that W(\o) = is decomposable. Since TV(Ao) is decomposable we can write it as 
W(X ) = P + Q Tb and thus W = P + Q Tb - X 1. Using the compactness of the 
set of PPT states V there will exist a state p such that Tr(pW(\o)) = 0. Now let 
p = (1 — p)p SC p + pd, with 5 an edge state, so that 

P Tr(W Psep ) + (1 -p) Tr(W5) = -A . 
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When p > this would lead to Tr(W(Xo)5) < 0, contrary to the assumption that 
W(X Q ) is decomposable. So p = 0, Tr(<5VF(A )) = and Tr(SW) = -A . It follows that 
Tr PS = TiQ Tb 6 = 0. The theorem follows by putting e = Ao, where the condition on 
e follows because W is an entanglement witness. □ 

This theorem can be generalise d st r aigh tforwardly to arbitrary convex sets and 



associated witnesses [39]. In Ref. 185l . 1 1861 ] it was shown how to optimise a non 



decomposable witness W (for instance constructed from an edge state). Optimising in 
this context means iterative subtraction of decomposable witnesses from W, yielding a 
stronger witness. We refer to the original references for details. 

With an entanglement witness W or alternatively the positive map Aw we can 
associate the following bihermitian positive form: 

FwixiiVi) = (xy\W\xy) = (x\A w (\y)(y\)\x) = ^W ij]k ix*y*x k yi. 

ijkl 

Suppose now that A\y has real coefficients in some basis, then for real Xi,yi, Fyy 
is a so-called real biquadratic form. Hilbert's 17th Problem is concerned with the 
question whether all real positive semidefinite forms can be written as the sum of squares 
(SOS) of forms (see Ref. 2451 ] for an excellent overview on Hilbert's 17th problem). 



Choi 



50 



52 



53|, 



2841 ] showed that non-decomposable maps with real coefficients allow 



to construct counterexamples for biquadratic forms. 

Theorem 1.2.9. (i) Let Fyy be a real biquadratic form which cannot be written as a 
SOS of linear forms, then the associated map Aw is non-decomposable, (ii) Conversely, 
if Aw is a non-decomposable with real coefficients, then Fw cannot be expressed as a 
SOS. 

Proof, (i) Note that decomposable maps acting on (real) symmetric matrices act as 
completely positive maps, so we can write 

A w (\ y )(y\) = J2Ai\y}(y\Af, 

i 

and with (x\Ai\y) = a\j we get 

Fw{xi,Vi) = {x\^w{\y){y\)\x) = ^2(x\Ai\y)(y\Af\x) = ^(a^Xiyj) 2 . 
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(ii) By reversing the argument of (i) it follows that if Fyy can be expressed as a SOS, 
then A\y is decomposable. □ 

Earlier in this section we have seen Choi's example of an non-decomposable map. 
The associated real form is given by 



Fc = XiVi + x 2 y 2 + x 2 y 2 - 2{xix 2 y 1 y 2 + x 2 x 3 y 2 y 3 + x^x^yi) + x x y 2 + x 2 y 3 + x^y x . 

By inspection Choi showed the positivity of this form and showed that it cannot be 
expressed as a SOS. This immediately proves the non-decomposability of the associated 
map. 

For a more abstract formulation of the separability problem, with various connec- 
tions to C*-algebra s the reader is referred to the series of papers by Majewski et al. 



196, 



1971 . 



19 



199]. 



1.2.4 PPT entangled states 

In this section we provide several ways of constructing PPT entangled states, with 
particular emphasis on the UPB construction and the range criterion. Based on this, 
we construct two new entangled PPT states with special ranks. For completeness, 
references will be listed to all other PPT entangled states kn own t o us . 



In the previous section we have seen St0rmer's example 



75, 



149 



278f | of a PPT state, 



constructed such that it is detected by Choi's indecomposable map. In general given an 



inde composable map, one can try to guess from it some state w 



lich is 



13a ] . see for example the PPT states constructed in Ref. 



18, 



36, 



'PT entangled 



2261 ]. Recently, it 



has been shown how to systematica lly construct PPT entangled states from a given 
indecomposable map 



112 



113 , 



114]]. Yet despite this advance, coming up with an 
indecomposable map is highly non-trivial and very much a matter of trial and error. 
Currently, the only systematic way of constructing PPT entangled states is based on 



unextendible product basis (UPB) [2; 



], Q. An UPB in a Hilbert space H = H a ®'Hb 



is a set of pure orthogonal product states which do not spa TL and such that there 
is no product vector orthogonal to all of them. This definition is readily extended 



3 Therefore, an unextendible product basis is in fact not a basis. 
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to multipartite systems. As an example, consider the following set of vectors 5 in 
n C 3 ® C 3 : 

l^ ) = -L|o)(|o)-|i)), | V , l) = _L(| )-|i))|2), 

|^ 2 ) = -L|2}(|1}-|2}) ; |^ ) = _L(|i)_|2))|0) ; 

|V 4 } = ^(|0) + |1} + |2})(|0} + |1) + |2}). 

Now if this set of vectors is extendible, then there has to exist a partition of this set of 
vectors 5", into two subsets S\ and 52, such that the local rank of the first subsystem 
for S\ is smaller than 3, and the local rank of the second subsystem for 5*2 is smaller 
than 3. In this case, one can construct a vector \ab) such that \a) is orthogonal to the 
first subsystem of S± and \b) is orthogonal to the second subsystem of S2 ■ But this is 
impossible since either 5i or 52 will contain at least three vectors, and it is easy to see 
that any three vectors, seen by either party span the full space of that party. We refer 



to the original papers and the review 71|] for a plethora of other examples. 

From any UPB one can construct a bound entangled state as in the next theorem. 

Theorem 1.2.10. Let the set of n vectors {\ipi)} be a UPB in TL = TLa ®T~Lb , with 
dimW = d, then 

1 - 

is a PPT entangled state. 

Proof. The state p has positive partial transposition as the set {|V ; j)(V'j| Ts } is again a 
UPB. It is easy to see that p is entangled and in fact an edge state as, by definition, 
there is no product state in the space orthogonal to the UPB. □ 

In Ref. 2851 ] it was shown how to systematically construct entanglement witnesses 
detecting UPB states. The generalisation of this construction for general edge states is 
precisely the one we outlined in the previous section (Theorem II. 2. 8h . Based on these 
witnesses, PPT ent ang l ed s tates in the neighbourhood of UPB PPT entangled states 



can be constructed 



15 



228 ] 



As we have seen, an edge state p contains no product vector \aibi) in its range such 
that |aj6*) is in the range of p Te . This condition is clearly a sufficient fo r a s tate to be 



PPT entangled. We can make this criterion slightly stronger as follows 



1461 ]: 
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Theorem 1.2.11 (Range criterion). If a state p is separable there must exist a set of 
product vectors {\aibi)} spanning the range of p such that {\aib*)} span the range of 



T 

Pb- 



Thus edge states violate the range cri terio n in the strongest way. The range criterion 
was originally devised by Woronowicz 3281 ] in the context of indecomposable maps. 



This criterion is remarkable simple and powerful, in particular for proving the PPT 
property of spar se matrices . Examples of st ates violating the range criterion have been 



published in Ref. 



38, 



53, 



appeared over the years 



97jjl4e 



11 



148, 



79 



12 



3281 1 . Various other PPT entang l ed st ates have 



129 



152 



1651 . 



2271 . 



3101 . 



322, 



3321 ] . In most 



cases the states are shown to be inseparable by showing that they can help performing 
some non-classical tasks. In Section l4~4"l we will say more about these so-called activation 
processes. 

We conclude this section with the presentation of two new PPT states 62j. The 



new types are (5, 5) and (6, 6) edge states. Here in (iV, M), N denotes the rank of the 
state, while M denotes the rank of p Te . The question of their existence popped up in 
the analysis of a conjecture on 2-positive maps between 3x3 matrices. We come back 
to this conjecture in Section [231 



A (5,5) edge PPTES 

Consider the following (5, 5) state: 



£(5,5) 
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and its partial transpose: 
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It is not so hard to verify that both operators are positive semi-definite and have rank 5. 
An analytical expression of their eigenvectors and eigenvalues is however quite complex. 
To show that pr^g) is an edge state we will show that it violates stro ngly the range 
criterion. For this we will use the 'divide and conquer technique' from [l4fj]. 

From the matrix representation above we see that every vector in the range of ^(5,5) 
can be written in the form 

V = (0, A, -E - F, C, 0, D, D, E, F), A, C, D,E,F € C. 
Now we have look at those vectors which can be written as a product 

V = (s, t, v) ® (x, y, z) = (sx, sy, sz, tx, ty, tz, vx, vy, vz). 

Taking these two conditions together we can therefore characterise all product vectors 

in the range of P(5 5)- 

From the condition sx = we can distinguish the following sub cases: 

1. x = 0, s 7^ 0, we have vx = D = = tz and therefore either t = or z = 0. 

Without loss of generality we can also put s = 1. 

1.1. t = 0, as v(y + z) = —z we have v = —z/(y + z). When y = —z, then 

z = = —y = x and we obtain the null vector. Thus the only case that remains is 



V = s 1,0 



y + z 



with y / —z. 
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1.2. z = 0, from £y = follows that £ = (otherwise y = 0) thus this case is already 
covered in 1.1. 

2. s = 0, from E = —F follows that vy = —vz. From ty = follows that either 
£ = or y = 0. 

2.1 £ = 0, since D = = -us we have that x = and jj^O, and thus z = — y. Thus 
we get 

V=(0,0,v)®(0,y,-y) 

2.2 y = 0, and thus vy = —vz = so either z = or v = 0. 

2.2.1 z = 0, and thus vx = £2 = 0. Because i/0we have that v = 0: 

V = (0,£,0) (g> (a:, 0,0) 

2.2.2 v = and thus vx = tz = 0. Because t / we have that z = and we arrive 
at 2.2.1. 

We now do the same analysis for P^ 5 y It is easy to see that every vector in its 
range can be written as 

V = (0, A, B, C, 0, D,E,A + 2B, —A — 2B + C + D + E), 

with A, B,C, D, E G C. We need to take the intersection of these vectors with those 
vectors which can be written as a product 

V = (s, £, v) ® (x, y, z) = (sx, sy, sz, tx, ty, tz, vx, vy, vz). 

From the condition sx = we can distinguish the following sub cases: 
1. x = 0, and since ty = we have that y = or t = 

1.1. y = and thus v y = and = — 2sz = 0. Since z = gives the null vector 
we get s = 0. Now = — sy — 2sz + tx + £z + t> x = tz and thus t = v, and thus 

F = (0,£,£) <g> (0,0, z) 

1.2. £ = 0, and since = — sy — 2sz = — vy either v = or 2 = — y 

1.2.1. v = and thus sy = — 2sz or since s = is trivial we get y = — 2z and thus 

V = (s,0,0) <g> (0,-2z,z) 
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1.2.2. z = —y, and thus — vy = —sy + 2sy = sy. Since y = —z = is trivial we get 
v = — s and 

V=(-s,0,s)®(0,y,-y) 

2. s = and thus vy = 0, we also have ty = thus either y = 0orv = t = 0, but 
this last case gives us again the trivial vector. From vz = —sy — 2sz + tx + ty + tz + vx = 
tx + tz + vx we have (v — t)z = (v + What remains is 

V = (0,v,t) ® (x,0,z), 

with either z = f^fx or x = (but only one of those when v = ±i). 

It is now straightforward to check that when v = (s, t, v) (x, y, z) belongs to the 
range of P(5,5) that v' = (s,t,v) (x*,y*,z*) does not belong to the range of pT§ 5 y 
This concludes the proof that P(5 ; 5) is a (5, 5) edge state. 

A (6, 6) edge PPTES 

Consider the following state 

1 0000000 -1 

20 -1 0000 

1 1 

-1 010000 1 

1 1 

000010 -1 

1 1 2 

0000 -1 01 

-1 0010000 3 



1 

P(6,6) - ^ 
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with partial transpose 
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It can be checked easily that both matrices are positive definite and have rank 6. 
Following a similar strategy as in the previous section we now proceed to show that 
jO(6,6) is an edge state. Every vector in the range of P(6,6) can be written in the form 

V = (A, B, C, D, E,F,C + E, —F, B + 2D- A), (1.10) 

rp 

with A, B, C, D, E, F £ C, whilst every vector in the range of takes the form 

V = (A, B, C, D, -A, E,E-C, D, F), (1.11) 

with A, B, C, D,E,F € C. The 'divide and conquer' method does not work so well here 
as we have no zeros to start with. Instead we will show directly that if the product 
vector 

V = (s, t, v) (g> (x, y, z) = (sx, sy, sz, tx, ty, tz, vx, vy, vz). 

belongs to the range of P(6,6)> then V' = (s,t,v) (8) (x*,y*,z*) does not belong to the 
range of p^f 6 ) • Comparing equations (jl.lUp and (|1.11|) with their product form we get 
the following set of equations: 

vx = sz + ty 
vy = —tz 

vz = sy + 2tx — sx 
ty* = —sx* 
vx = tz — sz 



1.2 The main road 



25 



We consider the following sub case 

1. vy 0, and it follows that all parameters are different from zero. Without loss 
of generality we can put x = 1 and get 

v = sz + ty 
vy = —tz 
vz = sy + 2t — s 
ty* = -s 

v = tz — sz 
vy = t 

We now use t = vy* and s = —v(y*) 2 to eliminate t and s, after which we can divide 
every equation by v resulting in 

1 = -(y*) 2 z + y*y 
y = -y*z 

z = -(y*) 2 y + 2y* + (y*f 
1 = y*z* + {y*) 2 z* 

From the first and the second equation we get yy* = ^, 
fourth equation we get 1 = y*(z* — y) or y = yy*(z* —y) = 
From the third equation we get 

3y* = -(y*) 2 y + 2y* + (y*) 2 

or 

3 = -y*y + 2 + y* 

and thus y* = | which contradicts yy* = \. 

2. v = y = 0, in this case it is straightforward to see that x = y = z = 0. 

3. y = and v / 0. From sx* = and tx* = follows that s = t = or x = 0. In 
the first case we get vx = vz = or z = x = by the assumption. In the second case 
follows from vz = that z = 0. 



while from the second and 
- y) and thus z = 3y* . 
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4. v = and from tz = and tx* = follows that t = 0orx = z = 

4.1 t = 0, from sir* = sz* = follows that either s = t = v = 0orx = z = = y, 
but y ^ by assumption. 

4.2 x = z = 0, from ty = sy = follows that s = i = v = 0. 
This concludes the proof. 

1.2.5 A numerical approach 

In this s ectio n we give an example of a numerical approach to the separability problem. 



Gurvits lldl | reformulated the separability problem as a weak membership problem 
(allowing some error in the decision) and showed that even in this scenario the problem 
is NP-hard. Yet, as demonstrated by the partial transpose criterion, this does not 
exclude the existence of efficient algorithms for low dimensional systems. There have 
been several algorithmic proposals, which can be categorised in three classes: 

From the outside Here a hierarchy of tests is devised which in every step detects 
entanglement of some states and provide a corresponding witness. Typically the 
first tests will detect highly entangled states, while further tests will detect more 
weakly entangled states. Ideally, the hierarchy should be complete, i.e. in the 
asymptotic regime every entangled^ state should be detected. Examples of such 
algorithms can be found in Ref. 



32, 



35,y,y, 3- 



From the inside In this approach, again a hierarchy of tests is devised, but this time 
able to deliver a certificate for separability. States detected in the first tests will 
be typically close to the maximally mixed state, while further tests will be able 
to detect states closer to the boundary of entangled states. The logic behind this 
is that for weakly (classical) correlated states there is more freedom in finding a 
convex decompositio n in t erms of product vectors. Examples of such algorithms 
can be found in Ref. 



159, 



325] 



Distance measure The starting point here is to take an entanglement measure E{p) 
which satisfies E(p) = for all separable states and E(p) > for all entan- 
gled states. Examples of such measures include the entanglement of f orm ation 



211 ]. and geometric measures such as the best separable approximation 1821 ] (see 
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Section [3.31 for an introduction to entanglement measures). The algorithm typ- 
ically w orks then by calculating E(p) to a certain accuracy. For examples see 



Ref. 



93 



1821 ] . 



An exhaustive review on algorithms for the separability problem, with particular em- 
phasis on comparing the different complexities of the algorithms, can be found in 
Ref. [162 ] . 



The most successful algorithms as of today use of optimisation theory, and in partic- 
ular semidefinite programming. Semidefinite programs (SDPs) are convex optimisation 
problems which can be written as the minimisation of a linear object_ function, subject 
to semidefinite constraints in the form of linear matrix inequalities 



3J: 



... t 

minimise c x, 



subject to Fq + xiFi > 0, 



where {x{\ is a finite set of real variables, c a constant vector and Fi some fixed Her- 
mitian matrices. Semidefinite programs have two appealing features. Firstly, efficient 
algorithms are available for solving SDPs in polynomial time with arbitrary accuracy. 
Secondly, with each (primal) SDP there is associated a dual SDP: 

maximise — Tr [Fq Z] , 
subject to Tr[FjZ] = q, Z > 

Here Z > means that Z is a positive semi-definite (Hermitian) operator. Now it is 
not so hard to see that c T x is always larger than — TrfFoZ], so that solving the dual 
problem gives a lower bound on the primal problem. It is often the case that these 
values coincide, in which case the SDP's have the so-called strong duality property. 
Suppose now that c = 0, then the primal problem becomes a feasibility problem. Then 
a positive value of — Tr[i<o^] shows that the primal problem is not feasible. In other 
words the operator Z gives a certificate for the infeasibility. Needless to say that such 
certificates are of great value in decision problems such as the separability problem. 

The two most important algorithms that use SDP, as of today, are the algorithms of 
Eisert et. al. [3] and Doherty et. al. [^J. The first approach is based on calculating the 
distance of a state to the set of separable states. In first instance this leads to a global 
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minimisation problem with polynomial constraints. Using Lasserre's method of semi- 
definite relaxations, this can then be reformulated as a sequence of SDPs. The hierarchy 
of tests obtained is complete, as every entangled state will be found at some stage in 
the hierarchy. Since a global optimum will be found at some point also separability can 
be detected 89]. 



Doherty's approach is based on the existence of symmetric extensions for separable 
states. It works from the outside and is also complete; furthermore, when a state is 
found to be entangled, the dual problem yields a corresponding entanglement witness. 
These entanglement witnesses turn out to have interesting algebraic properties, which 
is why we will say a bit more about (without proofs though). 

Suppose a state p on Ha <8> Ti-B is separable, then we know that it can be written 
as P = 'YliVi\^ ) i){ ( t > i \ ® l&K'&l- With p, we can associate the following state: 

P = /. 1 Pi |&)(<fti| ® ® \<l>i}{<Pi\ 

i 

on 7i.A<37~(-B<£>7~(-A- This extension is naturally symmetric under interchange of the first 
and third party. It also has positive partial transposition (with respect to any subsys- 
tem) and by tracing out the third party, we recover the original state. An extension for 
an arbitrary state satisfying these three properties will be called a symmetric exten- 
sion. Equally, we can consider a symmetric extension with 2 copies of the first party. 
Using de Quantum de Finetti theorem one can prove the following characterisation of 
separable states: 

Theorem 1.2.12. A state p on TLa®7~Lb is separable if and only if it has a symmetric 
extension to any number of copies of subsystem A. 



Apparently this theorem has been around for some time, see Ref. [75[ for its history. 
The new result by Doherty et. al. was to show that the feasibility problem of finding 
a symmetric extension to a fixed number of k copies can be cast into a semidefinite 
program. In this way, one obtains a hierarchy of tests, in which each test is stronger 
than the other. Indeed, when an extension to k + 1 copies exist, then naturally an 
extension to k copies exists. The first test is the usual partial transpose criterion, while 
further tests go beyond that and are able to detect PPT entangled states. In practice, 
it turns out that the algorithm detects all known PPT entangled states in the second 
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test of the hierarchy. The only major disadvantage is, that in practice, it is limited 
to density matrices acting on Ti with dimTi < 36 because of memory constraints of 
existing implementations. 

If the primal problem fails, the dual problem yields a corresponding entanglement 
witness. In Theorem 11.2.91 we have seen that when a positive map A\y with real 
coefficients is decomposable, then the associated real biquadratic form can be written 
as a SOS of linear forms. This result readily generalises to arbitrary positive maps. 
Now remarkably we have the following 751 ]: 



Theorem 1.2.13. Let E\y(x,y) be the bihermitian form associated with an entangle- 
ment witness W obtained from the (k + l)th test of the hierarchy. Then F(x,y) = 
(x\x) k Ey/{x, y) can be written as a SOS. 



751 ] the following characterisation of entanglement wit- 



With this one can prove 
nesses. 

Corollary 1.2.2. Let W be an entanglement witness in the interior of the set of en- 
tanglement witnesses (that is Ti(p sep W) > for all separable p se p)- Then there is a 
k > such that 

F(x,y) = (x\x) k E w (x,y) 

is a SOS. 

Note that a similar characterisation can be deduced from it for positive maps. As an 
example consider the St0rmer state a given by (|1.9p . detected by Choi's map Ac (|1.8p . 
The second test of the hierarchy yields just the entanglement witness Wc associated 
with Ac* and provides the SOS form for it: 

A{xy\W c \xy)(x\x) = \2x XQy - 2x 1 x^y 1 + x^yo - x 2 x* y 2 \ 2 + 3\x2xly 2 - x 1 x* 1 y 1 \ 2 + 
3\y xixl - x 2 y2Xo\ 2 + 3\x 2 x 2 y 2 - x xly \ 2 + ftxox^ - x 2 x yl - xix 2 yo| 2 + 
3\x x 2 yl - xix 2 y* \ 2 + \2x Q x* y 2 - 2x\_x\y\ + x 2 x\y 2 - XQX* 2 y^+ 
\2xQx\yv - 2x 2 x 2 yi + x 2 x*±y 2 - x^lyil 2 > 0. 
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1.3 Permutation criteria 



1.3.1 From cross norm over realignment to permutation criteria 

In this section we will introduce a class of separability criteria, known as the permu- 
tation criteria. We will introduce them in the way they have historically been found. 
Our star ting point is a characterisation of separable states known as the cross norm 



criterion 

as 

where \\A 
Rudolph 



249f| . The greatest cross norm || • || 7 is defined on operators in p € Cha ® £-Hi 



IpIL = inf 

Ui,Vi 



Ui <g> Vi 



TiKAA^) 1 / 2 ] is the trace norm (sum of the singular values). With this 



249j] proved the following theorem. 



Theorem 1.3.1. Let p be a density operator on TCa <8> TCb- Then p is separable if 
||p|L = 1 and entangled if ||p|| 7 > 1. 

The proof is technical and we omit it here. Note that the di ffere nce ||p|L — 1 satisfies 



2501 ] . Unfortunately the 



all the usual basic requirements for an entanglement measure 
greatest cross norm is in practice, with the exception of highly symmetrical states, very 
hard to calculate. 

One can equally well define other cross norms, by using a different norm in the 
definition. Let us denote by || • || s the cross norm induced by the Hilbert-Schmidt norm 
||yl|| 2 = [Tr^A 1 ")] 1 / 2 . From ||p|| s < \\p\\^n follows that a state p with \\p\\ g > 1 is 
e ntangled. It turns out that ||p|| s is easily computable, the construction is the following 



251 



254j . see also [5 



In Section [1.2.21 we have seen an isomorphism between vectors \ip) 6 TCa <8> TLb an d 
operators ^ E C(TLa,1~Lb) defined by 

The space C(TCa,T~(-b) is a Hilbert space, when furnished with the Hilbert-Schmidt 
inner product Tr^^) for \I', $ G £{TLa,'Hb)- It follows that the isomorphism is a 



This is because || • || 7 majorizes any unitarily invariant norm (hence the name greatest cross norm). 
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Hilbert-space isomorphism, i.e. isometric as (ip\<f>) = Tr(*& "!"<!>). Let TL = Ha ®Hb and 
p G CiTL). We can decompose p arbitrarily as 



with Ui G C{TLa) and G £(TLb)- Now we apply the above isomorphism to both it, 
and and obtain 

\p) = ^2 1"*) ® 

i 

with 6H® "H. Applying the isomorphism on |p) we obtain 

i 

The operator R(p) G is independent of the chosen decomposition of p and in 

general p ^ R(p)- This is because we can interpret \p) = |p)i234 G H® H as an element 
of (C d ) m = Hi 7i 2 ® H.3 H4, then with p acting on W x H 3 we have 

d 

13|24> 

ij=l 

as opposed to 

12|34- 

Here |VM-)i3|24 denotes the maximally entangled state with respect to the 13| 24 division 
(and analogue for |^+)i2|34)-_Th e next lemma is a well known identity from the theory 
of operator norms (see e.g. 2581 ]) 



Lemma 1.3.1. For any operator R G C(H) we have 



\R\\ = inf 

Ui,Vi 



^2 y/\{ui\ui)\ ■ \f\{vi\vi)\ 



With this lemma, we can prove that ||i?(p)|| = \\p\\ g and from that we obtain an 
easily computable necessary criterion for separability: 

Theorem 1.3.2 (Realignment criterion). For any state p G C(7i) the identity 

\\R(P)\\ = \\p\\g 



holds. If \\R(p)\\ > 1, then p is entangled. 
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Proof. Let p = ^ U{ (g) fj, then ^/| | = yTr(uJuj) = ||?ij||2 (and similar for Vi). 

The theorem then follows from the previous lemma together with the definition of || • || 9 
and inequality \\p\\ g < ||p||7 < 1 for separable states p. □ 

From the discussion above it follows that R{p) is nothing more then a rearrangement 
of matrix ent ries of p. For this reason, it is often referred to as the realignment criterion 



44 



2181 . 



252]. Some analytical properties of the realignment criterion have been studied 



m It turns out that the realignment criterion is independent of the partial 

transposition and hence can detect PPT entangled states. From our discussion it is 
easy to see that the realignment criterion corresponds to the operation R which acts 
on basis states as 

i^xv-ii ® * ® |0?>mi- 

Due to linearity, this operation is well defined for arbitrary quantum states. The partial 
transposition is also based on reordering of matrix entries. In the same notation, we 
denote it as the operation T which acts on the second subsystem as 

A state p is entangled if the trace norm ||T(p)|| > 1. Note that the usual partial 
transpose criterion says that a state is entangled if T(p) has some negative eigenvalues, 
but as T(p) is Hermitian this is equivalent to saying that ||T(p)|| > 1. This was first 
observed in Ref. [3051 ]. where it was shown that the quantity ||T(p)|| gives rise to a good 
entanglement measure (see also \^5\). 

This striking similarity between the two cr iteria led to the consideration of general 
permutations of multipartite density matrices 1441 ] . Consider an r-party state belong- 



ing to a Hilbert space 7i, whose subsystems have the same dimension d. A general 
state p E Ji. can be expanded as 

P= ^2 Phi2,i3i4, •• .,i2r-i»2r 1*1*3 •■■ *2r-l)(«2«4 •■ -hr\, 
il,i2,—,i2r 

where all indices run from 1 to d. Let S^r denote the symmetric group with (2r)! 
elements, that is, the group of the permutations of the set {1, 2, . . . , 2r}. We define for 
each permutation a G S2 r a corresponding map A CT on operators acting on TL as 

[^(P))i 1 i 2 ,...,i 2 r-li2r = ^Ml)M2)v,M2r-l)M2r) ' C 1 " 12 ) 
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We will represent permutations as [cr(l) u(2) . . . o~(2r)] or in disjoint cycles. With this 
notation the partial transpose criterion corresponds to the permutation [12 43] = (3, 4) 
while the realignment criterion corresponds to the permutation [13 2 4] = (2, 3). 



1441 ]). Let p be a r -party state 



Theorem 1.3.3 (Permutation criteria for separability 
acting on a Hilbert space p and a £ S2 r an arbitrary permutation. If 

IIMp)|| > i. 

with || • || the trace norm, then p is entangled (that is, not fully separable). 
Proof. Let p be a fully separable state on Tt, that is 

p = J^pM) W4l ® Iv4>(v4l ® • • • Ivi-i>(v4rl = J^PiPi, 

i i 

with V^fc-i = V^fc an d Pi the corresponding pure states. Then A a (p) corresponds to 

K( P ) = E^K(i))(C|)l ® ICf3))(C?4)l ® • • ■ ICf 2f -i))(C? 2r )l = E^(ft)> 

with = ip'jf if cr(i) and i have different parity and if} 1 ^ = t/ji otherwise. We need to 
show that IIA^p)!) < 1, but because of \\A + B\\ < \\A\\ + ||B[| it is sufficient to show 
that ||A CT ( / 9 i )|| < 1. But this follows easily from ||A(8)B|| < ||A||||5|| and || \(f>)(ip\ \\ = 1 
for normalised vectors \<j>) and \ip). □ 

As a corollary we will show here how to construct an entanglement witness for a 
r-party state p given that ||A CT (/o)|| > 1 for some permutation a. For th is we gen e ralis e 



the bipartite construction of Chen and Wu in Ref . [55f] , see also Ref . [KM Il45l . 13321 ] . 
Here and elsewhere permutations are evaluated from left to right. 

Theorem 1.3.4. Let ||A CT (p)|| > 1 for some permutation a, and define 

a = to~~ 1 t 

with t the global transposition and a^ 1 the inverse permutation of a. Then 

W = I - K{VU ] ) 

is an entanglement witness with Tr(pW) = 1 — ||A CT (p)|| < 0. Here U and V are the 
unitary matrices corresponding to the singular value decomposition of A a (p) = UDV^ . 
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Proof. It is straightforward to check that Tr(pW) = 1 — ||A CT (p)|| from the identity 

Tr[AA a (B)\ = Tr[(A ff -i {A T )) T B\ = Tr[A a (A)B], 

for arbitrary operators A and B. This identity can be verified by explicitly writing the 
trace as a sum over matrix indices. 

1 use 



127t], 



Now all that is left is to show that W is positive on separable states. We wi 
of the following variational characterisation of the trace norm (see for instance 
pl95): 

||A|| = max{| Tr(XAY)\ \XX ] = YY ] = 1}, (1.13) 

for arbitrary A. Let p sep be an arbitrarily separable state, then 

\\A a (p 8ep )\\=Tr(D')<l, 

with Ao-(p sep ) = U'D'V' the singular value decomposition of A -(p scp ). From (|1.13p 
follows that Tr[X A a (p sep )Y] is maximal when X = U' and Y = V and thus 

Tr(W Pscp ) = 1 - Tr(A a (Vtf) Psep ) 
= 1 - Tr(C/tA CT ( /9sep )y) 



> 1 - TiiU'A^p^V 



= 1-11^(^)11 >0. 

□ 

For an r-party state, there are (2r)! permutation criteria, but not all of them are 
inequivalent. For instance for bipartite systems it is easy to see that partial transpo- 
sition with respect to the first system is equivalent to the partial transposition with 
respect to the second. For bipartite systems all permutation criteria turn out to be 
equivalent to either the partial transposition or the realignment criterion 



144| ] . For 

more than two parties the classification of inequivalent permutation criteria is an ope n 



problem. First steps toward such a classification have been made in Ref. 



For three parties it was implied that there a re 9 inequivalent criteria 
four parties at most 34 inequivalent criteria [3241 ] . 



96 



96, 



144 



3241 ] . 



3241 ] . and for 
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In the next section we define exactly what we mean by independent criteria and 
review the work of Wocjan and Horodecki 3241 ] on so-called combinatorially indepen- 



dent permutation crit eria and introduce their graphical notation. Then we point out 
that both papers 



96, 



324] overlooked an obvious fact. We show that for three parties 
there are just 6 independent criteria and 22 for four parties. We prove that the class of 
the combinatorially independent permutation criteria contains some equivalent criteria 
which always occur in pairs. We completely classify these criteria and give a new upper 
bound of the number of independent criteria (this is Theorem II .3.7j) . It will become 
clear that most likely there are no more equivalences in the remaining criteria, although 
we were unable to provide a complete proof of this. A proof is presented for 2, 3 and 
4 parties and numerical evidence for up to 8 parties is given. These results were first 
published in Ref. [63l |. 



1.3.2 Definitions, graphical notation and first results 

Let a and /i be two permutations in Sir- We call the corresponding entanglement 
criteria a and \i dependent if and only if 

||A CT (p)|| = ||A>)||, (1.14) 

for all quantum states. Else, the permutation criteria are called independent. Observe 
that if two permutation criteria a and \i are independent by the above definition, it 
could still be that they detect exactly the same entangled states. This would occur if 
we had 

IIA^)!! > 1 & ||A„(p)|| > 1 

for all quantum states p. However, as we will see in the next section, the condition in 
(jl,14p is already good enough to reduce significantly the number of permutation sepa- 
rability criteria. Furthermore, we prove later that the criteria which are independent 
are also truly independent (in the sense that no criterion is strictly stronger than any 
other) for up to four subsystems. 

Our necessary conditi on e xtends the necessary condition of combinatorial indepen- 
dence introduced in Ref. 3241 ] . The definition of combinatorial independence is similar 



to the definition of independence. The decisive difference is that two permutations are 
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combinatorially dependent if and only if equality in (jl,14|) is achieved for all operators 
and not only quantum states. In this case, it is easy to see that the maps A CT and are 
related by a norm-preserving map A. That is A a (A) = A(A fl (A)) with ||j4|| = ||A(vl)|| 
for all operators A. Moreover, it can be shown that if A a and A^ are related via a 
norm-preserving map A, then A = A u must come from some permutation v. We call 
such a permutation a norm-preserving permutation. 

An example of such a norm-preserving permutation is the global quantum transpose 
(GQT), which transposes the complete system. It can be written as 

r = (l,2)(3,4)---(2r-l,2r). 

Another example of a norm preserving map A is the unitary transformation A(p) = 
UpV, where U and V are unitary operators. Reordering the different parties in the 
density matrix representation is an example of such a unitary transformation. Consider 
for instance the transformation 

p = y,p? ® pf ->p' ' = XX ® p^ 

i i 

This mapping can be imple men ted by means of multiplication on the left and on the 



right by the swap operator 



3241 ] 



To illustrate the necessary condition of combinatorial independence let us consider 
again a bipartite system. The operation R' induced by the permutation (1, 4) gives rise 
to the same criterion as R induced by the permutation (2, 3). Indeed, the permutation 
R't transforms [12 34] into [2413], which up to reordering of the parties is equivalent 
to [13 2 4]. This is just the transformation defined by R. 



3241 ]). The group T of norm 



Theorem 1.3.5 (Combinatorially independent criteria 

preserving permutations can be generated as 

T = ((2k, 21), (2k -1,21- l),r), 

where 1 < k, I < r and r as before denotes the GQT. The combinatorially independent 
permutation criteria correspond to the right cosets Sir IT of T in Sir- The number of 
independent permutation criteria is therefore not larger than | ( 2 r r ) . In this number, the 
class of trivial norm preserving criteria is also counted. 
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graphical representation 


corresponding permutation 


name 


k • — ► • I 


(2k, 21 - 1) 


reshuffle Rki 


k • < — • I 


(2k- 1, 21) 


reshuffle R' lk 


k 


(2k - 1, 2k) 


partial transpose 


k • 





identity 



Figure 1.3: Basic permutations, k < I (Table from Ref. 3241 ]) 



In the same paper, the authors also devised a graphical notion of the criteria, which 
leads to a way of selecting a simple representative for the right cosets. They decompose 
any permutation as a combination of 4 elementary permutations: the partial transpose, 
two types of realignment or reshuffling, and the identity. The corresponding graphical 
notations are loops from a subsystem to itself, arrows from one subsystem to another 
and free subsystems (no loops or arrows), as graphically depicted in Figure [L3l We 
call k the head and I the tail if there is an arrow from k to I. If there is a loop in m, 
then m is both head and tail of the loop. We call the support of an arrow or a loop the 
set containing its head and tail. A configuration of arrows and loops is called disjoint, 
if the supports of all pairs and loops are disjoint. 

The orem 1.3.6 (Canonical representation of combinatorially independent criteria 



3241 ]). The right cosets S2 r /T can always be represented by a disjoint configuration 
of arrows, that is, they can be decomposed as a combination of the 4 elementary permu- 
tations. It is always possible to represent a random permutation a £ S^r by a disjoint 
configuration by applying the following consecutive rules 

Rule 1: Pruning Writing a as a product of disjoint cycles, we can prune all cycles 
of a such that there are no adjacent even or odd numbers. 

Rule 2: Chopping We can chop up a into disjoint transpositions. 

Rule 3: Exchanging heads By exchanging heads we then obtain a disjoint configu- 
ration of arrows. 

Rule 4: GQT By reversing the direction of all arrows and replacing loops by free 
subsystems and vice versa we obtain an equivalent disjoint permutation. 



1.3 Permutation criteria 



38 



All criteria that can be represented by a disjoint configuration of arrows are combi- 
natorially independent up to a application of Rule 4, o,nd possible application of Rule 
3. 



3241 ] . 



This theorem differs slightly from the corresponding one in the original paper 
They claimed that any two permutations represented by disjoint configurations are 
independent, up to application of Rule 4. As the following example illustrates, this 
is not always the case: by exchanging heads of [13 245 768] we obtain [17645 32 8], 
which can be represented by a different disjoint configuration. And both configurations 
are not related via the GQT. 

1.3.3 Main theorem 

It was shown in Ref. [3241 ] that two permutations fi and a are combinatorially dependent 
if there is a norm preserving permutation such v that A ct (j4) = A u (A fJi (A)) for all 
operators A. One could assume that therefore Theorem 11.3,51 cannot be sharpened, 
that is, that all combinatorially independent criteria are independent. This is true when 
one considers the permutation acting on arbitrary operators. In quantum mechanics 
however, we deal with positive operators, which are Hermitian. The following theorem 
exploits this fact and counts the new upper bound to the number of independent criteria. 

Theorem 1.3.7. (i) The permutation criteria corresponding to the permutations a and 
to are dependent. 

(ii). Let Z := {e, r} be the subgroup of S2 r generated by the GQT r. Define the 
action of Z on the right cosets S2 r /T by multiplication from the left of the cosets, 
that is, e * o~T = eaT and r * o~T = tctT . The new upper bound on the number of 
independent criteria is the number of orbits under this action. It is given by [] 



2t \ IT 
r) +2r+ {r/2 ] ' <rr " (r) 



(1.15) 



where even{r) = 1 if r is even and otherwise. This number includes the trivial 
criterion given by the identity permutation. 



The integer seq uenc e generated by the number of independent criteria equals the integer sequence 
A045723 from Ref. [272 1. 
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Proof. To prove (i), let us apply the permutation ra on an arbitrary quantum state p. 
We have 



\KAp)\\ 



|A CT (A T (p))|| 
|A,(p T )|| 

IAj( 



K(p)\\ 

Mp)||. 



(1.16) 



Here we have used that p T = p because p is Hermitian and A cr (p) = A (7 ( / o) because A 
only permutes the entries of the matrix. 

Observe that multiplying a coset o~T by r from the left is the same as conjugating 
it by r because r is contained in 7", that is, we have 

tuT = raTr. 

It is readily verified that conjugation of a permutation by r corresponds to exchanging 
heads (always odd numbers) and tails (always even numbers). Therefore, the direction 
of all (true) arrows is reversed and loops and free subsystems are not affec ted. An 



32J, we call 



example is shown in Figure [T~I1 Following the four rules presented in Ref. 
this Rule 5. The new rule either glues two criteria together or does not change them. 
More precisely, the orbits under the action of Z have size 1 or 2. The orbit of size two 
is generated by the involution r: 

o~T —* tctT —* ttctT = o~T . 

(ii). There are ^( 2 r r ) combinatorially independent criteria. With the new rule, there 
are at most 

1 /2r\ 1 M /2r\ 

K 



1 2r 
2\r 



1 2r 
2\r 



independent criteria left. Here K denotes the number of criteria not affected by con- 
jugating according to the new rule. Now note that the only criteria (represented as 
disjoint arrow configurations) not affected by conjugation with r are 



1. criteria with no arrows and 
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a) 



b) 



c) 







Figure 1.4: Action of the new rule on a disjoint arrow configur ation : our new rule 



324i | which changes 



changes a) to b). Compare this with the action of rule 4 in Ref. 
a) to c). Our new rule only reverses the direction of true arrows, whereas rule 4 
additionally creates loops on free subsystems and removes loops. 



2. criteria containing only arrows and having no free subsystems. 

If r is odd, then situation (2) cannot occur. The number of these criteria is readily 
counted: 



r\ I r ' 
+ ' 1 + 2 1 + 



+ 



tr—1 



r 

l(r/2)\, 

where we have used an identity of binomial coefficients. So when the number of sub- 
systems is odd, the number of criteria becomes (including the identity) 



2r 



+ 2 r 



In the case r is even we need to take care of situation (2). Now this number equals 
picking r/2 heads from r choices, because exchanging tails does not matter. But we 
have to divide by two since exchanging all heads and tails does not matter either, so 
that the number of criteria satisfying (2) is given by 



r 
r/2 



We conclude that when r is even, the number of criteria is given by (including the 
identity) 

1 IV2r\ 

+ 2 f + 



r 
r/2 



To complete the proof, we have to show that the criteria in the orbits of size 2 are 
combinatorially independent. Let a be a permutation represented by a disjoint arrow 
configuration. Assume that there is an arrow from subsystem k to I in the disjoint 
configuration of a. Then there is an arrow from / to k in the disjoint configuration 
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describing tot. Loops and free subsystems are not affected. Using these observation 
we see that the configuration describing otot has a closed path from k to I and no 
loops. The closed path between k and I can be transforme d int o a loop on k and a loop 
on / with the help of Rule 3 (Exchanging heads) in Ref . 3241 ] . Now if we apply these 



arguments to all arrows of a we see that the permutation otcjt is not norm-preserving. 
Consequently, the permutations a and tot are combinatorially independent. This 
concludes the proof. □ 

Figure 11.41 shows an example of how the new rule glues two combinatorially inde- 
pendent criteria together. In general, we have that two permutations corresponding to 
disjoint arrow configurations are related by the new rule if and only if they have the 
same (true) arrow structure (up to exchanging heads) and a complementary loop/free 
subsystem structure. The latter means that if the first criterion has a loop on subsystem 
k then the subsystem k is free in the second and vice versa. 

1.3.4 Illustrations 

The different criteria resulting from our characterisation for two, three and four parties 
are shown graphicall y in Figure 11.51 In this figure, we have adopted the simplified 
graphical notation of 



3241 ] without arrows. Loops (QT) have been denoted by a little 
circle, (2/c, 21 — 1) permutations (R) as a solid line and {2k — 1, 21) permutations (R') as 
a dotted line. In order to construct this figure one first draws all disjoint configurations, 
and then eliminates those related by application of Rules 3, 4 and 5. 

Below we show that all listed criteria are truly independent, in the sense that no 
criterion detects strictly more states than any other criterion. A property we are 
going to use frequently is that the permutation criteria are insensitive to adding an 
uncorrelated ancilla. Let p be an r-party state and define the r + 1-party state p' = 
p<&pr+i- Suppose a € and define a' = a G S 2 r r +\-\, (the same permutation, but now 
acting on an enlarged set of 2 (r+1) elements). Then from ||A CT /(/9 / )|| = \\A a (p)(3p r+ i \\ = 
l|A,(p)|||k + i|| = \\A a (p)\\ follows that \\A a (p)\\ = HA^pOII- 
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Figure 1.5: Independent permut ation criteria for a) two, b) three, and c) four particles. 
The figure is adapted from Rcf. 324]. 
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Two parties 



For a quantum system consisting of two parties, there are only two non trivial in- 
equivalent permutation criteria: the partial transpose in one of the subsystems and 
reshuffling between the two subsystems. For the low dimensional systems 7i = C 2 ® C 2 
and TL — C 2 ® C 3 the po sitivi ty of the partial transpose is a necessary and sufficient 
con ditio n for separability 



1321 ]. while this is not the case for the realignment criterion 
2531 ] . For higher dimensional systems these criteria are truly independent, as there 
exist bound entangled states that are detected by the realignment criterion. As a cu- 
riosity, we have tested the realignment criterion on all known bound entangled states 
p £ H = C 3 (g> C 3 in the literature. The maximum value we have found for is 



7/6 and is achieved for a particular chess-board state 



3J: 
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(1.17) 



Three parties 

?or thre e parties we have proven that at most 6 criteria are independent, previous work 



96l . 13241 ] indicated that there were 9. The 6 criteria are the partial transposes (row QT) 
in the 3 subsystems and the 3 reshufflings (row R) between any of the two subsystems. 
On random states, it can easily be checked that the 6 criteria are indeed independent 
according to our definition; we now show that no criterion is strictly stronger than any 
other. 

To show that the criteria from row QT are truly independent from one another, it is 
sufficient to note that there exist tripartite NPT states which are separable with regard 
to two splits, but not to the third one. An example of these kind of states for qubits can 
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be found in Ref. [84]. Such states will obviously be detected by only one criterion from 
the row QT. To sh ow t he independence between the criteria in the row R, we use the 



property (see also [2531 ]) that the trace norm of the realigned density matrix remains 
invariant when an uncorrelated ancilla is added. Now taking a bipartite entangled state 
which violates the realignment criterion, and adding an uncorrelated ancilla, we find 
that it is only detected by only one criterion from row R (since it is separable in the 
other splits). 

All what is left is to show that there are no dependences between the rows R and 
QT. So let us take two arbitrary criteria, one from row R and one from row QT. Here 
we can distinguish two cases, the first is where the realignment works on different 
subsystems than the partial transposition. As an example let us take realignment R12 
between subsystem 1 and 2, and partial transposition QT3 of subsystem 3. Let p be 
a bipartite entangled state detected by both the realignment and the partial transpose 
criterion. Then the state p' = p P3, with p^ an arbitrary state, will be detected by 
R12, but not by QT3. Conversely the state p" = p\® p will be detected by QT3 but not 
by Ri2- The second case, is where the realignment works on a same subsystem as the 
partial transposition, such as R12 and QTi. To show their independence it is sufficient 
to consider bipartite states that are detected by either R or QT (but not both), and 
adding an uncorrelated ancilla. 

Note that the realignment criterion, in contrast to the partial transpose criterion 
can detect tripartite entangled states, that are sep arable under any bipartite cut. This 



144| with the tripartite bound entangled 



has been demonstrated experimentally in Ref. 
states from Ref. [23]. 

Four parties 

For four parties, there are at 22 non trivial independent permutation criteria. As for 
three parties, it easy to show the independence of the criteria within a fixed row, by 
starting from an entangled two or three party state and adding uncorrelated ancillas. 
This is probably hardest to see for the row R+QT, so we illustrate this with an example. 
Let us take the two criteria R12QT3 and QT1R23 and let /?23 be a bipartite entangled 
state which is detected by the realignment criterion. The state p\ (g> P23 <8> Pa (with 
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Pi,P4 arbitrary states) is detected by QT1R23, but not by R12QT3, as this criterion is 
equivalent to the criterion R'i2QT4. 

Let us now focus on the criteria which contain only partial transpositions. The rows 
QT and 2QT are easily shown to be independent by adding uncorrelated ancillas. As 
an example consider the criteria QTi and QT12. Let pu be an NPT state, and add an 
uncorrelated ancilla ^34; clearly this state will be detected by QTi only. Conversely, 
constructing a state as the tensor product of an NPT entangled state P24 and a separable 
state p\3 will give a state only detected by QT12. The independence between the criteria 
with only loops and the criteria with at least one realignment follows easily from their 
independence in the bipartite case. As before this is done simply by adding suitable 
uncorrelated ancillas. 

To show true independence within the set of realignment criteria (rows R, R+QT, 
2R and R+R'), let us first consider the rows R and R+QT. Again, we consider an 
example, whereas the other cases can be shown similarly. Consider the criteria R12 
and QT4R23. Let p be a bipartite entangled state which is detected by the realignment 
criterion and a a separable state. Then the state pu <8> 034 will be detected only by 
R12, while 0"i4 (8> P23 only by QT4R23. As for the criteria 2R and R+R', consider the 
non-trivial case R12R34 and Ri2R'34- Their independence can be shown by applying the 
criteria to the 4-qubit states V = (|1000) + |1101))/\/2 and V' = (|1000) + |1110))/\/2 
(the first is detected by R12R34 only, while the second by Ri2R'34 only). All that is left, 
is to show that the criteria from row R and R+QT are independent of the criteria from 
rows 2R and R+R'. Without loss of generality we can take one representative from 
the rows 2R and R+R', say R12R34 (the other criteria correspond to appropriately 
swapping the subsystems). It is easy to find states that are detected by the criteria in 
rows R or R+QT but not by Ri2R34- Consider for example the criterion R12QT3. One 
can verify that it will detect the entangled state p = p c 1/9; while this state is not 
detected by criterion Ri2R34- Finally an example to show that the criterion R12R34 
can detect states not detected by any criterion in the rows R or R+QT is given by 
p= (3/4)p c ®p c + (1/4)1/81. 
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Final remarks 

In Ref . 



144i | a powerful class of separability criteria was devised based on permutations. 
The class however contained many redundancies, and to give a com plete characteriza- 



324| | a graphical repre- 



tion of the independent criteria is an open problem. In Ref. 
sentation for permutations together with rules for simplifying them were introduced 
based on a sufficient condition for two permutations to yield dependent criteria. This 
equivalence meant that two combinatorially dependent criteria yield the same value of 
the trace norm on all operators. Combinatorially, density operators have a prominent 
Hermitian symmetry, that is a global quantum transposition together with a complex 
conjugate. We have exploited this symmetry and we have shown how this led to the 
dependence of particular combinatorially independent criteria. It is unlikely that there 
are more dependences in the criteria from Theorem 11.3.71 

Density operators differ from arbitrary operators also because they have positive 
eigenvalues. But since permutations only reorder matrix entries, it is unlikely that this 
positiveness would lead to more criteria to be dependent. We have verified the inde- 
penden ce of the criteria numerically on a random state (using the algorithm outlined 



in Ref. 



342| ) for 2 up to 8 parties. In Appendix lA.il we have reproduced the matlab 



program which implements a random permutation. In the same way as we illustrated 
for three and four parties, it is easy to see that the criteria with only loops (only partial 
transpositions) are independent. These criteria are independent from the ones having 
at least one realignment since those can detect bound entangled states. To prove the 
independence within the class of criteria with at least one realignment one could try to 
generalize the arguments from the previous section. 



CHAPTER 2 



Schmidt number for mixed states 



286], generali- 



In this chapter we introduce the Schmidt number for mixed states 
sation of the Schmidt rank for pure states. It can be thought of as a refinement of the 
entangled versus separable division, where the set of mixed states is decomposed into 
a whole series of nested convex sets. In th e firs t section, we study briefly the use of 
Schmidt witnesses to characterise these sets 2571 ] . In the succeedin g sec tions we derive 



95 



2861 ] and discuss a 



the corresponding characterisation in terms of positive maps 
conjecture on the decomposability of these maps in 3 ® 3. 

Robustness measures quantify the extent to which entangled states remain entan- 
gled under mixing. In the final part of this chapter, we generalise the notion of robust- 
ness of entanglement 



robustness 



304| and introduce the Schmidt robustness and random Schmidt 



611 ]. The latter notion is closely related to the construction of Schmidt balls 



around the identity. In particular we derive upper and lower bounds to the Schmidt 
robustness for pure states. We present two conjectures, the first one is related to the 
radius of inner balls around the identity in the convex set of Schmidt number n-states. 
We also conjecture a class of optimal Schmidt witnesses for pure states. 
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2.1 Introduction and Schmidt witnesses 
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2.1 Introduction and Schmidt witnesses 



In Ref. [2861 ] the notion of Schmidt rank S(ip) of pure bipartite states was extended 



to the Schmidt number S(p) of mixed states p as 

S{p)= min maxStyi). 

p=T,iPim)(m 1 

Thus p acting on Tt = TCa <8> 7~Cb is said to have Schmidt number n if there exists a 
decomposition of p = Y^i Pil^Pi) i^il with all vectors {IV'i}} having Schmidt rank at most 
n, and there exists no such decomposition with all vectors having a Schmidt rank n — 1 
or lower. This definition coincides with the Schmidt rank when p is a pure state. It 
is clear that separable states have Schmidt number one, whereas entangled states have 
a Schmidt number larger than one. Recall that we denoted the set of separable states 
as 5, and the complete set of states as T>. We will denote the set of density operators 
that have Schmidt number n or less by S n . Sets of increasing Schmidt number are 
embedded into each other as 

S = Si C S 2 C . . . C S d = V. 
The subsets Si are all closed and convex by construction. The Schmidt ran k of a pure 



1901], and 



state cannot increase under local operations and classical communication 
this holds true for the Schmidt number as well. Furthermore, the Schmidt number of 
a state is trivially invariant under local unitaries. This makes the Schmidt number k 
of a state p, or equivalently the logarithmic Schmidt number M(p) = log k a legitimate 
measure of entanglement. Most surprisingly, Terhal and Horodecki [286] showed that 
N{p) is not an additive quantity, that is there exists states p with Schmidt number k, 
such that 

N{p® n ) ^nlogk, 
so that in general all we can say is that 

M{p) < N{p® n ) < nM{p). 

For multipartite systems, there is no straightforward generalisation of the Schmidt 
decomposition. However there have been a number of proposals, both for pure and 



mixed states (see e.g. [40j, |90|): 
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In analogy to entanglement witnesses the s et of sta tes in can be characterised 
with the aid of the so-called Schmidt witnesses 



39 



2571 ]. We call a Hermitian operator 



W a Schmidt witness of class k (for short fc-SW) if and only if 

1. Tr(Wa) > for all a G S k _ u 

2. There exists a p £ Sk such that Tr(Wp) < 0, 

The existence of a fc-SW for a Schmidt number fc state follows directly from Theo- 
rem [TX2 

Corollary 2.1.1. A state p £ T> has Schmidt number k if and only if Tr(pWk+i) > 
for all Schmidt witnesses W^+i of class k + 1. 

In a completely similar way as the discussion in Section 11.2.31 on PPT entangled 
states and their witnesses, every Schmidt number k state can be decomposed as the 
convex sum of a Schmidt number k — 1 state and a Schmidt number k edge state 
(see Theorem I1.2.7P . Furthermore for every Schmidt number k edge state there is a 
canonical witness that detects the state and conversely, every fc-SW can be written in 
a canonical form with the aid of an Schmidt number k edge state (see Theorem ll,2,8p . 
Since the theorems and proofs are almo st identical we will not reproduce them here. 

Consider the following example 286| of a one-parameter set of states exhibiting all 



possible Schmidt numbers and a corresponding set of fc-SWs. 

Example 2.1.1. Let P + be a maximally entangle d state a cting on a Hilbert space 



7i = C d ®C d and the unnormalized isotropic states 



Mi, 



24 A] defined by 



pp = 1 + pp +i with - 1 < < 00. (2.1) 



The isotropic state pp has Schmidt number n if and only if 

d{{n - l)d - 1) < < d(nd - 1) 



d — (n — 1) d — n 

The operators 

W n = l-d/(n-l)P+, (2.2) 
are n-SW and detect in an optimal way the Schmidt number of the states pp. 
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Let us for a moment assume that W n is indeed a n-SW . Then it is not so hard to 
verify that when (3 > d(nd—l)/(d—n), pp has Schmidt number n+1. The example then 
follows from a continuity argument, if we can show that for j3 = d{nd — l)/(d — n) the 
state has Schmidt num ber n . This can be done by virtue of an explicit decomposition 

ts in 



in terms of pure states 



2861 ] . Such a decomposition can be deduced from our resu 



Section EH The fact that W n is a n-SW follows readily from the next lemma 

Lemma 2.1.1. Define the fully entangled fraction of a state p as 

F{p) =max(#|p|#), 

where the maximum is taken over all maximally entangled states, 
(a) For a pure state \ip) the fully entangled fraction is given by 

d 



131 



2m- 



|2 

"i=l 

with \{ the Schmidt coefficients of\ip). 

(b) For an arbitrary state with Schmidt number k or less we have 



Proof, (a) See Ref. [1311 ]. 

(b) Suppose we have a pure state \ip) with Schmidt rank k, so the problem is to 
maximise 

h n 

i=i 

subject to 

k 
i=l 

We solve this using the method of Lagrange multipliers and define 

k ,2 



i=i i=i 

An extremum can be found by solving 
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The first relation leads to \i = —k, and thus Y2i = kXj, squaring this equation and 
carrying out summation over j yields 

k 



i=l 



If the state is not pure, this inequality remains true, since a mixed state of Schmidt 
number k can be written as a convex combination of pure states of Schmidt rank k. □ 



In Ref. 1601 ] a connection was made between Schmidt witnesses and entanglement 
witnesses. The central idea in their construction can be deduced from the following 
lemma 



more genera, 
the years 83, 



123, 



60l | . of which special cases have appeared in the literature over 



3101 ] (see also the recent [227j |) 



Lemma 2.1.2. Let a be a positive operator with Schmidt number N > 1 acting on 
7~Ci = T~tAi G^Hbi and let rj be an operator acting on Ti. 2 = ~Ha 2 ® Wb 2 positive on states 
with Schmidt number KN. Then the operator a (g) 77 acting on Hi <8> H 2 — Ha <8> H.B is 
positive on states with Schmidt number K . 



Proof. It is clear that it is sufficient to prove the lemma for pure states a 
let 

JV 

\4>) = S ^(t>i\ai)A 1 \hi) Bl , 



. So 



and take an arbitrary Schmidt rank K state 

K 

|^> = ^2^j\ e j)A\fj)B- 
3=1 

Then we need to prove that Tr (\ip){ip\o~®r}) > 0. This trace operation can be composed 
of tracing out the first pair, and then the second, as Tr(-) = Tr2(Tri(-)). Then making 
use of the identity Tv 1 (C(A 1 B 2 )) = Tri(C(^i 1 2 ))B 2 we have 

Tr(|^)(^|(T ®rj) = Tr 2 (Tn(|V;>(V|(<x ® l 2 ))v) 
Now Tri(|^)(^|(<r (8) I2)) is the projector onto the pure Schmidt rank KN state 

|7)2 = (M) =Y,<lHM a i\ e 3)Wj) e ^2 
ij 

from which the lemma follows. □ 
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Let S be a given (k + 1)-SW and P\_ a maximally fc-level entangled state, then 
from the above lemma follows that the operator S' defined by 

S' = pf } ® 5, 

is an entanglement witness. Thus by enlarging the Hilbert space, any SW can be 
mapped onto an entanglement witness; a trick which will come back in Chapter [H 



2.2 Main theorem 

In this section we prove the main result on fc-positive maps: the operators corresponding 
to fc-positive maps via the Jamiolkowski isomorphism (see Section Tl. 2. 2 p are positive on 



Schmidt number k states and vice verse. This is a generalisation of Theorem ll.2.41 Then 
we derive a characterisation of in terms of positive maps, similar to Theorem ll.2.51 
We start this section by giving some alternative characterisations of positive maps. 

Recall that a map A : C-Ha ~~ *■ ^Hb ^ s a ^-positive map if and only if it the induced 
map 

Afc = Ifc <8> A : C C k <g> C Ha -> £ C fc ® £ Hb 

is positive for k > 1. The next lemma characterises ^-positive maps purely on their 
action on pure states with Schmidt rank at most k (see for instance Ref. 95]). 



Lemma 2.2.1. A Hermitian map A : C-Ha ~* ^H B * s k-positive if and only if 

(l fc ®A)(|V)M) >0, 
for all pure states \tp) 6C fc ® TLa with Schmidt rank at most k. 

Proof. Trivial. □ 
An improvement is made by the following lemma, which states that it is sufficient 



to test maximally entangled vectors of rank k. Ref. 



72 ] contains a proof for k = 2. We 



present here an alternative, much shorter proof for general k. 

Lemma 2.2.2. A Hermitian map A : £-^ A — > C-h b * s k-positive if and only if 

(i fc ®A)(|v}(VI) >o 

for all maximally entangled states G C k (g> TLa- 
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Proof. Let 6^8 Ha be an arbitrary state in Schmidt decomposition: 

fc-l 



A ;l e *)l/i' 



i=0 

In view of the previous lemma we need to show that 

D = (l k ®A)(\ip)(ifj\)>0 if and only if D = (l k ® A)(|^)(^|) > 0, 

with 

fc-l 

i$ = £i*>i/*>- 

But this follows immediately from 

L> = {L®t)b{L ] 

with L the diagonal matrix such that La = \. □ 
Now we are in the position to prove the following theorem by Uhlmann. The proof 



is from Ref. 



2861 ] 



Theorem 2.2.1. A Hermitian map A : Cn A ~* ^W B ^ s k-positive if and only if 

k 



y~] ^n^m{bn\A{\a n ){a m \)\b m ) > 



n,m=l 

for all orthogonal sets of vectors {|a n )}^ =1 and {\b n )}!^ =1 with X) re =lA l n = 1 ■ 

Proof. From the preceding lemma we have that a Hermitian map is fc-positive if and 
only if 

M(l fc ®A)(WMM >0, 

for all cj) and ip with Schmidt rank k and t/j maximally entangled. Using the Schmidt 
decomposition we set <fi = X^jWl *^) an d V* = ^2i\°i a i) (note that this is always 
possible, since ip is maximally entangled) and obtain 

k 

/W«m(crAi|lfc ® A(|c i Oj)(c : ,-a : ,-|)|c m 5 m ) > 

i ,j,n,m=l 

from which the theorem follows. □ 



With this we can prove the main theorem, first published in Ref. 



286] 
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Theorem 2.2.2. A linear map A\y : £t-Ca ~~ * ^-'Hb * s k-positive if and only if the 
corresponding operator W acting on Ha'&'Hb is positive on states p acting on 'Ha^'Hb 
with Schmidt number k or less (belonging to Sk)- 

Proof. Let 

W = {t®k w ){P + ), 
we have to show that fc-positivity is equivalent to 

TrW>>0, (2.3) 

for all p £ Sk- It is sufficient to take p as a pure state of Schmidt rank k, that is, 

k 

P= P-mPn\a m b m ){a n b n \. 

n,m=l 

With this (12, 3p can be rewritten as 

k 



^mPn{a n bn\W\a m b m ) > 0, 



n,m=l 

or 

PmPn{a n \i){a m \j){b n \K{\i){j\)\b m ) > 0. (2.4) 

i ,j,n,m=l 

Due to the linearity of A we have that A(\a n )(a m \) = YliijiAo-n) ( a m\j) A{\i) , or 

A'(\a n )(a m \) = J2(a n \i)(a m \j)A(\i)(j\), 

with A' = A o R ? where R maps the vectors \a n ) to |a*) where complex conjugation is 
performed with respect to the {\i)} basis. Therefore (|2.4p is equivalent to 

k 



^ /%Mm(&n|A'(|a n )(a m |)|6 m ) > 0, 



n,m=l 



and it follows that A' is /c-positive if and only if Tr Wp > for all p £ Sk- Now the map 
R corresponds to a unitary rotation which maps the basis {|a n )} to the basis {|a^)}. 
Therefore A is fc-positive if and only if A' is. □ 
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Thus in the language of Schmidt witnesses: A map is fc-positive if and only the 
corresponding operator is a k + 1 Schmidt witness. As an example consider the k + 1- 

SW (1221) 

W k+1 = t-±P+, 

the corresponding map is given by Af t (X) = Tr(A)l — kX. The positivity pro perti es 



of this map have been derived independently by Tomiyama [283] (see also Ref. [28 j) . 

but from our correspondence they follow directly. Note that for k = 1 we obtain just 
the reduction criterion. 

Theorem 2.2.3. Let A p be a linear map defined by 

A p (X)=Tt(X)1- P X, (2.5) 
be a linear map. Then A p is k-positive for 

11 , X 

<P<7> ( 2 - 6 ) 



k + 1 ~ k 
but not k + 1-positive. 

Next we give a characterisation of Sk in terms of positive maps, similar to Theo- 
rem [TT2T5]. For the proof one needs the following lemma. 

Lemma 2.2.3. A linear map A : C-n A — * C-^ B is k-positive if and only if A* : Cn B — > 
Cn A is k-positive. 

Proof. Suppose A is A>positive, then we have 

M(l fc ®A)(|$MM >0, 

where \4>) G Wfc ® "Ha and £ 7^ ® lis are arbitrary and both have Schmidt rank 
k. Applying the definition of an adjoint map yields 

M(l*®At)(MM)|$ >0 

for arbitrary |0) and of Schmidt rank k. □ 

Corollary 2.2.1. A density operator p acting on H.a'&'Hb belongs to Sk if and only if 

(1® A)(p) > 0, 

for all k-positive maps A : ~~ * ^Ha- 
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Proof. This follows from Theorem 12.2.21 and Corollary 12.1.11 in a completely analogous 



way as we derived Theorem 11.2.51 □ 

2.3 A conjecture 

The detection of PPT entangled states is a difficult problem and this suggests that 
PPT entangled states have generically a low degree of entanglement. In particul ar we 
expect PPT entangled states to have a relative low Schmidt number. In Ref. 8^, 257 ] 
evidence was given for the following conjecture. 

Conjecture 2.3.1. All entangled PPT states on TL = C 3 <S> C 3 have Schmidt number 
two. 

It is easy to see that this conjecture is equivalent to the statement that all 2-positive 
maps on C^3 are decomposable. A first piece of evidence is that the conjecture holds 



for PPT entangled states of rank 4 



257] 



Lemma 2.3.1. All entangled PPT states on H = C 3 (g> C 3 of rank four have Schmidt 
number two. 

Proof. Obviously it is sufficient to show this for edge states. So let p be a rank four 
edge state. We will show explicitly how to write p as a convex sum of Schmidt r ank two 
vectors, by iterative subtracting of such vectors (in this context, see also Ref. 182]). 



Since in C 3 (g) C 3 there exists a UPB of 5 product states there must exist a product 
vector |e/i) in the kernel of p. Since Tr(/9|e/i)(e/i|) = Tr(p Ts |e/j f )(e/*|) = and 
p TB > we have that \efl) belongs to the kernel of p Te . 

Constructing an orthogonal basis {]/»)} we have that {f*\p TB \efi) = = (fi\p\efi), 
so that for instance for i = 2, 

|^i>=p|e/2> = |c/2> + |4T3>. 

Since this vector has at most Schmidt rank two and belongs to the range of p, we can 
subtract it and obtain a new state 

P = P~ AlV'iXV'il, 
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with 



- = {ij) X \p Vi) = <"0l|e/2>. 



The state p has rank three and it is easy to verify that p\efi) = ple/2) = but 

|^2>=p|e/ 3 )Hc/2> + K/3>. 

This vector has again at most a Schmidt rank two vector, and we can again subtract 
it to obtain 

p = p- A|V>2> <-02 1 , 

with 

Y = (V^lp -1 ]^} = (^2|e/ 3 ). 
A 

The state p has rank two and p\efi) = ple/2) = p\efz) = 0. From this follows that we 
can write p as 

P = P + k\i) 2 ){ii)2\ +A|^i)(^i|. 

Since p effectively acts on C 2 ® C 3 it has Schmidt number two and therefore p has also 
Schmidt number two. □ 

In order to check the conjecture for arbitrary rank, first note that again it is sufficient 
to check it for edge states. Let W denote the set of all entanglement witnesses detecting 
these edge states. Then the conjecture is equivalent to the statement that for all W G W 
on C 3 <8> C 3 there exists a Schmidt rank two vector \ip) such that ( - 0|W / |V>) < 0. Now 



inC 3 ®C 3 the sum of the ranks of p and p Te are not larger t 



151 



ran 13 for edge states 



2571 ] gave an arg ume nt 



1521). 



1761 ]. For states with such a relative low rank Sanpera et al. 
that such a \ifj) always exists (for a generalisation to higher dimensions, see Ref. 
Their argument is based on counting the degrees of freedom versus the number of 
constraints in constructing such a \ip). Their argument is especially weak for ranks 
(6, 6) and (8, 5) since here only a finite number of solutions can be found, for which it 
is by no means guaranteed they will satisfy the constraints. When the conjecture was 
first published only PPT edge states of rank (4, 4) and (6, 7) were known. In view of 
the conjecture, there was special interest in constructing PPT edge states with ranks 



(6,6) and (8,5). In a recent paper Ha and Kye 



1131 ] managed to find edge PPTES 
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for all ranks except (5, 5) and (6, 6). It is for this reason we have specially constructed 
such states in Section 11.2.41 

For these states, we can construct witnesses in a number of different ways. First 
note that both states violate the realignment criterion, and thus using Theorem 11.3.41 
we obtain W\. Secondly, there is the canonical construction for edge states of Theo- 
rem [L2]8j which gives us Wi- Thirdly, one can use the algorithm [3] from Section fl.2.51 
yielding W3 . Numerically we calculated all three witnesses W\ , W2 and W3 for p^^) 
and />(6,6)- After normalisation we found that all witnesses are different and yield a 
different value on the states. In total this gives us six different non-decomposable wit- 
nesses. We have checked that these witnesses are negative on some Schmidt rank two 
state. This provides additional evidence in favour of the conjecture. In fact, we also 
verified this for the witnesses 

W t = Wi - (Tr(W iP ) + e)l, 

for arbitrarily small e > 0. These witnesses Wi barely detect our states and are therefore 
the best candidates to check the conjecture. Note, however, that this is only evidence 
that the states have Schmidt number two. A full proof would consist of an explicit 
decomposition in Schmidt rank two states. However, finding such a decomposition is 
in general a very hard problem and we did not attempt this. 



2.4 Schmidt robustness 



In this section we generalise the notion of robustness of entanglement (see also Sec- 
tion I3.3[) to the robustness of Schmidt number, which we anal yse f or pure s tates. We 



start by recalling the definition of robustness of entanglement 



117 



274 



304|] 



The fC-robustness of a state p, Rk(p) > is the minimal value of R such that 



1 



-{p + Rpk) 



1 + R 

is separable, for some state p^ G /C. Thus robustness measures measure how much 
mixing is required before a state becomes separable. With this basic definition we have 



The robustness R s (p) Here tC = S, the set of separable states. 
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The random robustness R r (p) Here K, = {1}, the maximally mixed state 

The generalised robustness R g (p) Here K = T>, the set of all states. 

From the definitions it is clear that R g < R s < R r . Both the robustness and the 
generalised robustness are entanglement monotones (see Section I3.3f) . It is easy to 
see from the definition that the random robustness is proportional to the so-called 
witnessed entanglement J^, defined as maxyKeA/[0, — Tr(pVF)] where M is the set of 
all normalised entanglement witnesses. The following theorem gives exact values for 
the robustness and the random robustness for pure states. 



Theorem 2.4.1 ( [1171 . 1274 304]). Let \ip) = Yli a i\^) be a pure bipartite state with 



ordered Schmidt coefficients <Zj > <2j+i. The robustness of ip is given by 
Rs{tp) = y^«jQj = C^ai) 2 - 1 = R g (ip) 

i^j i 

and equals its generalised robustness. 

A separable state that washes out the entanglement in if) most quickly is given by 

p s = p g = -jr^aiaj\ij)(ij\. 

The random robustness of tp is given by Rr{ip) = a±a2- 

An alternative proof of this theorem can be obtained as a corollary of our results 
(see Section [2.4.1|) . In what follows, we assume that all states act on a Hilbert space 
H = C d <g> C d . Throughout we also use local operations and classical communication 
(LOCC, see Chapter [3}, as they do not increase the Schmidt number. 

2.4.1 Generalised Schmidt robustness 

In this section we will extend the notion of generalised robustness to generalised Schmidt 
robustness. Analogously to the generalised robustness, the generalised Schmidt-n ro- 
bustness of a state p, R gn (p) is defined as the minimal value of R such that 

1 



1 + R 



(p + Rpg r , 



For simplicity we have chosen not to n ormalise 11 in the definition, in contrast to the original 
definition of random robustness in Ref. 304 1. 
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has Schmidt number smaller or equal than re, for some p gn . 

Let us now analyse the generalised Schmidt robustness for pure states. In Theo- 
rem E3?T] we have seen that for a pure state = a>i\ii) the state 

Pg = -p-^2 t Ha j \ij)(ij\, (2.7) 



Rg 

erases most quickly the entanglement present in It is plausible that this same state 
will also erase Schmidt number in an optimal way. We were able to prove this for the 
maximally entangled state. For the maximally entangled state = 7/jX)i K*)> Pg 
given by 12.71 reduces to p g = , (t — Z), with Z = ^ 

Theorem 2.4.2 (Generalised Schmidt robustness of the maximally entangled state). 
The state defined by 

M-Zffi whm Ps = 10 d (2.8) 



has Schmidt number n for 



d — n „ d — n + 1 

< < 



n n — 1 

The generalised Schmidt-n robustness of the maximally entangled state P+ is given by 

Rgn(P+) = ~ ■ 

Proof. Let S(j3) be the Schmidt number of p(/3) and let j3 n = We first show that 
S(f3 n ) < n. 

(i) We will give an explicit decomposition of the state p{(5 n ) in terms of Schmidt 
rank n states. Equivalently, we show how one can construct this state locally with the 
aid of Schmidt rank n states. In what follows we will often omit normalisation. Let us 
take a maximally entangled n-level state 



1 n 



n 



where S = . . . , i n } is a subset of {1, ... , d}. We can construct such a state in (^) 
possible ways, and clearly all these states have Schmidt number re. Now let us mix with 
equal weight the corresponding states \i^s){' t Ps\- Then for every i and j we will have 
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( _l) terms of the form and ( n _|) terms of the form i ^ j. Thus the 

resulting state will be proportional to 

(d-l)Z+(n-l)J2\n)(jj\- 

Therefore we have proven that the state (see Ref. {59] for n=2) 

A . = z + ^-i) 

a — n 

has Schmidt number at most n. It turns out that we can transform this state in the 
form (|2.8p by applying a certain partial twirl operation. Consider the following twirl 



operation 



131 



241 



309[] 



which maps any state p into one of the form 1 + aP + (an isotropic state). Here dll 
is the uniform probability distribution on the unitary group U(d). Note that the twirl 
can be implemented locally, both parties need to implement only one random unitary. 



Remarkably, it has been shown 



83 □ that the integral can be written as a finite sum 

n 



(for qubits this was first shown in Ref. [2l|). This will allow us to perform a partial 
twirl, just by considering a part of this finite sum. 

In the first step we apply the unitary transformation (T (g) T*)K(T (g> T*y , with 

d-l 



r a— i 

-^2^ ed b)< fc l. 



j,k=0 

which i s .jus t the quantum Fourier transform. Since T is unitary, it acts as the identity 
on P 4 



1311 ] . while it acts on Z as 

(T (8) T*)Z{T T*) f = Yl Yl 



i2ir(j — t-{-k — s)a 

e a \j)(t\®\s)(k\. 



a j,k,s,t 

The terms of the form \jj)(kk\ for s = j and t = k will give a contribution P + , while 
the for i ^ j (for j = t and k = s) will give a contribution of (1 — Z)/d, so that 



K' = d(T ®T*)K(T ®T 



1-Z+«$=&P ++ L, 
d — n 



2 In Ref. 



the finite decomposition of the U ® U twirl was given, but a similar decomposition can 



be deduced from it for the U ® U* twirl. 
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where L are terms not of the form \ii){jj\ or \ij){ij\. Now these contributions can be 
easily removed by repeated application of the following mixing procedure 

K" = -{U® U*)K(U ® U*) ] + -K'. 

First U is chosen to act as U\k) = e l7:Skl \k) for every I = 0, ...,d— 1. This defines 
d mixing procedures. Next U is taken to act as U\k) = e lnSkl ^ 2 ^ (another d mixing 



procedures). One can readily check 83] that these operations do not affect terms \ii){jj\ 
or \ij){ij\ but cancel out L completely. 

Thus S((3 n ) < n. Now for (3 n < (3 < the state p(f3) is a convex combination 

of p((3 n ) and and therefore S((3) < n. 



113]. For any 



(ii) For the second part, we generalise the trick introduced in Ref. 
state a, suppose t is a positive number such that P + + to has Schmidt number n. The 
operators W n = 1 — d/nP + witness Schmidt number n + 1 (see Example 12. l.ip . so that 
we have 

< Tr [(1 - d/nP + ){P + + ta)\ 

= 1 + t - d/n Tr [P + ] - d/n Tr [P+o] 

d — n 
< +t. 

n 

since TrP + cr > 0. Thus t > (3 n , and for a = p g it follows that S(f3) > n for (3 n < (3 < 

For general a, it follows that R gn > (3 n , but p((3 n ) = P + + j3 n p g has Schmidt number 
n, so that R gn {P + ) = (3 n . □ 

Note that the states (I2.8P constitute one of the very few examples of non-trivial one 
parameter states for which the Schmidt number is known. To our knowledge, the only 
other example is that of the isotropic states, Example 12.1.11 This theorem allows us 
to present non-trivial bounds to the generalised Schmidt robustness of arbitrary pure 
states. 

Corollary 2.4.1 (Lower and upper bounds for the generalised Schmidt robustness). 
The generalised Schmidt-n robustness R gn of a pure state \ip) = X^ a iN) satisfies 



1 < Rnn < Ra 



d — n 
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with R g = ^2i^j(iiaj = (X^ a «) ~~ 1 ^ e generalised robustness ofip. 

Proof. In Theorem 12.4.21 we have seen that p = ^-jrzdit - Z) + P + has Schmidt 
number n. Performing the filtering operation (A <8> A)p(A (g) Ay we hence obtain a 
state with Schmidt number at mos t n (this is because local filtering cannot increase 



the Schmidt number of a state 



28a]). If we take A = y/ak\k){k\ we obtain that 



has Schmidt number n. Here p g = ^- Yli^j a i a j\^j)(^j\ as before. This gives the upper 
bound. The lower bound can be readily proven using exactly the same trick as in part 
(ii) of Theorem EMI □ 

The lower and upper bound only coincide when ip is the maximally entangled state 
or when n = 1. Note that the lower bound can be negative. The upper bound depends 
on the dimension of the Hilbert space in which the state is embedded, and hence will 
in general not match the value of the generalised Schmidt robustness. 

2.4.2 Random Schmidt robustness 

We define the random Schmidt-n robustness of a state p, R rn (p) as the minimum value 
of R such that 

' (p + Rt) 



l + R 

has Schmidt number n. 
Upper bounds 

As we have seen from Example 12.1.11 for p = P + we have R rn = (d — n)/[d{nd — 1)]. 
For general pure states a (weak) upper bound to the random Schmidt robustness can 
be obtained as follows. We know that 

T n = (d- n)l + (nd - l)dP+ 

has Schmidt number n. Local filtering (A® A)T n (A® Ay cannot increase the Schmidt 
number. So that with A = Y2k V^I^K^I we obtain 

T' n = (d- n)p A ® pa + (nd - l)\ip)(ip\, 
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with pa the reduced density operator of \ip)(tfj\, and because 1 — pA ® PA is a separable 
state we get as an upper bound for the random Schmidt robustness R rn (ijj) < (d — 
n)/(nd — 1). The following theorem presents a non-trivial upper bound. 

Theorem 2.4.3 (Upper bound to the random Schmidt robustness). The random 

Schmidt-n robustness R rn of a pure state \tp) = 'Ylnd^ii) satisfies 

Rr{d-n) 
an — 1 

with R r = a\a2 the random robustness of ip. 

Proof. In this proof, we work again with unnormalised states. Note that, if we add 
two unnormalised Schmidt number n states together we end up with a state of at most 
Schmidt number n. Equivalently, this can be shown by mixing the normalised states 
with different weights. 

Analogously to the construction in Theorem l2.4,2l let us take a maximally entangled 
n-level state 

n 



\C)S = ~ V". 

n 

i&S 

where S = {it, . . . , i n } is a subset of {1, ... , d}. Again, we can construct such a state 
in („) possible ways, and mixing the corresponding (unnormalised) states \ips){^s\ 
together, we end up with a state proportional to 

Y^a*\ii)(ii\ + ^^mM, 

i 

which has at most Schmidt number n. Now, in Corollary 12.4.11 we have seen that the 
state 

d — n 

has Schmidt number at most n. Adding these two states together we find that 

dn 



~W){^\ + ^aiaj\ij){ij\ 



d — n 

has Schmidt number no more than n. Mixing this state with the separable state 
^2ij( a i a 2 — a j a j)\'ij){'ij\, we obtain finally 

dn — 1 
(a — n)a\a2 

□ 
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Conjecture 2.4.1 (Schmidt balls around the identity). Consider the unnormalized 
mixture 

pp = i + Pp, with - 1 < < oo, (2.9) 

with p an arbitrary normalized state acting on C d ® C d . The state pp has Schmidt 
number at most k for all d and p when 

P<2(2k 2 -1). (2.10) 

Note that it is a generalisation of the fact that 1 + 2p is separable for all normalised 



states p 304]. The conjectured value (|2.10p can be obtained as follows. Starting from 
the identity, it is natural to assume that we can go fastest to a higher Schmidt number by 
mixing with some maximally entangled state. Taking p = Pf (the maximally entangled 
cf -level state) we have that p < d'(kd! — l)/(d' — k), which reaches its maximum for 



d' = Vk 2 — 1 + k or for integer d' = 2k. Substituting this expression in the expression 
for p gives the upper bound from the conjecture. Thus it looks like, starting from the 
identity, we can get most quickly to a Schmidt number k + 1 state by mixing with the 
maximally entangled state in 2k (g> 2k. Upper bounds on the random Schmidt number 
two robustness are particularly useful as a distillability criterion can be deduced from 
it, which we will discuss in Section I4T21 

Lower bounds 

Lower bounds on the random Schmidt robustness of a pure state ip can be obtained from 
any Schmidt witness that detects ip. Indeed, suppose Wn+i is a normalised Schmidt 
number n + 1 witness such that Tr(|^)(^|W n _|_i) = —a, with a > 0. Then it is easy to 
see that R rn (ip) > -p. It follows that 

d 2 R rn (^) = - min TrGV^lWn+i). 
Tr W n +i=l 

The question is now to determine this class of optimal Schmidt witnesses. We have 
conducted a large numerical test and have found a class which we conjecture to be 
optimal. The construction is the following: 

Consider the Schmidt witnesses W n +i = 1 — d/nP+ from Example l2.1.1l Performing 
the filtering operation (^4 (%> B)W n +\(A By we obtain again a Schmidt number n + 1 



2.4 Schmidt robustness 



66 



witness. Let us consider the particular case where A = B = Ylk V^I^K^I diagonal 
and such that ^ a? = 1; we obtain 

T?/ "EyOiajlu'Xy'l -Ey a i a iR0'j'l , 011 v 
= n(E^-l • ( 2 - U ) 

For a given pure state = Yli bi\ii) we therefore conjecture: 

Conjecture 2.4.2. The random Schmidt robustness R rn of a pure state \ip) = bi\ii) 
is given by 

d 2 R rn (iP)=- min TV(|^)^|W n+1 ) = - nun ^ .* » ^ * t} , (2.12) 



280[]. The 



One can solve the minimisation problem with Lagrange multipliers 
solution is given by 

min Tr(|V>)(V'|W / n+i) = -A min , 

TrW„+i=l 

where A m i n is the minimum eigenvalue of 

»-;£t)<" d - b >' 

with J the all-ones matrix, D the diagonal matrix with Da = bf and I? the matrix 
with coefficients Bij = bibj. 

A first step in proving this conjecture would be to show that the class of witnesses 
(A® B)W n +i(A l g) By with A,B arbitrary matrices, is no more powerful than the class 
of witnesses (A® A)W n +i(A® A)^ , with A diagonal. We have numerically verified this 
for n = 2, 3 and d = 3,4. In Figure 12.11 we have illustrated the minimization of (|2.12p 
for a particular set of pure states (see caption). 

In conclusion, we have presented strong upper and lower bounds for the generalised 
and random Schmidt robustness for pure states. The problem of finding exact values is 
very hard, as in the end, one has to come up with an explicit convex decomposition in 
Schmidt rank n states on the one hand, and on the other with a construction of optimal 
Schmidt witnesses. We hope that our results may stimulate further work, especially in 
proving or disproving Conjecture 12.4.11 




Figure 2.1: Lower bounds on R r 2(ip) as a function of a\, where = ai|00) +a2(|ll) + 
1 22)). The graph was obtained using the witnesses from Conjecture 12.4.21 



CHAPTER 3 



Entanglement distillation: 
Introduction and classical results 



Distillation is the process of converting with some finite probability a large number of 
mixed entangled states into a smaller number of maximally entangled pure states. It 
relies on local collective manipulation of the systems together with classical communi- 
cation between the parties (LOCC). It has become apparent that characterizing and 
quantifying distillation is important in understanding the nature of entanglement from 
a physical point of view. Considerable effort has been devoted to the characterisation 
of distillable states, but even for bipartite systems the matter is not settled. In the 
next two chapters we analyse entanglement distillation in the bipartite scenario in great 
detail. 

This chapter is meant as an introduction, and most of the material presented can 
be considered as literature review. In the first two chapters entanglement was studied 
as static material, however entanglement distillation involves manipulation of quantum 
states, so that we start this chapter with a study of quantum operations (Section l3.ip . 
In the following section, this is then applied to the study of possible transformations 
between pure states both in the finite and asymptotic regime, with particular emphasis 
on entanglement concentration. In Section \3. 31 we review the most important measures 
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of entanglement and some of their properties. These divide natural into two classes, 
the geometrical entanglement measures, such as robustness measures already encoun- 
tered in the previous chapter, and operational measures. An example of the latter is 
the entanglement of distillation, which will be the main focus of our attention. We 
also show that in general, entanglement is not conserved in the asymptotic regime for 
mixed states. We discuss two examples of such irreversibility: the maximally correlated 
states and the bound entangled states. We end this chapter with a sketch of the most 
important distillation protocols for mixed states in the asymptotic regime. 



3.1 Quantum operations 
3.1.1 General quantum operations 

A quantum operation A is a physical operation which transforms a state p on Ti\ 
into another state p' on Ti.2- The aim of this sect ion i s to give a unified mathematical 



77 



1781 . 



2151 ]) of any quantum operation, 



representation (the so-called Kraus formalism 
independent of how this operation is to be implemented physically. Basically, there are 



two approaches to this 



138, 



214l | . The first is the axiomatic one: we impose on the 
class of maps A : £(TCi) — * £(^2) some physically motivated constraints and define the 
remaining operations as quantum operations. The second approach is constructive and 
the set of quantum operations is defined as the one which can be obtained by combining 
operations from a certain set of elementary operations. Both approaches turn out to 
be equivalent; we will sketch the proof of the equivalence here. 

The standard quantum mechanical formalism allows us to carry out four elementary 
operations on quantum states: 

(01) Unitary transformations: p — > UpW, with U a unitary operator. 

(02) Adding an uncorrelated ancilla: p — > p (8> cr, with a a density operator. 

(03) Tracing out part of the system: p — ► Tr^'(p), with Ji\ =H! ® H" . 

(04) Projective measurements with possibly postselecting: p — ► Y2i=i PipPii with Pi, 
pairwise orthogonal projectors such that YJi=\ Pi = ^ with k < re. 
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In general a quantum operation will only be successful with a certain probability 
p = Tr(A(p)), so that p' = A(p)/p, and A does not increase the trace. It is clear that 
the operations (01)-(03) can be represented as trace preserving completely positive 
(CP) maps. As the composition of CP maps is CP, we see that any quantum operation 
composed of operations (01)-(03) can be represented by a trace preserving completely 
positive map, these operations will also be called non-measuring operations. The reverse 
turns out to be true: any trace preserving completely positive map can be implemented 
by means of the quant um operat i ons (Q1)-(Q3). This result follows from the Stinespring 



dilation theorem 



77 



138, 



174 



276l | . Now we invoke the Choi representation of CP 
maps (in this context also called the Kraus representation), which says that any CP 
map A can be written as A(p) = LipL\ for all p € (see Theorem II .2.3j) . The 

operators Li are commonly called the Kraus operators. For a trace preserving CP we 
have that Yli^i^i = 1- If we allow operations of the type (04), we have to relax this 
constraint to LiLi < II and renormalise the density operator after application of 
A. Another way is to keep the constraint ^ i L\Li = 1 but allow postselecting on the 



particular outcomes of i. Thus we obtain (see for instance Ref. [771]): 



Theorem 3.1.1. A quantum operation A can be decomposed into operations of the form 
(01)- (04) if and only if A acts as a completely positive trace non-increasing operator: 

Hp) = ZiLipLl with Zi4 L i < £ - 

From a physical point of view, the requirement that a quantum operation is CP 
and trace non-increasing, is obvious. Thus we could equally have started this section 
by demanding this, and we would also have arrived to the elementary operations (01)- 
(04). 

3.1.2 Bipartite quantum operations 

Consider now a bipartite system on Ha ^T~Cb- We will review some important classes 
of local qu antu m operations where increasing degrees of communication are allowed 



22 



58, 



77 



242] 



<0> Local operations. This is the class of operations generated by the operators 
which have Kraus operators of the form A^(g)l and l(g>-Bj with ^ i A\Ai = Y2i B\Bi = 1. 
In this restricted class of operations, the parties are not allowed to communicate and the 
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operations are non-measuring. Combining them yields for an arbitrary local operation 

1-local operations. Suppose we allow local operations with one-way communi- 
cation from A to B. So let us suppose Alice performs a generalised measurement on her 
subsystem, with Kraus operators Ai with Yli^l^-i = ^- ^ the result of her operation 
is i the operation on the state acts as (notice that this is a trace decreasing operation 
as it will occur only with a certain probability) 

Kf(p) = {A l ®t)p{A\®t). 

She will communicate Bob that she found result i, and depending on that outcome 
he will perform a certain trace preserving operation defined by Bji operators, with 
^ ■ Bj { Bji = 1. Here the index i denotes that the operation Bob implements depends 
on the result that he got from Alice. Thus the state will change as 

Now if this same procedure is done on many particles the total ensemble will change as 

^ AB {p) = Y,{A i ®B ji )p{4®B\ i ). 

Of course postselection (for certain i) can occur in one or more terms of this expression. 

<^> 2-local operations. This class is also known as the class of local operations 
with two-way classical communication (LOCC). This situation is somewhat more com- 
plicated, since in this case, it is useful to do alternating measurements and communica- 
tions. Following a similar strategy as for 1-local operations, we focus on one particular 
outcome of each measurement, and do the summation at the end. Alice starts off with 
her measurements, and finds il as result, therefore the relevant operator is An <g> 1. She 
communicates her result to Bob and he will then perform 1® Bji(il) and communicate 
his result to Alice who will perform Ai2(il,jl) <£> 1. All these operators satisfy similar 
normalization properties as in the previous paragraph. In that way we get at the end 
of the series (for the total ensemble) 

Hp) = J^(An ® B m )p{A m ® B m )\ 

m 
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with m = (il, i2, . . . , in,jl,j2, . . . ,jn) and 

A m = A in (il,i2, ... ,i(n- 1), jl, j2, . . . , j(n - 1)) . . . A i2 (il, jl)An, 
B m = B jn (il,i2,.. . ,in,jl,j2,.. . ,j(n - 1)) . . . B j2 (il, %2, jl)Bji(il). 

Here also, postselecting can be done for particular choices of il, . . . ,in,jl, . . . ,j(n — 1). 

(} Separable operations. Separable operations are defined as any operation which 
can be written as 

Hp) = ^2(Ai ® Bi)p(Ai ® Bi)\ 

i 

where each A{ and Bi are arbitrary operations (including measurements). The normal- 
ization condition is ^2i(Ai ® Bj) '{A4 <8> Bi) = 1. It is obvious that any LOCC is also 



221 ] . However, we do have the following 



separable, but the reverse is generally not true 
theorem [58]: 

Theorem 3.1.2. It is always possible to simulate a separable operation using only 
LOCC, but with probability possibly smaller than 1. 



Proof. So, [57| suppose we have two sets of operators {A{\ and {B{\ such that Y^i(A%® 
Bi)\Ai ® Bj) = 1 and a state p. In general we will have that ^2nA\Ai > 1 and 
Y^iBjBi > 1. We renormalise A'j = aAi and B'i = aBi so that Y^i^\A\ < 1 and 
B'\B'i < 1. Then we let Alice and Bob perform the generalized measurements 
corresponding to these operators. They will obtain two measuring outcomes k\ and k 2 
and communicate their result to each other. Their strategy is to keep the state only if 
the outcomes of their measurement coincide, if k\ = k 2 . Forgetting which outcome was 
obtained, we see that p transforms as 

p^N^Pi = 2_^\Ai ® Bi)p(Ai ® Bi) T , (3.1) 

i Pi i 

with pi = Tr[(A'i (g) B'i)p(A'i (g> B'^y] and TV an appropriate normalisation factor. □ 
<0> PPT preserving operations. Literall y, th ese are operations which map PPT 



states to PPT states. Or with Rains' definition 
map 

A : p ^ A( P Tb ) Tb 



2411 ] operations A such that the induced 
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is co mpletely positive. All separable operations are PPT preserving, as follows from 
138 ] 

[{A®B)p{A®B)^ TB = {{A® (B^f)p TB (A^ B t )]. 

PPT preserving operations are the first class of operations which can be regarded as 
non-local, as for instance creation of an entangled PPT state is a PPT preserving op- 
eration. The key result on PPT preserving operations is that they can be implemented 
using LOCC when both parties share a specific entangled PPT state. This was first 
realised in the seminal paper by Cirac et. al. [5a]. They characterised completely gen- 
eral quantum operations (and in particular entangling operations) on bipartite systems 
by means of a generalised Jamiolkowski isomorphism. Because of its importance, we 
briefly sketch the isomorphism in the next section. 

3.1.3 Entangling operations 

The Jamiolkowski isomorphism (Section II. 2. 2|) is an isomorphism between a linear map 
from an input space to an output space and an operator defined over the tensor product 
of these two spaces. As we have illustrated, this correspondence is useful when the map 
acts on one part of a bipartite system. For maps acting on compound systems it is not 
so easy to interpret this duality. We now introduce an extension of the Jamiolkowski 
isomorphism with a very nice physical interpretation [58] . 

Suppose we have two systems A and B, each consisting of two d-level subsystems 
A = Ai, A2 and B = B\,Bi respectively. We will establish an isomorphism between a 
CP map A = A^Bj : C(H.a 1 b 1 ) - > £(J~(-A 1 B 1 ) and an operator D = Da 1 a 2 ,b 1 b 2 acting 
on the total system Ha ®7~Lb- The isomorphism is given by 

D = 1 A2 b 2 ® A AiBi (P Ai a 2 8> Pb 1 b 2 ), (3.2) 

where A is understood to act only on the systems Ai and B\. Here Pa\A 2 is the 
projector on the maximally entangled state \i(j}a 1 a 2 = ■^\H)a 1 A 2 (and analogously for 
PbiB 2 )- Equivalently we can write 

A(pAiBi) = d 2 TV A 2 B 2 {D P T A 2 B 2 ) = d 4 TrA 2 A 3 B 2 B 3 (DpA 3 B 3 PA 2 A a PB 2 B 3 )- 

From (13. 2|) it follows that D is the result of the action of A on two systems A\ and 
B\ which are prepared in a maximally entangled state with two ancillary systems. 
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The second form has an equally simple interpretation: if both parties share a state 
D, then they can implement the map A on a certain state pa 1 b 1 by simultaneous 
projecting pa 3 b 3 and D on the maximally entangled states of A2A3 and B2B3. From 
this isomorphism one can easily deduce the following correspondences: 

1. A is a separable operation if D is separable, and conversely A is an entangling 
operation if D is entangled (with respect to A and B). 

2. A is a PPT preserving operation if and only if D Ta > 0. Thus any PPT preserving 
operation can be implemented locally with the aid of a shared PPT state. 

This isomorphism has a large number of applications (see the original reference and 



8l|) which would lead us too far to go into here. A simple application is the classi- 



ficat i on of glo bal unitaries according to their entangling power on composite systems 



314 



336 



3371 ] . Entangling power of a unitary is defined as the average entanglement 



a unitary creates on product states. When A is a unitary map, D is a pure state and 
there exists a closed expression of the entangling power in terms of the entanglement of 
this state vector. Using this, we were able [641 ] to show that the unitaries with highest 
entangling power are special kinds of permutations, with possible exception for d = 6. 
The isomorphism then immediately yields a construction of generic examples of 4-qudit 
maximally entangled states for all dimensions except for 2 and 6. 



3.2 Pure states and the reversibility of entanglement 

In this section we will study the distillation of entanglement of pure bipartite quantum 
states, which is often referred to as entanglement concentration. An arbitrary pure 
state if) can be represented by a vector \ip) in Ha ® Hb = C d ® C d . Without loss 
of generality we have taken the dimensions of the two Hilbert spaces to be equal. A 
way of quantifying the entanglement of ip is by looking at the local density operator 
PA = Trg(|?/))(?/)|). For a separable state = \iP)a ® \'<P)bi the reduced density oper- 
ator equals pa = \iP)a(iP\, and is thus a pure state. The reduced density operator of 
an entangled state will have a rank larger than one. If the reduced density matrix is 
proportional to the identity matrix the state is maximally entangled, as is the case for 
the EPR singlet. For such states there is no information available locally, although the 
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global state is known exactly [26CJ, 12951 ] . A way of quantifying the purity of a local den- 
sity operator pa is by calculating its von Neumann entropy E(ip) = — Tv(pa log 2 pa)- 
We will call this quantity the entropy of entanglement or simply the entanglement of ip. 
Note that the entanglement is a fully additive measure: E{ip\ ® 1P2) = E(tpi) + E(ip2)- 
Writing ifi in the Schmidt basis 

d 

iv>> = E A * 



i=l 



with Aj > Aj + i > and Yli = 1> the reduced density matrix can be written as 
PA = 2i^!I*)a(*| an d the entanglement of ip as i£(^>) = — S« -^i 1°§2 -^f • Maximally 
entangled states in a d®d system have Schmidt coefficients \ = X/\fd and entanglement 
E = log 2 d. 

Entanglement concentration is the process of converting partially entangled states 
to maximally entangled states by LOCC. The reverse process is called entanglement 
dilution. We will show how in the asymptotic regime n arbitrary states ip can be 
concentrated into nE(ip) singlets for n — * oo. Conversely, nE(ip) singlets can be diluted 
into n states tp for n — ► oo. From the above it follows that the entropy of entanglement 
quantifies entanglement of pure states in the most natural way since it is preserved in 
the asymptotic regime: a collection of n states if) can be converted into m n states (p with 
LOCC and the ratio m n /n approaching E{ip) / E(<j)) for n — > oo. The entanglement of a 
singlet is given by E = 1 and serves thus as a proper unit of entanglement, the so-called 
ebit Q. 

In the first part of this section we discuss some elementary entanglement concentra- 
tion protocols, and present a dilution protocol. Next, we list some celebrated theorems 
of general pure state transformation in the finite and asymptotic regime. 

3.2.1 Pure state distillation protocols 
Local filtering — the procrustean method 



The procrustea: 



arf 



tion protocol 



method is conceptually maybe the simplest entanglement concentra- 



19]. It works even when both parties only share one single copy of the 



Procrustes is a giant from Greek mythology. Allegedly he would capture people and tie them to 
an iron bed. He would then hack off their legs or stretch them to make them fit the bed. 
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state. So let Alice and Bob share the state \ip) = cos 0|OO) + sin 6\11) with cos 8 > sin 9, 
then we can define 

sin 9 , . , . , . , . 

4~oo + ii 

cos 9 



a x = y i - a\a . 

Now the protocol consists of Alice performing the measurement corresponding to these 
operators and communicating the result to Bob. The protocol will be successful when 
Alice obtains the outcome corresponding to A$, in this case, occurring with probability 
2 sin 2 9, the output state will be a maximally entangled singlet. This protocol is optimal 
for single qubits and we will later see an extension of this protocol which is optimal for 
any pure bipartite state. Such a one-sided generalised measurement, where the state is 
conditionally kept will be called filtering. This is a technique which will be useful for 
mixed states as well. 

An analogous protocol works for entanglement dilution: if Alice and Bob start out 
with a maximally entangled singlet and Alice applies the operation defined by 

Aq = cos0|O)(O| +sin0|l)(l|, 

A x = sin0|O)(O| +cos0|l)(l|, 

then they obtain the state \ip) with probability one. 

Schmidt projection 

Following we illustrate this method for qubits. As the method does not work for one 
copy we suppose that Alice and Bob share n copies of the state \ip) = cos 0|OO)+sin 0|11), 
so that the combined system is described by the state vector 

= \^)® n = ®" =1 [cos0|OO)i + sin0|ll);]. 

The binomial expansion of this will have 2 n terms of which n + 1 will have different 
coefficients and the n + 1 corresponding subspaces will be orthogonal. Let us label 
these subspaces by k = 0, • • • , n corresponding to the power in which sin0 appears. 
Now Alice can perform a von Neumann measurement corresponding to these subspaces. 
The probability of projecting upon the fc-th subspace is given by 

Pk=(fj cos 2 ^ fc )0sin 2fe 0. 
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She can then communicate her result of the measurement to Bob or Bob can perform a 
similar measurement (which will have the same outcome). In either case they will end 
up with a maximally entangled state in a known 2 (^) -dimensional space. After distilling 
many such states, they can be converted to singlets without loss of entanglement. In 
the asymptotic regime this protocol can be shown to be optimal: the ratio of singlets 
that can be distilled per input state tp approaches E(ip) for the number of input states 
going to infinity. To achieve this, the Schmidt projection protocol has to be applied N 
times on batches of n states, with both TV and n going to infinity. 

Entanglement swapping 



334 



3381 ] was first studied in the context 



The phenomenon of entanglement swapping 
of coincidence count rates of entangled pairs of photons. In the experimental setup, 
two different sources emit pairs of entangled photons and two halves of the pairs are 
being measured simultaneously. The remaining halves then become entangled. Thus 
entanglement swapping is t he ph enomenon of two particles getting entangled without 

for an experiment d emo nstrating this effect). This 



29 



3351 ] . but has to our knowledge 



direct interaction (see Ref. 

phenomenon also occurs in a multipartite setting 
only been demonstrated for qubits. 

The idea of using entanglement swapping as an entanglement concentration proce- 
dure was first introduced in Ref. .30]. Consider the following two states of photons 



\cf)) 12 = cos<9|00) +sin0|ll), 
|0)34 = cos0|OO) +sin0|ll). 

Suppose |0)i2 is shared by both parties Alice and Bob, and in addition Bob has the 
state 10)34. When Bob measures particles 2 and 3 in the Bell basis {|00) + |11), 1 00) — 
|11), 1 01) + 1 10) , 1 01) — 1 10) j- , the resulting state for particles 1 and 4 will be one of the 
four states 

1^)14 = 1 =(cos 2 0|OQ) ±sin 2 6>|ll)), 

Vcos 4 9 + sin 4 9 

1^)14 = ^(101)1110)), 
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the first two occurring with probability (cos 4 9 + sin 4 9)/2 and the last two with proba- 
bility sin 2 9 cos 2 9. In the first two cases, the entanglement is lower than in the original 
pair, but in the last two cases Alice and Bob end up with a maximally entangled state. 
In any case, however, the procedure can be repeated. Repeating the whole protocol 
over and over again, yields on one hand a group of maximally entangled pairs, and 
on the other, completely separable states. The fraction of maximally entangled pairs 
achievable in this way turns out to be 2cos 2 #. This number just equals the maximal 
probability of the initial state to be converted into a maximally entangled state. In 
this sense, the entanglement swapping protocol is optimal. Note that since the second 
pair belongs to Bob, he has complete control over it, and he could in principal use 
any two-qubit state. It turns out that when he uses a state which is less entangled, 
the protocol does not concentrate entanglement but destroys it. When the state he 
uses is more entangled, no entanglement is destroyed, but the concentration is less ef- 



ets disti 



158 



lable from 



2641 ] which do 



ficient. (Here the entanglement is measured as the amount of sing 
the state). Extensions of this protocol have been proposed in Ref. 
not rely on such an iterative procedure since the second step of their protocol is just 
the procrustean method. 

With this we conclude our discussion of elementary pure state entanglement con- 
centration protocols. Various other protocols have bee n pr oposed, as for instance in 
Ref. [la ]. Noteworthy is the protocol outlined in Ref. [221! ] which demonstrates that 
entanglement concentration is possible using the effects of quantum statistics of indis- 
tinguishable particles only. 

3.2.2 Finite regime 

In this section we will discuss the more general question of feasibility of arbitrary 
transformatio ns b etween pure states. The first work on this subject was done by Lo 
and Popescu 1901 ] . in particularly they proved the following: 



Theorem 3.2.1. The Schmidt number of a pure bipartite state cannot increase under 
local quantum operations and classical communication. Entanglement manipulations of 
pure bipartite states with only one way communication are as powerful as those with two 
way communication, but they are more powerful than those with no communication. In 
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other words, we can replace any LOCC protocol on pure states, no matter how complex, 
by a single local generalised measurement. 

The groundbrea king work on pure state entanglement transformation was carried 



out by Vidal in Ref. 



299l | , where he introduced the notion of an entanglement monotone . 



His work was sign ificantly extended by several authors in a series of papers 



17(1 



213 



301 



116 



169 



3061 ] . Recently it became apparent that their proofs can be si mpli fied 



2951 ] to 



by making use of the quantum steering theorem as recognized by Verstraete 
which we refer the reader for full proofs. An entanglement monotone \x is defined as 
a quantity associated with a quantum state which on average does not increase under 
local transformations. That is, if an LOCC maps a state p with pro bability pi into 



states pi then p,(p) > YliPi^iPi)- The following theorem by Vidal 



2991 ] is crucial: 



Theorem 3.2.2. The restriction of an entanglement monotone \i to pure states tp is 
given by a unitarily invariant concave function of the local density operator fi(\ip)) = 

f{9 B )- 

This theorem implies that any entanglement monotone acting on pure states 

d 

iv>> = v 7 -^) 

i=l 

can be expressed as a concave function of its Schmidt coefficientay Aj. An important 
class of entanglement monotones Ef~, for k = 2, . . . , n is given by 

d fc-1 
i=k i=l 

where as usual Aj > Aj+i. The concavity of this function can 
variational characterization of eigenvalues of a Hermitian matrix 

Clearly we have that if we can convert <j> to ij) with certainty then we need to have 
E(ip) < E(<p) for every entanglement monotone E, a nd in particular this has to be the 



3C 



proven using the 



126]. 



case for E = Ej.. A surprising result by Nielsen 



213] is that the converse is true: 



Theorem 3.2.3. A necessary and sufficient condition for the transformation <j) to tp 
to be possible with certainty is that E^{ip) < E^((j)) for k = 2, . . . , n. Alternatively, in 



2 For convenience Schmidt coefficients in this section denote the 'normal' Schmidt coefficients 
squared. 
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terms of majorization, this can be expressed as Xj, y X^, where X^ stands for the vector 
consisting of the Schmidt coefficients of if) in decreasing order. 

For anyone who is familiar with unitarily invariant norms, this result will be rem- 
iniscent of Ky Fan's dominance theorem for unitarily inva riant norms and the dual 
representation in terms of sy mmetric gaug e functions 



much stronger can be proven 



116 



169 



3011 ] : 



25 



1271 ] . In fact, something 



Theorem 3.2.4. The probabilistic transformation from a pure bipartite state (f) into 
the ensemble {pi,4>i} can be accomplished using LOCC if and only if PiX^ i >z X^ or 
equivalently J2iPi^k(ipi) < Ek((j)) f or a ^ k. 

Note that from the theorem follows that the maximum probability p of transforming 
(ft into if) can be written as p = min^ or >- \^ The non-trivial part of these 

theorems is to show that there is an actual protocol achieving the transformation. The 
reader is referred to 295] for an elementary proof. It is there shown how the protocol 
achieving the transformation can be built from a few simple generalised measurements. 

Let us now return to entanglement concentration. The above theorem gives us 
the optimal probability of converting a pure state to a maximally entangled state of 
arbitrary dimension. When we allow the concentration of entanglement into different 
maximally entangled states (of different dimensions), we can ask how much entangle- 



ment can be distilled on average. This is answered in the following theorem 



lift 



1691 ]: 



Theorem 3.2.5. The optimal entanglement concentration procedure for a single pure 
bipartite state ifi with Schmidt coefficients X± > . . . > > is one that produces a 
maximally entangled state 4>j of j < d levels with probability pj = j(Xj — Aj+i). The 
corresponding optimal average distilled entanglement is 

d 

{E) m ax = X]( A i ~ A i+l)i l0 §2J- 
k=l 

The proof of this theorem uses techniques from linear optimization, showing that 
there exist optimal pj compatible with Theorem 13.2.41 The protocol achieving these 
probabilities can be obtained by going through the proof of this theorem. However, in 
view of Theorem 13.2.11 there exists a single local generalised measurement which does 
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the job. The measurement operators corresponding to Alice's generalised measurement 
are 



4 = £ 



i)A(i\, 



i=l 



(3.3) 
(3.4) 



indeed it is easy to verify that Y2j=i AjAj < 1 and 



Aj®l\l/)) = y/Pj\<t>j). 



(3.5) 



This protocol, a generalised procrustean method, thus turns out to be o 
Another application of Theorem 13.2.41 is the idea of catalysis 17Q . 



ptimal. 



290 ], where a 



pure state is used to help converting one state into an other, such that the pure state 
is recovered at the end of the procedure. Mathematically this means that there exist 
states (f), ip and \ such that 



but 



Finally note that one can broaden the discussion about pure state transformations by 
including questions about the approximate transformation of one state into another, 



see Ref. 



306]. 



3.2.3 Asymptotic regime 



In the introduction we have a 
which we now state formally 



ready s 

3,y, 



tetched the main result in the asymptotic regime, 



236J: 



Theorem 3.2.6. In the asymptotic regime n arbitrary states ip can be concentrated into 
nE(ip) singlets for n — > oo. Conversely, nE(ip) singlets can be diluted into n states ip 
for n — > oo. Thus the entropy of entanglement is conserved in the asymptotic regime 
and hence quantifies entanglement of pure states in a unique way. 



169], or using the Schmidt 



One can prove the distillation rate using Theorem 13.2.51 
projection distillation protocol []j]]. The dilution rate can be proved using the t< 'Im- 
portation idea ljj. Both the dilution and distillation rate can also be proven using 
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Nielsen's theorem (Theorem I3.2.3p , The crucial part of the latter exploits the asymp- 
totic equipartition theorem which implies that \ip)® n can be approximated by 2 nE ^' 
terms in the Sch midt decomposition. 

In Ref. 1911 ] Lo and Popescu showed that pure state entanglement is a truly in- 
terconvertible, in the sense that the amount of classical communication needed in the 
asymptotic regime vanishes. This is trivial for the entanglement concentration, as the 
Schmidt projection method works without communication. The dilution step is the 
nontrivial part, and it can be shown that nE(ip) singlets can be diluted into n pairs 



tjj with 0(y/n) classical bits of communication. Recently other asymptotica 



protocols have been proposed, some of which are universal (see for instance 

The entanglement concentration protocols we studied in the previous subsection 
are all of probabilistic nature, with the notable exception of Nielsen's result. As a 
result, the protocols reaching the entanglement yield E(ip) in the asymptotic limit 
are all probabilistic in a finite number of copies and only become deterministic in 
the asymptotic limit. The following theorem gives an upper bound to the asymptotic 
yield of singlets we can distill from a pure s tate "d) if we only use protocols that are 



ly optim al 



119 



205J]). 



deterministic for any given number of copies 



211 



2121 ] 



Theorem 3.2.7. Let tp be a pure state with Schmidt decomposition = Yli=i VM\^)j 
then the maximum number of singlets that can deterministically be concentrated is given 
by Ed e t(ip) = — log 2 Ai. When tp is a two-qubit state this bound can be obtained by a 
series of two-pair manipulations. 

Note that E^et just like E is a fully additive measure. The upper bound follows 
from Nielsen's theorem. The surprising result here is the situation for qubits where 
two-pair manipulations are sufficient. This is in sharp contrast with the probabilistic 
concentration protocols which generally need to be applied on a very large number of 
pairs. The values of E^ et and E do not coincide, although both yields can be reached 
with probability one in the asymptotic limit. One can generalise this and ask for the 
distillable entanglement as a function of an error exponent, representi ng the ra te of 



decrease in failure probability as the number of copies tends to infinity 



na . 



1201]. 



3 Sometimes also called blind protocols as neither prior knowledge of the state nor any form of 
quantum tomography is required. 
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3.3 Mixed states and the irreversibility of entanglement 

As we have seen in the previous section, for pure states there are essentially a finite 
number of entanglement measures in the finite regime and a unique measure in the 
asymptotic regime. There are two obvious ways of quantifying mixed state entan- 
glement, operational and geometrical. Operational entanglement measures are directly 
related to the physical applicability of entanglement. Geometrical measures can loosely 
be described as those qua ntifying t he distance from a state to the set of separable states. 



As excellent reviews 



54 



1301 . 



230] have appeared recently in the literature, we do not 
attempt to give an extensive review of the subject here. In particular the geometrical 
entanglement measures will be treated in brief, while for the operational measures we 
focus mainly on the entanglement of distillation. For simplicity, entanglement measures 
here are defined simply as entanglement monotones and thus all presented measures 
vanish on separable states and are on average non-increasing under LOCC, unless oth- 
erwise stated. 

3.3.1 Geometric measures 

The fir st rigorou s study of entanglement measures based on distance was carried out 



in Ref. 



293, 



2941 ] . Intuitively, the quantity 



mm D(p\\a), 

er£6 

with D a suitable distance measure, can be regarded as a measure of entanglement. 
The most important measure of this ki nd arises when D is the quantum relative entropy 



D = 5(p||<r) = Tr(plog 2 p — /olog 2 a) [255l . l26ll . 12911 ] . The associated measure will be 



called the relative entropy of entanglement and denoted Er. Note that for pure states, 
Er reduces to the entropy of entanglement E. Although the relative entropy of entan- 
glement is a distance measure, it stretches much further as it has several operational 
interpretations (see als o ne xt section). 



A related measure (293j can be defined by using as a distance measure the Hilbert- 
Schmidt distance D = \\p — a\\2- Although this quantity h as re ceived considerable 
interest, it was shown not to be an entanglement monotone 2191 ]. We will call this 



measure the best separable approximation in the Hilbert-Schmidt norm analogous to 
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the best separable approximation BSA(p) in the trace norm 1821 ] . which can be shown 
to be a proper measure of entanglement. It is defined as 

BSA(p) = min{||p - p s \\ \p > p s , p s G S}. 

or equivalently 

BSA(p) = mm{X\p = (1 - A)p s + Xp e , p s £S,p e £V\ S}. 
The geometric measure of entanglement (GME) is historically the first ex amp le of 



26^ for 



an entanglement measure. It was originally introduced by Shimony in Ref. 
pure states only as 

GME(V) = ^min{|| ty) - \</>)\ || 2 ,0 G S} = 1 - max{|<<^>| 2 | 4> G S}. 

Thus the GME is proportional to the max imu m overlap with a separable state. This 
measure was extensively analysed in Ref. 3171 ] . including for multipartite and mixed 



states. The mixed state generalisation of the GME was done by the usual convex roof 
construction, that is as a minimisation over all convex expansions of p: 



GME(p) = min I ^^GME(^ 



Dual to the distance measures are the robustness measures. While distance mea- 
sures quantify the distance of a set to the set of separable states, robustness measures 
quantify to which extend entangled states remain entangled under mixing. Recall (Sec- 
tion 12. 4p that the robustness Rk relative to a set 1C is defined as 



Rk(p) = min < R 



P + Rpk G5i/OfcG/C)i? > 



l + R 



Taking 1C to be either S, {1} or T> we obtain respectively the robustness R s , the random 
robustness R r and the generalised robustness R g . No te that the random robustness R r 
was shown not to be an entanglement monotone |ht| . However, for a given dimension, 
the maximum random robustness over all entangled states gives rise to a lower bound 
on the volume of separable states 304 ] . 



All the previous entanglement measures are very hard to evaluate for a general 
mixed state. For a NPT entangled state p, one can envisage that the sum of the negative 
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Measure 


Definition 


Key references 


Relative entropy 


min{Tr(jolog 2 £)\a G S} 


[123, 231. 293, 294, 309J 


BSA in trace norm 


min{||p - p s \\ \p > p s ,p s G S} 


[151, 171- 182. 193. 3201 


BSA in HS-norm 


min{||p - p s || 2 |ps G S} 


[24, P92, 219, 293, 3231 


Robustness 


mm{R > 0\p + Rp s G S, p s G 5} 


[112,304] 


Gen. robustness 


min{i? > 0\p + Rp d G S, p d G V} 


[34, ii, ill, 274] 


Ran. robustness 


min{i? > /9 + Rl G 5} 


[32, 117. 3041 


GME (pure) 


l-max{|(0|^)| 2 |<^G«S} 


[31, 41, 266, 317, 318, 3191 


Log. negativity 


lo g2 IIP Ts ll 


[229, 305] 



Table 3.1: Geometrical measures of entanglement. 



eigenvalues, the so-called nega tivity M(p), quantifies to some extent the distance of p 



to the set of PPT states 



3421 ]. And indeed, it can be shown to be an entanglement 



mon otone. However, it is often more convenient to work with the logarithmic negativity 



3051 ] . defined as 



E N (p) =log 2 ||p TB ||. 

The reason is that En is fully additive by construction and also has considerable oper- 
ational meaning, some of which we will outline in the next section (see also 9|). 

With this, we conclude our roller coaster through the main geometric entanglement 
measures. Note that the y also can be approached in a more unified way using either 



base norm Ref. 



230 



3091 ] or entanglement witnesses Ref. 33]. 



3.3.2 Operational measures 

Having a unique measure for pure states in the asymptotic regime, it would be natural 
to ask if one can generalise this to mixed states. We have seen that for pure states, this 
measure quantifies both the amount of entanglement that can be distilled from them in 
the asymptotic regime and the amount of entanglement needed to construct the state 
in the asymptotic regime. Initially, the idea was to define two different measures for 
these quantities for mixed state: the distillable entanglement Ed and the entanglement 
cost Eq respectively. 
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A first attempt to the entanglement cost was made in Ref. 2l[ as the convex roof 
of the entanglement 

E F (p) = miny] pjE(\ipj)) 

i 

where the minimum is taken over all convex decomposition p = This 
quantity is now known as the entanglement of formation. The evaluation of Ep for a 
general mi xed s tate is very hard, and analytical results are only known for high symmet- 
ric states 



287, 



309[| (see also 



461 ] and reference therein). The notable exception is for 



2 qubit sy stems, where an exact analytic formula is known in terms of the concurrence 



326, 



327 ]. 



t is however not clear that Ep has the operational meaning we want it to have 
1221 ] . In particular Ep(p) quantifies the entanglement needed to construct p via a 
specific procedure. Namely, given a decomposition p = ^iPityi) one chooses with 
probability pi a particular \ipi). Then, one constructs tpi out of singlets (costing in 
the asymptotic regime E(ipi) ebits) and forgets the particular ipi chosen. In this way, 
we obtain one copy of p. Repeating this process over and over, then to construct p 
on average we need Ep(p) ebits. Other procedures might be more efficient, such as 
constructing more than one copy at the time. In particular it is easy to see that 

Ef{p)= lim F(P ] < Ep{p). 

n—>oo n 

This quantity Ep°(p) is known as the regularised entanglement of formation. 



In Ref. 



1221 ] it was rigorously proven that E F °(p) equals the asymptotic value of 



m/n such that m qubits can be converted into n copies of p. Th is as ymptotic yield, 
we call the entanglement cost and is formally defined as (see Ref. 2301 ]) 



E c (p) = inf < r 



lim finf||A(Pf 







where the second infimum is taken over all trace preserving LOCC A and P+ is the 
maximally entangled singlet state. We have used the trace norm in this definition, but 
some other measures work as well. Note that the equivalence E F °(p) = Ec(p) is highly 
non-trivial as Ep° gives rise to a protocol with a variable number of input singlets, 
while Eq works with a fixed number of input singlets. Now even more surprising, it is 
conjectured that the entanglement of formation is additive, in which case we would have 
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207, 



Ec = Ep. This additivity is o ne of the major open problems in quantum information 
theory (see Ref. h3 ; 
cost is very hard pi 
prove the additivity 



2681 ] and references therein). Calculating the entanglement 
cost is very hard pr oblem, and its value is only known so far in cases where one can 



140 



206, 



207 



265 



307|, 



3331 ] . Both the entanglement of formation 
and the entanglement cost can be proven to be convex 



n 



241 



2421 ]. quantifying 



We now move on to the entanglement of distillation [211 . 
how much pure state entanglement we can extract from a state p. Analogously to the 
entanglement cost we define 



E D (p) = sup<r 



lim [infllP? 

n— >oo L A 



Hp 



OS n \ 







(3.6) 



where the second supremum is taken over all trace preserving LOCC A and P+ is the 
maximally entangled singlet state. Note that just like in the definition of Ec, the 
maximisation is over all protocols where the output has a fixed n umb er of singlets. 

242I ] that allowing 



But similar to the entanglement cost, it has been proven rigorously 
stochastic protocols does give rise to the same value for Ed. To elaborate a bit more 
on this point, suppose we have a protocol that from k copies of a state p produces 
n singlets for half of the time and no singlets at all for the other half. However, the 
above definition only allows deterministic protocols, that is in our case we would like to 
replace this stochastic protocol by a deterministic protocol which produces exactly n/2 
singlets. The equivalence of both definitions says that in the asymptotic regime this 
can always be done. The proof of the equivalence exploits the law of large numbers, as 
when a stochastic protocol is repeated a large number of times, it can be approximated 
with a deterministic protocol on a larger input. 



78 



1351 ] 



From the above, it will be clear that Ed < Ec- The following theorem 
pushes this observation further, stating that under some reasonable conditions the 
entanglement cost and the entanglement of distillation are extremal measures. 

Theorem 3.3.1. Let E be a positive functional on the set of states T>(H.a (8) Wb) 
satisfying for all states p G T> 

1. Normalisation: For any d-level maximally entangled state Pf, E(P+) = \og 2 d. 

2. Weak LOCC monotonicity: For any trace preserving LOCC map A; E{p) > 
E(A(p)). 
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3. Continuity: Let (p n ) n£ N and (cr n )„ e N be any two sequences of states on (7^) n( =N 
and (7in)n€N respectively. If \\p n — cr n || — > then 

E{p n ) - E(a n ) Q 
1 + log 2 dimHn 

Then: 

(i) IfE is additive (E{p® n ) = nE(p) for allpeV) then 

E D (p) < E{p) < E c (p). 

(ii) If the regularisation E°°(p) = lirriyj^oo E(p® n )/n exists for all peV then 

E D {p) < E°°(p) < E c {p). 

In some sense this theorem is a weakened version of the asymptotic reversibility 
for pure states. The above theorem can be used to provide both lower bounds to the 
entanglement cost and upper bounds t o th e entanglement of distillation. The conditions 
for the latter can be slightly relaxed [135 ]. 



Theorem 3.3.2. Let E be a positive functional on the set of states V(Ti.A ® ~Hb) 
satisfying for all states p G T> 

1. Weak LOCC monotonicity: For any trace preserving LOCC map A; E{p) > 
E(A(p)). 

2. Partial sub-additivity: E(p® n ) < nE{p) . 

3. Continuity for the isotropic state ( QOP 



p(P d ,d) = -^—(l + (3 d P + ). 
d z + (id 

That is, if we have a sequence of isotropic states (p(Pd, d))d — ► P+ when d — > oo, 
then we require that 

lim 1. 

d^oo log 2 d 

Then E D (p) < E{p) for all p. 
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Thus for an entanglement measure to be an upper bound to the entanglement 
of distillation it is enough to be continuous in the neighborhood of the maximally 
entangled state. We proceed by giving some examples of entanglement measures that 
are strictly larger than Ed- The relative entropy o f entanglement Er was the first 



example of an upper bound to Ed 



135 



241 



293 



294f | , this follows relatively easily from 



the previous th eore m. Another early upper bound to Ed is given by the logarithmic 
negativity En 3051 ]. and again one easily verifies that it fulfills the conditions of the 



previous theorem. Note that both Er and En do not satisfy the requirements of 



2381 ] that the previous two measures can be 



Theorem 13.3.11 

Rains showed in a seminal paper 
combined as 

E Rains (p) = min (S(p\\a) + log 2 \\a TB \\) , 

providing an upper bound to the distillable entanglement. Note that here the minimum 
is taken over all states, not just the separable ones. The Rains bound gives an upper 
bound to the PPT distillable entanglement, that is the entanglement that can be dis- 
tilled via PPT preserving protocols. Since these protocols include LOCC operations, 
the Rains bound provides also an upper bound to Ed- Despite this, -ERains is probably 
the strongest known upper bound known on Ed- Note that when a = p in the defini- 
tion, we simply recover En- A weaker bound can also be obtained by minimising just 
over the PPT states, so that the last term vanishes. The resulti ng quanti ty is known 



the relative entropy of entanglement with respect to PPT states 



239 



2411 ] 



E 



PPT (p) = mmS(p\\a) 
<rev 



Measures that satisfy the requirements of Theorem 13.3.11 are the regularised version 
E^ of the entropy of entanglement Er and the regularisation of Eft PT , Eft PTo °. There 
is some evidence that -E^ams = -Er PTo °, based on the equality for W erner states, 
although all that is known is that for qubits Eft PPo ° < E/R, a ins . 1164 1 . 

Lower bounds to the entanglement of distillation are much harder to obtain. One 
way of obtaining a lower bound is by providing an explicit distillation protocol (see 
Section I3.4H and calculating the rate. The most important lower bound up to date is 
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given by the hashing inequality [2l|, |70j, 1 137T ] 



E D {p) > E D ^ (p) > max{0, S(p B ) - S(p), S(p A ) - S(p)}. (3.7) 

Here Ed^ is the one-way distillable entanglement, that is the entanglement of distilla- 
tion where classical communication is only allowed in one direction (for instance only 
from Alice to Bob). Note that this relation also has important implications on the 
theory of quantum channel capacities. For some states the value of E R coincides with 
this bound, and thus allows for exact calculation of Ed. Examples will be given in the 
next section. 

From our treatment, it is easy to see that (see also 
on entanglement ordering) 

TP ^ TP TP -> TP°° TP PPT -> rpPPToo \ p 
> Er, CjQ > E R , ti R > tL, R > tip 



230] for excellent treatments 



and 



EN, E r PT > -ERains > Ejj. 



As an example, in Figure 13.11 we have plotted kno wn values of entanglement measures 



for the Werner states (d = 3). The Werner states 138l . 13211 ] are states of the form 

1 



' d 2 + (3d' 

with — 1 < (3 < 1, which are entangled (3 < -1/d or -1 < Tr(pF) < 0. We will 
study this class of states extensively in the context of distillability in Section 14.31 see 
also Section 14.11 Because of their prominent symmetry the value of some important 
entanglement measures_can be c alculated exactly. In particular we have E R =E R PT 



and -ERains = E R PToc [3, 0, l310l ]. Despite this, an analytical expression for both Eq 
and Ed is still out of reach. It is not even known whether the distillable entanglement 
for all entangled Werner states is strictly positive (see Section I4.3[) ! In Section 13.41 we 
redraw the graph for d = 2 for which the situation is somewhat simplified and for which 
good lower bounds are known on Ed- 
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0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 



Tr(pF) 



Figure 3.1: Entanglement measures En, Ep, Er = E P . PT and ^Rains = E^ PToc as a 
function of the parameter Tr(pF) for the Werner states. 
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3.3.3 Irreversibility 
Bound entanglement 

For pure states, entanglement is reversible in the sense that Er)(ip) = Ec{ip) = E(tp) 
for all pure states tft. In this section we will show that this is not necessarily the case for 

ed s tates, 



13J. The 



mixed states. The first indication of this was the discovery of bound entangl 
that is entangled states from which no entanglement can be distilled at all 
class of bound e ntan gled states includes the set of PPT states. We present the direct 
proof from Ref. [138l |. 

Theorem 3.3.3. For any PPT state p the distillable entanglement is zero. PPT states 
are therefore bound entangled. 

Proof. As we have seen in Section f3.1l the set of PPT preserving protocols includes all 
LOCC transformations. Therefore LOCC cannot map a PPT state to an NPT state. 
Now we prove that PPT states are bounded away from the maximally entangled state. 
Let p be a d-level PPT state. We have 

Tr pP + = Tr p Ts P+ B =^Tt P Tb F, (3.8) 

now since p is PPT, p = p TB is also a legitimate state and thus 

TrpP+ = ~TrpF<^, (3.9) 

where the inequality follows because the swap operator F has eigenvalues ±1 only, and 
therefore has a mean smaller than one on any state. Thus whatever LOCC protocol 
we apply, PPT states can never achieve a higher overlap than ^ with the maximally 
entangled state. Note that -j can be achieved by a product state. □ 

We come back to the question of the existence of NPT bound entanglement in 
Chapter [H Now the above theorem would imply irreversibility if we could prove that 
Ec{p) > for PPT entangled states p. However likely this might seem, it is not 
straightforward to prove this formally. This is because, all we know for entangled states 
is that Ep(p) > 0, but in the asymptotic regime, it might be possible to con struc t PPT 



entangled states with an arbitrary small amount of entanglement. In Ref. 



3021 ] Vidal 
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and Cirac obtained the following lower bound on Eq (the proof is straightforward and 
we omit it here). 

Theorem 3.3.4. Let P denote the projector onto the range of a state p on Ha ®T~Lb- 
If 

(ef\P® N \ef) < a N , 
for all product vectors \ef) £ TC^ N (8> T~L% N , then Eq > — log 2 a. 



This lower bound is relatively easy to calculate for some states 



302, 



303], in par- 



ticular it was shown to b e pos itive for some bound entangled state. The breakthrough 



however came with Ref. 



3311 ] which showed that Ec{p) > for all entangled states. 



Note that this was the one crucial piece missing in the original definition of entangle- 
ment (see Section ll.2p . The proof is relatively elaborate and proceeds in two steps. 
First a functional G is derived which is positive on all entangled states and zero on all 
separable states. Although this quantity G is by no means easy to evaluate, it serves 
the purpose as it is a lower bound to Eq. 



Generic irreversibility: maximally correlated states 

Irreversibility in entanglement transformation is not limited to bound entangled states. 



23J, 



2411 ] which 



Consider for instance the class of maximally correlated states (MCS) 
can be written as 

P = ^2 a ij\ ii )(jj\- 

A necessary and sufficient condition for a state p to be equivalent to a MCS after 



a possible local unitary transformation is given in Rcf. 



125(]. The tag MCS reflects 



that when both parties perform the same measurement, they will obtain the same 



result. The 



239, 



lains bound for this state is relatively easy to calculate and one obtains 



241 



3291 ] 



£Rains(p) = S(PA) - S(p) = E R (p) = E P R PT (p) = E PPT °°(p), 



238, 



2401 ] that this 



and by virtue of an explicit distillation protocol it can be proven 
value equals the PPT distillable entanglement. Note that the relative entropy is additive 
when acting on maximally correlated states. In view of the hashing inequality we have: 
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Theorem 3.3.5. The distillable entanglement of a maximally correlated state p is given 
by 

D{p) = S{ PA )-S(p) = E R (p). 

Next we derive a sufficient condition for the additivity of the entanglement of forma- 
tion, which allows for an exact expression of the entangl ement cos t Eg of some MCS. 



142 



256, 



2671 ]. this is a map 



Central is the notion of an entanglement breaking map 
M : C(H B ) -> C(H B >) such that for all p £ C{H A ) ® C(H B ) the state 1 B tg> M(p) is 
separable. A map A4 is entanglement breaking if and only it can be written as 

M(a) =^Tr(A>yi, 

i 

where Xi > and Y2i Xi = 1 and {<7j} is some set of density operators. With this we 



have 



3071 ] the following theorem. 



Theorem 3.3.6. Let V be a subspace of TCa 1 ®H Bl - Define a map Ai : C(V) — > 
C(H Al ) by 

M{p v ) = Tr Bl p v 
for all pv G £>(V). If M is entanglement breaking, then 

E f (pv <8> a) = E f {p v ) + E f (a) 

for all pv € jC(V) C C(Ti.A 1 8> T~Lb x ) and all a G C{TLa 2 ® Hb 2 )- ^ n particular, for all 
pv G jC(V) we have Ec(pv) = Ef{pv)- 



140 



3Q7J, 



3121 ] . Any MCS 



Let us now apply this to the maximally correlated states 
is supported on the span of so let = ^ \^i}c\H)A 1 B 1 , then 

Ha ® A<(|V)<^|) = 1VB a (|^><V|) = X) I^X^ilc ® NXiUi. 

i 

In other words, the map A4 associated to V is entanglement breaking, and thus E B (p) = 
Ec(p) for all MCS p. The entanglement of formation can be calculated exactly for 
qubits, and comparing the entanglement cost with the entanglement of distillation, one 
obtains generic irreversibility for all two-qubit MCS. In view of these results it would be 
worthwhile investigating whether a closed expression for the entanglement of formation 
can be obtained for MCS of arbitrary dimension. There is no reason however to suspect 
that irreversibility in higher dimension would not be generic. 
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Reversibility and Distinguishability 



The previous section raises the question if there are any reversible mix ed states. Trivial 



91 



12 A 



270 



3121 ] of the form 



examples constitute the so-called quasi-pure states 

i 

A simple measurement of both parties in the {|i)} basis will result in the collapse 
p — ► <8> | 1 1 therefore Ec{p) = En(p) = ^2 i PiE(ipi). Note that the process of 

identifying the tag is reversible in the asymptotic regime, hence the reversibility. 

More general quasi-pure states take the form 

p = ^PiPi ® Wi){i>i\, 

i 

with {pi} a set of separable and orthogonal states which are locally distinguishable. 



This latter condition is neces sary as t 



be distinguished locally 



22 



141 



lere are separable orthogonal states that cannot 



3301 ] with probability one. It seems to be an open 



ques tion whether there are reversible truly mixed states which are not quasi-pure (see 



also 



lsdi). 



Distingui shab ility has far deeper implications on distillation. We start with an 



example 



48, 



10l| . Let {\Bij}} denote the set of the four Bell states ±(|00) ± |11)) and 



|(|01) ± 1 10>) , and define 



1 1 

Pn = \ B ij)( B i. 



i,j=0 



We focus on the case where n is even. The state p2 is separable and therefore pf m is 
also separable. The relative entropy for n = 2m is bounded by 

E r (p 2m )<S(p 2m \\pf m ) = 2m-2. 

This last step follows by direct ve rifica tion. Now the set of two copies of all four Bell 



states can be locally distinguished [3131 ]. so that by sacrificing two copies, one can distill 
2m — 2 ebits from p2m- Since the entanglement of distillation is bounded by the relative 
entropy we have E r (p2 m ) = Ec{p2m) = 2m — 2. Other exam ples o f using explicit LOCC 
distinguishing in a distillation process can be found in Ref. (lpj . 

The link between distinguishability and distillation can be form alised as follows. 



First we start by recalling the Holevo bound (see for instance 



2141 ]. Chapter 11 and 



3.3 Mixed states and the irreversibility of entanglement 



96 



12). Suppose Alice prepares a state p x with probability p x and hands it over to Bob. 
He performs a measurement, with measurement outcome y. The random variable X 
with values x is characterised by the entropy S(X) = — ^2 x Px ^°S2Px an d similarly Y 
is characterised by S(Y) = — YlyPy ^ &2Py The joint entropy is then characterised by 
S(X, Y) = — ^2 xy p(x, y) log 2 p(x, y). With these definitions the mutual information is 
defined by S(X : Y) = S(X) + S(Y) — S(X, Y) and quantifies how much information 
X and Y have in common. The accessible information I acc is the maximum mutual 
information over all possible measurement strategies of Bob. As an example, suppose 
that Alice prepares a state in d possible orthogonal states, then clearly Bob can always 
access this information completely, and therefore / acc = log 2 d. When Alice prepares 
non-orthogonal states, this value will be smaller. Now, the Holevo bound is an upper 
bound to the accessible information: 

S(X : Y) < I acc < S(p) - Y,P* S (P X ) < N 

X 

with p = ^2 x PxP x , and N = log 2 dim'H. In particular it follows that N qubits can 
be used for information transmission of at most ./V bits. Now we generalise this in the 



scenario that Charlie prepares a 
who can only perform LOCC 



j fpar t ite s tate and gives a half both to Alice and Bob, 



12 



103, 



1431 ] . Note that the local accessible information 



is bounded by the accessible information /acc CC — ^acc- A Holevo-like bound can be 
obtained in this way (we present a strongly simplified version): 

Theorem 3.3.7. The local accessible information of a quantum state p E C(TCa^>T~(-b) 
is bounded by 

lL° c CC (p)<N-E F (p)-E F (p), 

with N = log 2 dim'H. Here p is the residual state, the state left-over by Alice and Bob 
after a particular measurement process. The amount of entanglement Ep(p) that may 
be extracted during the process satisfies E 'distilled < Ep{p). 

Thus to some extent, the entanglement (both initial and final) and local accessibility 
are complementary quantities. In our example above, the local accessible information 
J a° CC (/°2m) = 2, while N = 4m and E F (p 2m ) = 2m, so that ^distilled < E F (p) < 2m-2 
in the process of information extraction. The bound is here saturated; as in the process 
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of information extraction, we distill 2m — 2 ebits and the extraction of locally accessible 
information is in this case the best distillation procedure (see above). The upper bound 
of local accessible information hence gives an upper bound on the e ntan glement of 
distillation that can be obtained through local distinguishing (see also 262]]). We will 



encounter distillation protocols based on local distinguishing in the next section. 

3.4 Some important distillation protocols 

In this section we review some important distillation protocols. We analyse them in 
some detail for bipartite qubit systems, and only point to the relevant literature for 
higher level systems as already in the bipartite qubit scenario the main principles of 
distillation can be captured. From a practical point of view, it is also probably the 
most interesting. Note however that we discuss distillation protocols from a theoretical 
point of view, with disregard for practical limitations and operational simplicity. The 
foundations of all current distillation protocols have been set out in two brilliant papers 
Ref. [20, 



211 ] which we will follow closely. It is somehow remarkable how relatively little 



progress has been made since these were originally published. 

3.4.1 Recursive protocols 

The recurrence protocol (IBM protocol) 

Being the first protocol 



20 



211 ] . the recurrence protocol is also the simplest. Alice and 
Bob start with a large pool of identical qubit states with fidelity Y = Tr(pP + ) > 1/2. In 
every round of the protocol two pairs p® 2 {Y) get consumed, and with some probability, 
a single pair p'(Y') is produced with Y 1 > Y. Doing this a large number of times, the 
original pool of states p(Y) will be transformed into a smaller pool of states p'(Y'). 
Repeating the protocol, this time on the new pool of states, they can increase the 
fidelity further, and so on (hence the name recurrence protocol). 

So let Alice and Bob take two identical pairs with fidelity Y . In the first step of the 
protocol, they will bring both pairs in the standard isotropic form (see also Chapter [2] 
and Section I4TT1) 
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States of this form are invariant under local bilateral unitaries p' = (U ®U*)p' {U (&U*)^ 
and it follows that we can bring any state into the isotropic form by performing a 
random local unitary (this is the so-called twirling operation) without changing its 
fidelity. 



The unitary XOR (or CNOT) [21J defined by 



C/xor = |00)(00| + |01)(01| + |11)<10| + |10)(11|, 

behaves just like a classical XOR, flipping the second bit, conditionally on the first. 
Now the second step of the recurrence protocol consists of both parties implementing 
the local unitary XOR on their two halves. For the sequel we will call the first pair 
the source pair and the second pair the target pair. This local unitary operation called 
BXOR (bilateral XOR) entangles both Alice's and Bob's local halves, while it leaves 
the global entanglement between Alice and Bob unaltered. 

Finally Alice and Bob measure the second pair in the basis {|0), |1)}. If the out- 
comes of their measurement coincide, it is not so hard to show that the fidelity of the 
first pair will change as 

yi 9 y 2 + (i-y) 2 

9Y 2 + 6Y(1 -Y) + 5(1 - Y) 2 ' 

Figure I3T21 illustrates this relation, and it can be seen that the fidelity strictly increases 
when Y > 1/2. The first pair will no longer be isotropic, so that for further rounds in 
the protocol, the first step (twirling) will have to be repeated. When the measurement 
outcomes do not coincide, the first pair has to be disregarded, and the protocol was 
unsuccessful. It is important to note at this point that the exchange of the measurement 
outcomes requires two-way communication. 

As remarked at the start of this section, this round of the protocol will have to be 
repeated many times, until a sufficiently large pool of states with fidelity Y' has been 
obtained. Then a new round can be started, using the new pool of states. In this way 
the recurrence protocol is able to transforms states with fidelity Y > 1/2 into states 
with fidelity arbitrary close to one. The probability of doing so however tends to zero 
as Y — ► 1, as in every round at least half of the states are disregarded. 
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Y 



Figure 3.2: One iteration in the recurrence protocol. The curved line is Y'(Y), while 
the staircase curve illustrates the evolution of the fidelity through repeated recurrence 
rounds, starting with Y = 0.65. The black dots denote attractive fix-points, while the 
black square denotes a repulsive fix-point. This figure is an adaption of a similar one 
published in Ref. [a]. 



3.4 Some important distillation protocols 



100 



The QPA protocol (Oxford protocol) 

The QPA protocol {§9, 



194 ] was introduced in the context of quantum cryptography 
(Quantum Privacy Amplification) and is an improved version of the recurrence protocol. 
It achieves a certain fidelity with higher probability and does not require twirling. In 
fact it works in an identical fashion as the recurrence protocol, except that the twirling 
operation is replaced by a rotation. For completeness, the protocol consists of the 
following steps: 

1. Denote \<f>) = ~7|(|0) — and define a unitary operator R = \4>}(0\ — 

Then Alice implements R on each of her qubits and Bob implements R) on his 
qubits. It can be verified that this amounts to ordering the coefficients pij (see 



below), such that poo > Poi > Pio > Pn [fi£| • 



2. The parties implement the unitary BXOR on their two halves. 

3. Alice and Bob measure the second pair in the basis {|0), |1)} and keep the first 
pair if and only if their measurement outcomes coincide. 

To see that this is a distillation protocol, let us write an arbitrary state in the Bell 
basis 

i l 

P= Pij\Bij)(Bij\ + X] Pij,kl\ B ij)( B kl\, 

i,j=0 i,jy^k,l 

with 

\Boo) = -^(|00) + |11)), |B 10 ) = -J=(|00> - |11)), 
\B 01 ) = -±=(|01) + |10», \B n ) = ^(|01) - |10». 

Here the first index will be called the phase index, and the second the shift index. 

In this way, we can describe one step in the protocol as the mapping {pij,Pij,ki} ~^ 
{p~ij,Pij,ki}- Now by going through the algebra it can verified that the coefficients 
only depend on p^j, more precisely: 

Poo + Pn ~ 2p iPio 

Poo — p , Pio — — p — , 

- Pl+Pio ~ 2pooPii 

Pm = p ) Pn = — p — ! 
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with P = (pnn 



be proven 
Poo -> !)• 



+ Pn) 2 + (Poi +Pio) 2 - It is not completely straightforward, but it can 



194i | that when poo = r ^ I (pP+) > 1/2 the protocol converges to P+ (that is 



Local permutations 



Sev e ral impro vements of the early recursive protocols have been proposed 



209, 



210 



1571 . 



201 



21a]. It has also been realised that the recurrence and the QPA protocol 



can be interpreted as coming from a quantum stabilizer code^ and how to convert an 
arbitrary quantum stabilizer code to a distillation protocol 



208 



315f | . As it turns 



out, all these different sche mes can to a large extend be captured in the so-called local 



permutations schemes 



68, 



155]. 



The main idea is to start with n pairs every round, and at the end measure m pairs. 
Such a protocol will be called an [n, m]-protocol, the QPA and recurrence protocols are 
thus [2, 1] protocols. Naturally, for n > 2 the XOR operation will have to be replaced 
by some multipartite unitary entangling operator. One way of doing this stems from 
analysing the BXOR operation in a bit more detail. Note that the XOR operation can 
be written as 

EWK)|j) = j), 
with © addition modulus 2. From this follows that the BXOR operation acts o n two 



Bell states as (note that as usual the two Bell states are shared by both parties) [311 ] 



U B xon\Bij)\Bki) — \BiQk t j)\Bkj^i). 

Thus although BXOR is an entangling operation, it happens to be non-entangling on 
Bell states 1. In view of the above we write ijkl B 2¥^ i'k'j'l' and interpreting ijkl as a 
bitstring, we can attach a number from 1 to 4 2 to it. It is then not so hard to show 
that in this sense the BXOR corresponds to a permutation of 4 2 elements. Generalising 
to n pairs, one can consider permutations of 4 n elements, such that the corresponding 
unitary is local. The aim of working with this restricted set of unitaries, rather than say, 
all local unitaries or even all permutations is to make optimisation easier. Furthermore, 



4 To see that BXOR, or thus XOR is entangling it is sufficient to consider an example such as 
t/ XO R.(|00) + |10)) = |00) + |ll). 
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an elegant closed expression for these special permutations is known [68(]. The problem 
of finding an optimal [n, m]-protocol is highly non-trivial, as a trade off has to be 
found between the increase in fidelity and the probability at which this happens. Also, 
in general, the choice of the permutation will be state dependent. Once a suitable 
permutation has been found, the distillation protocol is identical to the QPA protocol, 
except that in the second step the BXOR need to be replaced with the corresponding 
local unitary and in the third step, m pairs need to be measured. 



Note that recurrence protocols can be readily generalised to higher level states 



p wit 



i Tr( pP+) > 1/d by replacing the XOR by some multilevel entangling unitary 



28 



131 



202 ]. 



3.4.2 Hashing and breeding 

The hashing protocol 



21 



269l | . unlike the recurrence schemes has a strictly positive 
yield for Bell diagonal states with S(p) < 1. Somewhat contradictorily, the protocol 
does only need 1-way communication. In this section we present a simplified version, 
the breeding protocol [^(J, which assumes that Alice and Bob share initially S(p) pure 
singlets (we call them here pre-distilled singlets) per pair of p. However as it turns out, 
both protocols give rise to the same yield. 

In the breeding protocol Alice and Bob share an arbitrarily large amount n of 
identical Bell diagonal states: 

p ®n = ^ Pkx,h-k n ,ln P kxh ® • • • ® P k n l n , 

with P a b = \B a b)(B a b\. Note that any state can be made Bell diagonal by a certain 
discrete twirl operation (see Section [4. 1|) . The main idea of the breeding/hashing pro- 
tocols is now to interpret the above as Alice and Bob sharing P^h ® ■ ■ ■ <8> Pk n i n with 
probability Pfci,/ r --fe n ,z n - This way of reasoning is perfectly valid as the density operator 
contains all physical observable information, independently of the preparation proce- 
dure, and the state p® n can be thought of as prepared in this way. The goal of Alice 
and Bob is then to identify the specific sequence (k,l) = (ki,-- - , k n , l±, ■ ■ ■ l n ), after 
which Bob can simply apply local unitary transformations to obtain P^ 1 . 

The breeding protocol uses pre-distilled singlets as target states for BXOR opera- 
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tions to gain information on the source states characterised by the sequence (k,l). If 
Alice and Bob implement a BXOR with a pure target \Bo m ): 

U BX oB.\Bki)\B 0m ) = \Bk t i)\B 0t i 9m ), 

then the shift index of the target pair toggles when the source pair has shift index 
1 = 1, while the source pair remains unchanged. So when Alice and Bob perform BXOR 
operations on each pair of a subset of n' pairs with a pure singlet as a fixed target state, 
they can identify the number of source pairs with shift index 1, by measuring (locally) 
the shift index of the target state. Thus, this procedure allows them to identify the 
parity of the shift indices in the subset. Similarly, by replacing BXOR by a different 
local unitary, they can determine the parity of the phase indices in the subset. 

The question is now how many parity checks are needed to gain complete informa- 
tion of the sequence (k,l). This is important, as each parity check consumes 1 pair of 
their initially pre-distilled singlets. A total random sequence of n pairs needs 2n par- 
ity checks to uniquely identify, so in doing so we would lose more entanglement than 
we gain. However, for n — > oo the sequence will belong to the class of 2 nS ^ typical 
sequences with probability going to one. It can be shown that the probability that two 
different sequences have the same parities after r checks on random subsets is smaller 
than . So typically, after every parity check the list of typically sequences reduces by 
a factor ^ and thus r = nS(p) parity checks will be sufficient to identity the sequence in 
the asymptotic regime (see Ref. [20| for a rigorous analysis). The yield of the breeding 
protocol is therefore 

lim n ~ nS{p) =l-S(p), 

n^oo n 

which agrees with the hashing inequality (|3.7p . Both hashing and breeding can be 



generalised to multilevel systems [3 111 ] . As a side remark, note that there is a close cor- 
respondence between 1-way distillation protocols and quantum error correcting codes. 
The most important implication of this correspondence is the quantum noisy coding 

recall that this inequality gives a 



21 



2d, 



1371 . 



2691 ] . 



theorem, which follows from the hashing inequality 
lower bound on the 1-way distillable entanglement) 

In Figure [331 the value of distillable entanglement by the hashing protocol (denoted 
H) has been plotted for the Werner states. As can be seen from the figure, the hashing 
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Figure 3.3: Entanglement measures E F and E R = 



PPT 



ERa 



E 



PPToo 
R 



as a 



function of the parameter Tr(pF) for the Werner state s. H and E^ are the distillation 



rates of the hashing protocol and the protocol of Ref. 
obtained courtesy of E. Hostens). 



1561 ] respectively (the data was 
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protocol works only for Ti(pF) < —0.62. This can be bypassed by preceding the hashing 
protocol by a few recurrence steps, increasing the fidelity above this threshold. While 
the recurrence protocol exploits two-way communication, only the hashing protocol is a 
proper asymptotic protocol. By introd ucing two- way communication on the asymptotic 



1561 . 



level, small improvements can be made 
rate for the Werner states has been plotted in Figure 13.3 



3081 ] . The currently best known distillation 



CHAPTER 4 



Distillability, activation and NPT 
bound entanglement 



Having encountered PPT bound entangled states in the previous chapter, this chapter 
deals with the more general question of distillability of general quantum states. We 
call a state distillable if and o nly i f the disti llable entanglement is strictly larger than 



zero. In three seminal papers 



131 



133, 



134i | the Horodeccy derived the fundamentals 



of the theory of distillability. Namely they showed that (i) any two-qubit entangled 
state can be distilled; (ii) any general d-level state p is distillable if and only if there 
exists a two-qubit subspace on which p® n is entangled for some n; (iii) any state vi- 
olating the reduction criterion can be distilled. The second statement has a number 
of nontransparent reformulations. Here we present a unified presentation of these and 
other distillability theorems, using tools borrowed from the separability problem. This 
allows us to systematically construct positive maps, with the same properties as in (iii) 
of the reduction criterion [59(. We give several other applications, including a numerical 
method for checking distillability which we apply to estimate the volume of distillable 
states. (Sections 

In the next section we address t he questi on of NPT bound entanglement in detail. 
After reviewing previous results 



72, 



83 



3161 ]. we show that this question is equivalent 



106 
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to a specific Schmidt number problem and to a separability problem. In particular, we 
show it is equivalent to the existence of a certain class of PPT entangled states. We 
end this section with a set of one parameter states which we conjecture to exhibit all 
forms of distillability (the rainbow states) 601 ] . 



As bound entangled states are not directly useful for typical quantum information 
tasks such as teleportation, the question arises whether they have any non-classical 
use at all. We will see that bound entanglement can be activated. That is, bound 
entanglement can be released with the help of other ( poss ible bound) entangled states. 
A first activation-lil 
obtained in Ref. 



87 



i e ef fect was presented in Ref. but the main results were 

203| | . This is the topic of Section [Ol 



4.1 Classic results 

4.1.1 States symmetric under local unitaries 

In previous chapters we have encountered the isotropic states and the Werner states. 
These states are invariant under certain local unitaries, and because of their symmetry 
they have been exhaustively used in the study of distillability. In this section, we review 
briefly the general theory of states symmetric under local unitaries, as developed in 



Ref. 



309[]. 



Let Go be a subgroup of the unitaries, or possibly the whole group of the unitaries. 
Then we can define the group G of the unitaries consisting of all pairs of the form 
U ®U' acting on a Hilbert space TL = TLa 'Si'Hb, where U £ Go and U' is some given 
function of XI . The set of bipartite states left invariant by the group G will also be 
denoted as [/[/'-invariant instead of G-invariant. An arbitrary state can be projected 
onto an [/[/'-invariant state by twirling it 

T G (p)= [ dU(U ®U')p{U ®U')K 
Jg 

Here the integral is performed according to the Haar measure on Go- In particular, for 
a finite group, the integral can be replaced by the sum 

T g(p) = tL'£ i U®U'p(U®U')1 

1 1 UeG 
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The commutant G' of the group G is defined as the set of all operators invariant under 
G-twirling, or equivalently the set of all operators A such that [U, A] = for all U G G. 
It follows that the set of bipartite states left invariant by G is just the intersection of 
the state space with the commutant G' of the group G. Note that there always exist 
a finite basis of operators B(G') spanning the commutant G' . Now for any operator C 
it is straightforward to check that 

Tr{T G {p)C) = Tr(pT G (C)). 

This relation has a number of important consequences: 

1. A G-invariant state p is uniquely determined by the expectation values Tr (pA) 
for all A G G' . In particular, if B(G') is finite, then the set of G-invariant states 
can be characterised by a finite set of parameters gA = Tr(pA), for all A G B{G'). 

2. For an arbitrary state p, the parameters gA = Tr(pA), for all A G B(G'), are 
invariant under twirling. From the previous point it follows that in order to 
calculate the image of p under twirling it is sufficient to calculate these parameters. 

3. Suppose we need to calculate an extreme for some G-invariant state p such as 
mine Tr(/?G), where the operators G belong to some set C. It is then sufficient to 
minimise over the (smaller) set) Tq{C). 



For the Werner states 



1 and F. The isotropic states 



321 1 , which are [/[/-invariant, the commutant is spanned by 



131 



24=11 ] are [7[/*-invariant and a basis is given by 11 and 
P + . Note that for qubits, the Werner and isotropic states are related via a local unitary, 
while this does not hold for higher dimensions. As we have seen in the previous section, 
some distillation protocols make extensive use of twirling. As remarked in Section P2.4.1l 
for the Werner and isotropic states, it can be shown that the twirl can be written as a 



finite sum ([2l|, |83J). In practice however, twirling can be done simply by performing 
a single random UU or UU* transformation and forgetting which specific unitary was 
used. 

The Bell diagonal states encountered in the previous chapter are invariant under 
the group 

G = {1 (8) t, -a x (8) a x , -a y <g> a y , -a z (8) a z }, 
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with as usual a x ,a y , o~ z the Pauli matrices. A generalisation of Bell diago nal states and 



311 1 



their symmetry group for arbitrary dimensions is presented in Ref. 

From basic symmetry groups one can generate others using tensor products. Con- 
sider the case where we have two symmetry groups G = {U (g) U'} and K = {V (g> V'} 
acting respectively on Hi = TLa ® ~Hb and H.2 = H'a ® Ti! b- Let Bq and Bk be 
a basis for the UU'- and the yi/'-invariant states respectively. Then a basis for the 
t/t/'VV'-invariant states acting on 7i\ ®7i2 will be given by Bq®Bk- In what follows, 
we will number the systems belonging to party A with odd numbers and party B with 
even numbers. As an example, a basis for the i/i/VV-invariant states is given by the 
operators {I12 <8> 134,-^12 <8> 134,112 ® ^34,-^12 ® -^34}- Imposing the extra condition 
that the coefficients of Fyi <8> I34 and I12 ® .F34 should be the same, we end up with the 
so-called UUVV F-\xivai\axit states as introduced in Ref. [30^], where they were used as 
a counterexample to the additivity conjecture for the relative entropy of entanglement. 
We return to this set of states in Section 14.31 

4.1.2 Fundamental theorems 

In this section we present the thre e basic th e orem s (Theorems I4.1.R 14.1.21 and 14.1.3ft 



on distillability by the Horodeccy 



131 



We also present a first result 



133, 



1341 ]. as mentioned in the introduction. 



1311 ] on NPT bound entanglement. Because of their 



importance we present full proofs. 

Theorem 4.1.1. All entangled two-qubit states p are distillable (Eo(p) > 0). 

Proof. The proof consists of locally transforming p into a state p", which has a singlet 
fraction larger than a half. In Section 13.4.21 we have seen that for such states the 
distillable rate is strictly positive. 

Since p is an entangled two-qubit state it has an negative partial transposition 
p TB ^ 0. We use the same technique as in (]1.7p to transform the state such that 
it is detectable by a suitable entanglement witness. Let us denote the eigenvector 
corresponding to the negative eigenvalue of p Ts with \tp), so that (tp\p TB \ip) < 0. Now 
any vector can be written as = A <8> 1]^+) (see (II. 3ft ). With this we get 

Tr[(A f t)p(A ® t)Pl B ] < 0. 
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Now let Alice perform a filtering operation corresponding to the Kraus operators A. 
Then the state becomes 



, (A* <g> l)p(A ® 1) 
P ~ 



Tr[(At®l)p(A(8)l)]' 
In this way we have transformed the state such that Trp'F < 0, as P+ B = F. Now in 
£(C 2 (g) C 2 ) the swap operator can be written as F = 1 — 2P_, with P_ the projector 
upon the state |^_) = (|01) - |10))/\/2. So that we have 

TV p>P„>\ 

a simple local rotation can then map p' onto a state p" such that Tr p" ' P + > I. □ 
The next theorem is the fundamental theorem of distillation, we present a slightly 



different proof as compared to Ref. 



134]. 



Theorem 4.1.2. A state p onTC = T-La^'Hb is distillable if and only if there exist two 
dimensional projectors P : W® n — > C 2 and Q : 7i^ n — > C 2 such that for some n > 1 the 
state 

p" = (P ® Q)p® n (P ® Q)t (4.1) 

is entangled. Since the resulting state acts on C 2 (8> C 2 ; i/iis is equivalent to p" having 
a negative partial transposition. 

Proof. The condition is clearly sufficient, as when such projectors can be found, many 
copies of p" can be produced, and from the previous theorem it follows that these can 
be distilled. Suppose now that the state p is distillable, and thus Ejj(p) = r > 0. Let 
us fix arbitrarily small e' > 0. From the definition of distillable entanglement (13. 6p . it 
follows that for any e > there exists a number n and a trace preserving LOCC A such 
that 

\\K{p® n ) - Pf rn \\ < e, 

with P + the maximally entangled two-qubit singlet. This implies that, by choosing e 
small enough, and disregarding all but one approximate copy of P+, we can find an 
LOCC protocol A' such that 

\\A'{p® n )-P + \\ <e', 
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We write this LOCC protocol in the separable form (with N a normalisation factor) 
p< = A'(p® n ) = 1 ^(A, ® Bi)p® n (^ ® Bi) f - (4-2) 

i 

Since this is approximately P+ we know that p' is entangled and we can assume that all 
(Ai ® Bj) project onto the same two-qubit space. Thus every A{ and Bi can be written 
as 

Ai = |0}(Vu<| + |l)(0Ai| ; 5< = |0)(VBi| + |l)(0Bi|. 

Clearly, since p' is entangled, at least one term in this sum (14.2)) must be entangled, 
for instance for % = 1. We denote the projector on the space spanned by {<t>Ai,ipAi} 
as P and the projector on {</>bi,^>bi} as Q. We have thus shown that the following 
two-qubit state 

/ = {A x ® Si)p® n (Ax ® Bx)t = (Ax ® ® Q)p 0n (P ® Q)t(Ai 

is entangled. It follows that also the two-qubit state (P (8> Q)p® n (P (8> Q)t must be 
entangled. □ 

If the condition of the theorem is satisfied for a particular number n, then we call 

. This result has 



the state pseudo-n-copy distillable or in short n-distillable 



72, 



proven to be extremely useful and is often used as a defi nitio n of distillability. For 



completeness, we state the following corollary (see Ref. [72. Il3ll ]): 



Corollary 4.1.1. A state p on Ti = TCa ® 7~(-B is n-distillable if and only if there exist 
some local operators A and B mapping onto C 2 such that 

Tr(F(A ® B)p® n (A B) ] ) < 



or 



Tr(P + (A®B)p® n (A®B)^) 
Tr((A <g> B)p® n (A <g> B)t) > / ' 

with P+ and F acting on C 2 ® C 2 . 

Proof. This corollary follows directly from the previous theorem and the fact that any 
entangled two-qubit state violates both the partial transpose and reduction criterion. 
Via local filtering, see equation (|1.7j) and the subsequent discussion (see also the proof 
of Theorem I4.1.1jl , one can then transform these states such that they are detected by 
the witnesses F and II — 2P + . □ 
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Theorem 4.1.3. All states violating the reduction criterion are distillable. 

Proof. In a similar fashion as in the proof of the previous corollary, a state p violating 
the reduction criterion can be locally transformed into a state p' , which can be detected 
by the witness 1 — dP + , or 

Tr(p'P+) > \. 

As remarked in Section 13.4.14 for such states a recurrence protocol is available. To get 
a finite yield, one can then use a multilevel generalisation of the hashing protocol, see 

Section [3X3 
Another way 



138] of showing this is to C/[/*-twirl the state p', after which it will 
become isotropic. The fidelity remains invariant under UU* -twirling so that the state is 
still entangled. Now both parties apply on their half the operation P = |0)(0| + |1)(1|. 
This projects the state p' onto a two-qubit isotropic state p" with Tr(p"P + ) > ^, and 
such states are distillable, as we have seen in the discussion of distillation protocols 
(Section E3D. □ 

The second proof implies that all states violating the reduction criterion are distil- 
lable if and only if all entangled isotropic states are. A dual result can be obtained for 
the partial transpose criterion and the Werner states: 

Theorem 4.1.4. All entangled NPT states can be distilled if and only if all entangled 
Werner states can be distilled. 

Proof. Since entangled Werner states have a negative partial transposition, one direc- 
tion is trivial. To prove the other direction, note that as before, any NPT state p can 
be locally transformed into a state p' such that 

Tr(p'F) < 0. 

Now twirling this state, we obtain an entangled Werner state p" (because Ti(p'F) = 
Tr (p"F) < 0). If we are able to distill all Werner states, then also the original state is 
distillable. □ 



In Section 14.31 we will present strong evidence that some of the entangled Werner 
states cannot be distilled, and hence for the existence of NPT bound entanglement. 
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We end this section on the fundamental theorems with a remark on finite copy 
distillation. We have seen that any distillable state is (pseudo)-n-copy distillable. The 
prefix pseudo reflects that if we project upon a two dimensional subspace, we are only 
half way through our distillation process. Indeed, in the next step we would like to 
repeat this procedure m times on batches of n copies of p, giving us m copies of the 
qubit pair p'. Finally we can use existing protocols to extract maximally entangled 
singlets from p'® m for m — > oo. Thus in general, even when a state is n-distillable we 
might need an infinite number of copies to get a finite yield of distillation. 

It is then natural to ask whether there are any states which can be distilled using 
a finite number of copies, and in particular if there are any states which are single 
copy distillable. Clearly, pure states are single copy distillable and so are the quasi- 
pure st ates enco untered in Section [3.3.31 For generic qubit mixed states, it has been 
proven 



172 



189l | that this is not possible. For qubits the best strategy to incr ease the 



entanglement of formation is transforming the state into a Bell 
The explicit form of this transformation was found in Ref. 142 



di agon al state 
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2981 ] . See Ref. 



147] 



for bounds on single copy distillation of arbitrary level systems. Despite that in general 



single-copy distillation is not possib 
which they called quasi-distillation 



e, t he Horode ccy s howed the existence of an effect 



1361 ] (see also [2981 ] ). In this scenario, it is possible 



to extract a pure maximally entangled singlet from some (specially constructed) single 
mixed state, but with probability going to zero, as the fidelity of the output state 
approaches this maximally entangled singlet. 

The problem of finite copy distillation for more than one copy of the state has been 



solved for qubits in Ref. 



1681 ] . It was shown that n > 2 copies of a two-qubit mixed 



state p can be distilled into a pure entangled state if and only if the range of p cannot 
be spanned by product states. Surprisingly, this result is independent of n. This result 
is important as it shows that for entangled qubits, which are all 1-distillable, we need 
typically an infinite number of copies to distill pure entangled states. The reason is 
that although we can get a state with very high fidelity given a large number of copies 
(for instance by means of the recurrence protocol), in general this final state will still 
be mixed entangled. Only in the asymptotic limit of the number of copies, it is then 
possible to obtain a pure state. 
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4.2 Uniform treatment of distillability 



This section contains an original contribution first published in Ref. [59(. The motiva- 



tion behind this work was an attempt to answer the following questions: 

What is the special role the reduction criterion plays in the story of distil- 
lation? Are there any other positive maps with the same property? Can 
we find a decomposable map which is stronger than the reduction criterion, 
but so that all states detected by it are still distillable? 

4.2.1 Borrowing tools from the separability problem 

We have seen in Chapter [T] that separable states form a convex set, from which fol- 
lows that there exists an entanglement witness for each entangled state. Using the 
Jamiolkowski isomorphism between operators and maps, this can be translated in 
terms of positive maps. The idea is to apply a similar reasoning to the set of dis- 
tillable states. The crux is that 1-undistillable states also form a convex set. Indeed, 
from Theorem 14.1.21 and the linearity of the partial transpose, it follows that mixing 
1-undistillable states can never yield a distillable state. 

The question is now to determine the corresponding witnesses. Notice that from 



Theorem l4.1.2l it follows [83f] that distillability is equivalent to the existence of a Schmidt 
rank two state \ip) = C\\ai, b\) + C2|a2, 62), with {|ai), {0,2}} two orthonormal vectors in 
7Y® n and I&2}} two orthonormal vectors in H^ n , and some n such that 



Thus D = \ip)(ip\ TB , with \ip) Schmidt rank two, can be interpreted as a distillability 
witness. The following observation generalises this. 

Observation 4.2.1. A state p is 1-undistillable if and only if 



for all operators D = o B with a a Schmidt number two state. 

This observation is compatible with the fact that PPT states cannot be distilled 
as D is decomposable by construction. Note that to study the phenomenon of distilla- 




(4.3) 



Tr(Dp) > 
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tion it is sufficient to characterise the 1-distillable states; n-distillable states p can be 
characterised by looking at one copy of p® n . 

Now in a similar fashion as the separability problem and Schmidt number problem, 
we have to look at map A associated with the operator D as given by the Jamiolkowski 
isomorphism 

D = A)P+. 

Let us start with the simple case where D = |^)(^| Ts . As D is a decomposable 
operator, the map A is decomposable and of the form 

A = ToA CP , with |^)(^| = d(t ® A CP )P + . 

What more can we say about this completely positive map A CP 1 Recall that the 
general form of such a map is given by (see Theorem 11.2.31) 

A CP (A) = Y,V l AV l j , (4.4) 

i 

where Vi are arbitrary operators. Now let \ip) = Xa=i c i\ a ibi) be an arbitrary Schmidt 
rank k state, we obtain after some algebra that the associated map3 is given by 
A CP (A) = VAV\ with 

k 
i=l 

In our case k = 2, so that V will have rank 2. 

Let us now move to the case where D = o Tb , with a an arbitrary Schmidt number 
two state. A similar line of reasoning shows that the associated map is given by A = 
ToA CP , with A CP a CP map with rank two Kraus operators Vi. We call maps A defined 
in this way two-decomposable^. Note that the adjoint map of a two-decomposable map 
is two-decomposable. It is now an easy exercise to derive the necessary and sufficient 
condition of one-undistillability in terms of two-decomposable maps. For completeness, 
we present the derivation. 

Theorem 4.2.1. A state p is 1-undistillable if and only if 

(i® A)(p) > 



1 This fact was independently obtained in Ref. 
2 This term has different meaning in other contexts, such as in Ref. 



179 [ 
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for all two-decomposable maps A. 

Proof. Suppose p is 1-distillable, so that there exists a D = (i/j\ Tb , with a 
Schmidt rank-2 vector such that Tr(Dp) < 0. Using the associated map this can be 
written as 

Tr[(l ® A)(P+)p] = TV[(1 ® A r )(p)P+] < 0, 

and since P+ is positive it follows that 1 (g) A^(p) ^ 0, i.e. we have found a two- 
decomposable map that detects the state. To prove the converse, let A = A CP o T be a 
two-decomposable map such that 1 <g> A(p) has a negative eigenvalue. Denoting by \4>) 
the corresponding eigenvector, we get 

® A CP (p TB M = Tr[(l ® At)(|0)(0|)p T s] < o. (45) 

It is sufficient to consider A CP (A) = VAV\ with V a rank-2 operator, therefore 1 <S) 
At(|0)(0|) = \(f/)((f/\, with 10') Schmidt rank 2. Thus Tr(p Ts \^'}((/)'\) < 0, and p is 
1-distillable. □ 

When D = |-i/')(^| Ts the associated map is of the form 

p^ {t®V)p TB {t®V) ] . 

It is instructive to compare this to Theorem 14.1.21 The above mapping is basically 
a projection of the partial transposed state on a C n ® C 2 (this special case was also 
derived in Ref. [3]). Theorem 14.1.21 characterizes distillability in an analogous manner, 
but instead projects onto C 2 (g> C 2 . Thus regarded in this way, the above theorem is 
both a generalisation and a simplification of the main theorem of distillability. 

Note that given a two-decomposable map A, the above result implies that undistil- 
lable states must satisfy both (1 ® A)(p) > and the dual criterion (A (g) !)(/>) > 0. 
As the reader might have guessed, the reduction criterion is just an example of a 
two-decomposable map, which we will show below. We also give numerous other ap- 
plications of our Theorem 14.2.11 and Observation 14.2. H starting with an elementary 
derivation of some alternative formulations of distillability. 
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4.2.2 Application 1: Derivation of other distillability theorems 

In this section we rederive two distillability theorems. The original derivations were 
quite lengthy, and as we will show, both theorems follow straightforwardly from our 
Observation I4.2.T1 and some basic theorems on Schmidt numbers from Chapter El The 
follo wing theorem was first derived in Ref. [3] , see also the related characterisation in 



Ref. 



2571 ] 



Theorem 4.2.2. Let p be a state and S be the completely positive map defined by 

p=(l®S)P+. 

Let T denote the transpose map. Then p is 1-undistillable if and only if A = T o S is 
two-positive, and more generally n-undistillable if and only if A® n is two-positive. 

Proof. For undistillable states p we have that 

Tr (pD) = Tr ( P Tb D Tb ) > 0, 

for all D Tb having Schmidt number two. From the isomorphism between states positive 
on Schmidt number k states and A;-positive maps (Theorem l2.2.2p . we deduce that the 
map associated with p Ts is two-positive for undistillable states. This map is just given 
by A as defined in the statement of the theorem. Thus p is 1-undistillable if and only 
if A is two-positive. The second part follows trivially by applying the first part on p® n 
and some elementary algebra. □ 



1771 ]. 



The second theorem is a slightly simplified version of the one presented in Ref. 

Theorem 4.2.3. Let P<i be the projector onto a maximally entangled state acting on 
T~ti = HAi ® Hfij = C 2 (g> C 2 . Then for an arbitrary operator X acting on Ti.2 = 
Ha 2 ®H.b 2 w e can define 

W x = P2®X Tb . (4.6) 
A state p acting on Ti.2 is n-undistillable if and only ifW p ®n is an entanglement witness. 

Proof. From Observation 14.2.11 we have that p is 1-undistillable if and only if p Ta 
is positive on Schmidt number two states. The theorem then follows directly from 
Lemma I2.1.2| and from the fact that p is n-undistillable if and only if p® n is 1- 
undistillable. □ 
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4.2.3 Application 2: A plethora of reduction-like maps 

In this section we show that the map associated with the reduction criterion is two- 
decomposable, which immediately proves that all states violating the reduction criterion 
can be distilled. This greatly simplifies the original proof which relied on a series of 
protocols: filtering, twirling and distillation of isotropic states. Let us construct the 
operator a from the Schmidt rank-2 vectors = — \ji), i ^ j- Adding the 

corresponding projectors we get 

i<j i^j 

From this follows that D = o Tb = 1 — dP + is positive on undistillable states. The 
associated map is given by Ai(j4) = Tr(A)l — A. This is just the map used by the 
reduction criterion for entanglement. Thus if a state p satisfies 1 (g) — p^Oor 
PA <8> 1 — p then it is 1-distillable. In this way we have also proven that the map is 
decomposable 1311 ] . and using our results one obtains readily the explicit Choi-Kraus 
form of the map. Next, we construct some other elementary examples by combining 
different Schmidt rank-2 vectors: 

<^> Taking our vectors = \ij) + we obtain 

° = E (1*^X^1 + ly>0»l) = i + f-2z, 

with Z = Y2i Therefore D = o Tb = 1 + dP + — 2Z and the associated map is 

given by A2(^4) = Trvll + A — 2diag A The map diag maps all off-diagonal elements 
to zero and leaves the diagonal itself invariant. 

(} For Schmidt rank-2 vectors of the form = \ii) + \jj), we get: 

° = E (i">< ii i + = dP + + ( d - 2 ) z > 

and thus D = o Tb = F + (d — 2)Z; the associated map is given by A3 (.A) = A T + {d — 
2) diag A 

<3> If we take Schmidt rank-2 vectors of the form = \H) — \jj), we get: 
a = Y, - = -dP + + dZ, 
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and thus D = o Tb = —F + dZ; the associated map is given by A^^) = — A T + ddmg A. 

<0> In fact, every operator with Schmidt number 2 gives us a strong distillability 
witness. In Chapter [2] we have seen that the special isotropic state 

a = (d- 2)1 + (2d - l)dP+ 

allows a Schmidt rank-2 decomposition. We find D = cj Tb = (d — 2)1 + (2d — 1)F and 
for the corresponding map A 5 (A) = (d-2) Tr At + (2d - l)A T . 

It is worth noting that although these maps are just like the reduction map two- 
decomposable, they lack any physical interpretation. The reduction criterion on the 
other hand, has several intuitive prop erties, linking with other notions in quantum 



information (see for instance 43, |231| |). As an example, note that the hash inequal- 



ity only gives a non-trivial lower bou nd on the entanglement of distillation for states 
violating the reduction criterion 



1371 ] . We believe therefore that the presented maps 
deserve further investigation. 

There are however a few general statements we can make about the introduced 
maps. As we have seen for the separability problem, given an entanglement witness 
D, the dual positive map A will detect much more states (see the discussion after 



Theorem II .2.5f) . In a similar fashion as is shown there, the map A corresponds to the 
class of witnesses (A (g> t)D(A^ <g> 1), for arbitrary A. It also implies that the criteria 
1 <S> A(p) > and A <g> l(/t?) > are insensitive to local transformations by one of 
the parties (thus to filtering operations). Indeed, suppose Alice performs a general 
measurement, with measurement operators A\ and A2 satisfying j4{^4i + A 2 A2 = 1. 
Then the state p will be transformed into pi = (Ai (g) l)p(A| <S> l)/pi. It follows that, 
if the original state did not violate the criteria, then the transformed state does not 
either. The map corresponds to the operator witness, together with all possible local 
filtering operations by one of the parties. 

We have shown how to construct maps that detect distillability, applicable in arbi- 
trary dimensions, which can be easily evaluated on states. There is however a catch. 
As can be seen from the reconstruction of the reduction criterion, the witness from 
which the map is derived, is a convex combination of many witnesses. So it is possible 
that one of those witnesses (and the associated map) detects a state, while the sum 
does not. In other words, the map could be weaker than expected at first sight. The 
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sum of the negative eigenvalues of D, or the entanglement of D Tb , can be seen as a 
measure for the strength of a certain map. 

The basic requirement for a distillability witness D to detect a two-distillable state 
p which is not 1-distillable is that D Tb must be entangled with respect to the first and 
second pair. Indeed suppose D = \ip)(tp\ TB with \tp) = \tpi) ® 1^2) separable, then 

Tt(D Tb P ® 2Tb ) = Tr(]ih)(ih\p)Tr(\ih)(ih\p) > 0, 

since the vectors \ipi) have at most Schmidt rank 2. Consider the witness for the 
reduction criterion on two pairs: 

D Tb = 1 - F = li ® 1 2 - F 1 ® F 2 . 

But from Ref. 



1381 ] we can write F = P s - P A and 1 = P s + P A with P s and P A 
projection operators onto the symmetric and antisymmetric subspaces respectively. 
Substitution yields 

D Tb = Psx ® P A 2 + P A \ ® Ps2, 

so that D Tb is separable with respect to the different pairs. Note that this property is 
not adhered by the map since {A® 1)D Tb (A^ ® 1) could be entangled with respect to 
the two pairs, even if D itself is not. So it is not at all obvious that collective application 
of the reduction criterion (for instance on p® 2 ) does not yield a stronger criterion. 



It was proven in Ref. [13 it ] that it is not the case, namely, if 1 (8> Ai(/9® 2 ) ^ then 
1® Ai(p) ^ 0. It is readily verified that all introduced maps share this property. In the 
case of the reduction criterion, the reverse is also true. This can be proven as follows. 

n 

From Ref. [1311 ] we have 

1 A <8> A B (pi ® Pi) = [Iai ® pm] ® [tA3 ® A B2 ](p 2 ) + [Iai ® A B i](pi) 8> P2- 

Now suppose 1® A(/9j) is not a positive operator, so that there exists a vector \ip) with 
negative expectation value. From the above expression follows that ® A B (pi ® 

/O2) | x I r ) < 0, with |^) = | -01 } ® l^)- So applying the reduction criterion to one or more 
pairs is completely equivalent. A similar statement and proof applies for A4, but for 
the other maps, this reverse statement is not true. 
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4.2.4 Application 3: Random Schmidt-2 robustness and distillability 

The following proposition, first presented in Ref. [61] uses upper bounds to the random 
Schmidt-2 robustness as defined in Section [2,4.21 to construct a distillability criterion: 

Proposition 4.2.1. Let p be an arbitrary bipartite state, such that p Ta has negative 
eigenvalues. Let be the eigenvector corresponding to a negative eigenvalue A and let 
R r 2{ip) be an upper bound to its random Schmidt-2 robustness. Then p is 1-distillable 
if 

A < -Rr 2 (ip). 
Proof. From the definition, we have that 



W 



has Schmidt number two. Now from Observation 14.2.11 follows that p is 1-distillable if 
Ti(Wp TB ) < 0. This can we rewritten as 



+ Rr 2 (lp) = \ + Rr2(lp) < 0. 



□ 



This proposition provides an important reason to find an exact analytical formula for 
Rr2 (VO • Note that this distillability criterion only depends on the minimum eigenvalue 
and the corresponding eigenvector. It is easy to see that it is independent of the 
reduction criterion as we show by means of an example. Consider the state 

1 00 000000 

20 -1 00000 

02 000200 

0-10 200000 

00 010000- (4.7) 

00 002020 

02 000200 

00 002020 

00 000002 



'=16 
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It is easy to check that the reduction criterion is not useful here. The partial transpo- 
sition p TB has an eigenvector = -^=(|00) + |11) — |22)) and corresponding eigenvalue 
— g. Since —R r 2{ip) = —j^ > p is distillable. 

An unconditional distillability criterion, which does only depend on the minimum 
eigenvalue of p Te would follow from Conjecture 12,4.11 Pending on the proof of this 
conjecture, we have that p is distillable whenever p Te has an eigenvalue A < —1/14. 

4.2.5 Application 4: Volume of 1-distillable states 

In this section we will outline a powerful algorithm for the detection of distillable states 



(first published in Ref. 6CJ]). As before, it is clear that we can confine ourselves to the 
study of 1-distillability. A first numerical method for this problem was devised by Diir 
et al. 



starting from the basic minimisation over Schmidt rank two vectors 



min {tp\p TB \ij) < 0. (4.8) 
ip£SR2 

They convert this problem into a minimisation of the minimum eigenvalue of certain 
matrices. Their method involves also calculating an inverse of a matrix in every step. 

In view of our discussion above, we will perform this minimisation not over the 
distillability witnesses, but over the much stronger associated two- decomposable maps. 



Thus for a given state p, we can reduce the problem (see also again Ref. 



721]) to checking 



whether there exists an operator P = |0)(a| + with (b\a) = such that 

? 

(1® P) P Tb (1® P) j ■£ 0. 
One way of doing this is to parametrize a countable s ubse t of vectors which is dense 



within all vectors, such as the one introduced in Ref. 
the form G = Aj|£)|(Ai, A 2 , • • • , Ad) G Gn} with 



1591 ] . Explicitly their set takes 



, -, i . Pi 2iri r -± P2 2wi^ 

G N = { \ —e s i , , — e s 2 . 

qi Q2 



for < pi < qi < N; < rj < s, < N. Thus for every N, we can construct sets of 
pairs | a) and \b), and taking N increasingly large we will detect all 1-distillable states, 
except those arbitrary close to the boundary of the convex set of 1-undistillable states. 
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Figure 4.1: The probability of detecting 1-distillability as a function of the test number 
for random states drawn from T> x U. For visual purposes only the first 100 tests are 
depicted and the points have been smoothened out to curves. 



In practice two improvements can be made which greatly enhance the performance. 
First note that the countable subset above will yield vectors \a) and \b) not necessarily 
orthogonal. Furthermore, clearly such a countable subset will in general not pick vectors 
uniformly distributed according to the Haar measure. One way of overcoming this is to 
take for \a) and \ b) two columns of a random unitary. Here orthogonality and uniformity 
are automatically guaranteed. The associated algorithm works well for low dimensional 
density matrices. For higher dimensions, a local optimisation of the \a) and \b) yielding 
the minimum after a certain cut-off number of tests turns out to work well. The explicit 
algorithm has been reproduced in Appendix IA.2I 

Let us apply this to give a numerical estimate of the volume of 1-distillable states 
for low dimensional qua ntum sta tes. A similar numerical estimate has been carried 
out for entangled states 



339 



342]. When talking about volumes on the set of density 
operators, it is clear th at t he results will depend on the measure. We choose the 
measure applied in Ref. [3421 ] as it seems very natural. 

A general bipartite quantum state p acting on TCa ^Ti-B can be expanded by virtue 
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Figure 4.2: Probability of a random NPT state drawn from T> x ti to be 1-undistillable. 
The curve has been drawn to guide the eye. 



of the spectral decomposition as 



UDU\ 



with D = (da) diagonal and U unitary. The measure we are going to use is the 
product measure T> x U. Here T> represents the uniform distribution of the points on 
the manifold given by Yli da = 1- A simple method for generating such a distribution 
from independent uniformly distributed random numbers chosen in the interval (0, 1) is 



outlined in Appendix A of Ref. 



3421 ] . Similarly, U is chosen to be the uniform measure 



on unitary matrices (the Haar measure). To g enerate ra ndom unitaries according to this 



2371 . 



3411 ] which relies on a decomposition 



measure one can use the algorithm from Ref. 
of a general unitary in two-dimensional unitary transf orma tions . A simpler method 
for generating Haar unitary matrices is as follows 



86, 



2751 ]. Take a random matrix 



A, whose entries are complex numbers which are independently normally distributed 
with mean zero. The polar decomposition A = PU, with P positive, then yields U 
distributed according to U. 
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Figure 4.3: Probability of a random state drawn from T> x U to be 1-undistillable. 



We have tested the our method on 10 5 density matrices acting on C d (g> C d for 
d = 3,4, 5, 6, 7. The estimated probability for proving the distillability of a state as a 
function of the test number is displayed in Figure l4~Tl A striking feature is that the vast 
majority of the distillable states were detected in the first few tests. This is especially 
true for low dimensional states. For example for d = 3 about half of the NPT states are 
found to be distillable in the first test. For d = 7 the probability of finding a distillable 
state in the first test is about 1/6. 

In Figure I3~2l the probability of an NPT state being 1-undistillable is plotted versus 
the dimension d. To obtain sufficient precision we carried out 10 5 random tests per 
state, and in addition 10 4 <i optimisation steps seeking for a local minimum. Of course, 
this method does not guarantee to detect every 1-distillable state, but we obtain an 
upper bound of the number of undistillable states. The figure suggests that the volume 
of 1-undistillable states drops to zero for high dimensions. In Figure FOl the same graph 
is drawn, but now PPT states are included. 

It is tempting to conclude from this numerical evidence that bound entangled states 



4.3 The question of NPT bound entanglement 



126 



are primarily a phenomenon present in low dimensional quantum systems. In high 
dimensional systems most undistillable states are therefore situated in the immediate 
neighbourhood of the set of separable states. Our results are consistent wit h the fact 
that bound entanglement for continuous variables is a rare phenomenon [153]. in 
particular it was shown that the subset of undistillable states is nowhere dense in 
the set of bipartite continuous variable states. From this it follows that the set of 
undistillable continuous variables states does not contain any open ball, an argument 
which was made explicit for separable states in Ref . 92J . Given an infinite dimensional 



separable state one can construct sequences of closer and closer states all of which are 
entangled. Following the same methods one can explicitly construct a distillable state 
in a ny e- neighbourhood in the trace norm of any state [89]. However, note that in 



Ref. 



1211 ] a parametrized family of measures on states was introduced, which in some 



region yields states primarily 1-undistillable. 



4.3 The question of NPT bound entanglement 
4.3.1 Early results 



As we have seen in Theorem 14. 1.4} the question of the existence of NPT bound entangle- 
ment is equivalent to the existence of undistillable entangled Werner states. Numerical 
evidence seems to suggest that some Werner states are indeed undistillable. However, 
the problem is still open, and constitutes one of the major challenges in quantum infor- 



mation theory. What has been shown 
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is that there are entangled Werner states 



that are n-undistillable for every n, however the range of states for which one is able 
to prove this, vanishes as n goes to infinity. We will prove this partial result following 



the method presented in Ref. 



831 ] . We start with the following key lemma. 



Lemma 4.3.1. Let P be the projector onto the maximally entangled d-level state and 
Q = 1 — P, then for any Schmidt rank two vector \ip) we have: 



1. (ip\l® N - k ® P® k \il;) < 



2. {lP\Q® N ~ k ® P® k \lP) < Jr 



3. WP N Wi > (i - z) N 
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Proof. (1) For k = 1, this is clearly the case, as the overlap of the maximally entangled 
state and a Schmidt rank two vector is at most 2/d (see Lemma l2.1.1j) . For k > 1 this 
follows from the case k = 1 and Pf k = P d k . 

(2) This follows from Q < 1. 

(3) For N = 1, this follows directly from (1) and the definition Q = 1 — P. We 
prove it for arbitrary N by induction, so let it be true for a certain N > 1. Now 

Writing this trace as Tr(-) = Tr(Tri(-)), where we first trace out over the first pair, we 
get 

TMMMQ ® Q®"] = Tr[Tr 1 (|^)(V'|Q)Q 0jV ]. 

Now since Q is separable (see (|2.ip ). it is easy to verify that Tri(\ip)(ip\Q) can be 
written as Y2i CilV'i)^! with ^ Schmidt rank two. Also from the case iV = 1 follows 
that Tr\EiCi\il>i)(il>i\] >1 -2/d. □ 

This lemma is practically all we need to present the main result of this section. 

Theorem 4.3.1. The Werner states 



with —1 < (5 < 1 are entangled for (3 < —1/d or —1 < Tr(/aF) < 0. T/jey are 1- 
distillable if and only if [3 < —1/2 and n-undistillable for > —1/d + e n , for some 
e„ > 0. 

Proof. The first part follows from the fact that positivity of the partial tra nspo sition 



is a necessary and sufficient condition for separability of the Werner states [309(. The 
Werner states are 1-distillable for (3 < —1/2, as follows from the overlap of the partial 
transposition with a Schmidt rank two vector \tp) 

<V|p TB W = 1 + PWdP+ty) > 1 + 2/?, 

where we have used Lemma 12.1.11 For n copies let us put a = 1 + @d and 1 = P + Q 
so that 

p T B ®n = g®n + aR 
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Here R consists of all the other terms in the tensor product expansion. Note that all 
these terms contain at least one P. From the previous lemma follows then that for any 
Schmidt rank two vector 

MfP^W)>(l-l) N + aW\R\1>)- 
Applying again the previous lemma, we also have that we can calculate a lower bound 
for a(tp\R\ip). Thus we can always choose a sufficiently small so that {il)\p TB ® n \ip) > 
for which p will be n-undistillable. We refer to the original reference and our next 
section for qualitative results on the value of e n . □ 

From this theorem it was conjectured that all Werner states with (3 > —1/2 are 
undistillable, although e n — > for n — > oo. Apart from some numerical evidence, the 
main reason was that at the time these results appeared, not a single example of a 
state was known which is 2-dis tilla b le b ut 1-undistillable. In the meantime, several 
such states have appeared 
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3101 ] . These first examples arose in th e co ntext 



of activation (see Section l4,4p . Important progress was made by Watrous 3161 ] who 
constructed a one-parameter set of distillable states which are n-undistillable in some 
range. We will discuss these states in the next section. Watrous's result does not seem 
too unlikely, as the notion of n-distillability, as we explained before, is really as much a 
mathematical as a physical notion. A distillation process will typically require a very 
large number of copies, even if the state is 1-distillable. 

4.3.2 The UUVVF-invariant states 

In this section we will discuss the distillability properties of a certain class of highly 
symmetric two-parameter states, which includes the Watrous states and two copies of 
the Werner states, and hence is ideal for the study of distillability. We will see that over 
a wide range of the parameters, the states are probably undistillable. Other examples 
of conjectured NPT bound entangled states can be found in Ref. |l4 |. 

The set of states we are going to study are the C/C/V^-F-invariant states introduced 
in Section 14.11 A convenient parametrisation of the set of UUVV F-mv&riant states is 
given by 

ed - 1 . , 1 - 2ed + 5d 2 
p = l 12 ® 1 3 4 + —j— (I12 ® F 34 + F12 ® I34) + j 2 F 12 (8) F 34 . 
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The set of density operators is restricted by the following inequalities 



(d - l) 2 + 2ed{d - 1) + 5d 2 > 
d 2 - 1 + 2ed - 5d 2 > 
(d + l) 2 - 2ed(d +l)+5d 2 >0 



This set of states includes the Werner 



3211 ] st ates 



pair m d 2 ® d 2 for e = 1/d) and the Watrous |31Jl sta tes (1 - 2ed + 8d 2 = d 2 ). The 



(2 pairs in d ® d for 5 = e 2 and 1 



309]). The states are entangled for 



PPT states are just the separable states (see Ref. 
e < or 5 < which follows from 

P Tb = Q\2 ® Qu + de(P 12 ® Q 3 4 + Q12 ® P34) + 8d 2 P 12 P34, (4.9) 

with as usual Q = 1 — P. 

Let us now investigate the distillation properties of these states. To prove distilla- 
bility, all we need to do is to find Schmidt rank two vectors ip such that (ip\p TB < 0. 
We present three such vectors, which we conjecture to detect all 1-distillable states: 

W)A = |00)i 2 ® (|00) + |11)) 3 4 J 

= |01)i2<E> (|00> + |11)) 3 4, (4.10) 



Wic = \ u )A\jj)B + + 1))aU0' + 1))b- 

ij 

which detect distillability for the states respectively satisfying 

d 2 + 3d(ed - 1) + 2(1 - 2ed + 5d 2 ) < 0, 
1 1 

, 1 1 

6< ^~2- 

This set of states is shaded dark grey in Figure 14.41 

Now let us derive results in the other direction, namely which states are undis- 
tillable? First let us show that for arbitrary n we can find NPT states which are 
n-undistillable. We will use the inequalities from Lemma 14.3.11 From (|4.9p it follows 
that all potentially negative terms contain at least one P term. But we will also have 
a term (\I/|Q® 2 " |^) > (1 — f) 2Af , and this term can always dominate when we choose 
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e and 8 small enough. Thus for each n, as long as we choose e or 5 small enough (but 
one or both negative), we obtain an n-undistillable NPT state. 

Next, let us derive states that are definitely 1-undistillable. From (|4.9fl and the 
given inequalities it follows straightforwardly that the states satisfying 

(1 - 2/df + min(4e, 0) + min(2<5, 0) > 0, 

are 1-undistillable. Note that this set does not touch the set of 1-distillable states. The 
reason for this is that although the inequalities are sharp, the sum of such inequalities 
are not. In some regions we can get a better bound by rewriting (|4.9f) as 

P Tb = I12 <8> 1 3 4 + (ed - l)(Pl2 ® Q 3 4 + Q12 ® P 34 ) + {Sd 2 - l)Pi 2 ® P34, 

and we find that all states such that 

d 2 + min(4d(ed - 1), 0) + min(2((5(2 2 - 1), 0) > 



are 1-undistillable. This second set of states is depicted in Figure | 

To conc lude our discussion of the distillation properties we will discuss the Watrous 
states 316] of which the general form is given by 



p = ll, 2 <S> 1 3) 4 + ^-^(ll,2 ® ^3,4 + Pl,2 ® 13,4) + Pl,2 8> P 3 ,4, 

with 1 + 1/d > e > l/d—l. The states are entangled if and only if e < and definitely 
1-distillable if e < l/d — 1/2. We will show that all entangled Watrous states are 
distillable. Suppose now we have two pairs of entangled Watrous states, the second 
pair having indices 5, 6 and 7, 8. Using the identities 

Tr((P 1)5 ® P 2 ,6)(li,2 ® F 5>6 )) = l/d, 
Tr((P lj5 ® P 2 , 6 )(Pi,2 8> 1 5 , 6 )) = 1M 
Tr((P lj5 P 2 , 6 )(i ? i,2 ® P 5 ,6)) = 1, 

it can be verified that projecting upon ^ ® P2 5 6 will yield a new Watrous state p' 
with parameter 

, ( 2(ed + d 2 -l) \ 
6 ~ e \d 2 e 2 + d 2 -l )' 

Since e' < e for > e > l/d — 1, the state will be more entangled than the state 
we started from for each e < 0. One can repeat this protocol on many pairs until 
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Figure 4.4: UUVV ^-invariant states. All states satisfying e > and 5 > are separa- 
ble. Shaded areas mark out distillable states (1,2 and do- copy distillable). The states 
lying above and to the right of the dotted line are states for which we were able to show 
1-undistillability. 



4.3 The question of NPT bound entanglement 



132 



finally e < 1/d — 1/2, at which point we obtain a 1-distillable state. More generally 
the protocol can be applied to the whole set of states; a straightforward but tedious 
calculation leads to 

, _ e{d 2 5 + d 2 - 1) 
6 ~ d 2 e 2 + d 2 -l ' 
e 2 {d 2 -l)+d 2 5 2 
d 2 e 2 + d 2 -l ' 

The states which are 2-distillable with this protocol are depicted in Figure 14.41 Re- 
peating the protocol recursively, it is not hard to show that all states satisfying 

3d 2 + Ad - 8 1_ 
> 2d(d-2) £+ d?' 

are distillable (see again Figure 14. 4|) . In Section 14.2.31 we have seen that the basic 
requirement for a distillability witness W, to detect 2-distillability of a state p (or thus 
1-distillability of p® 2 ) which is not 1-distillable, is that W Tb be entangled with respect 
to the two pairs. The Watrous states allow us to illustrate this. Indeed from the above 
protocol we can deduce that Wi on p® 2 will detect 2-distillable states that are not 
1-distillable, with 

W 2 = (P 15 ® P 2 T 6 ® P B ,3478) T2468 . 

Here T is the normal transposition, T2468 is the transposition with respect to the sub- 
systems 2, 4, 6 and 8 and Pb is the projector onto b defined in (|4.1U|) . The projectors 
-Pi 5 <S> -P26 are responsible for the entanglement between the first and second pair. 

Using the numerically method outlined in Section ^. 2. 5l we have checked numerically 
the distillability of the C/c^^y-F-invariant states over the complete range of parameters 
for 1 and 2 copies for d = 3. We easily recovered the proposed boundaries for 1- 
distillability. For two copies, the states act on a C a ® C a Hilbert space and numerical 
matrix manipulations in a space of this magnitude seem very hard. Fortunately, the 
states are sparse and our method only requires minimisation of the minimum eigenvalue 
of a 2d 2 x 2d 2 matrix. We were able to detect distillability for the Watrous states in 
e > 1/d — 1/2 readily and exhaustive testing suggest that also the proposed boundary 
for 2-distillability is correct. This for the first time provides strong evidence that the 
Werner states are 4-undistillable for j3 > —1/2. 
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4.3.3 As a Schmidt number problem 

As we have seen, the NPT bound entangled conjecture is equivalent to the statement 
that there exists no Schmidt rank two vector \ip) such that 

(V|/9 Ts0n |V> < 0, (4.11) 

for all n > 1 and p Ts = 1 — |P. In other words, it seems that affirmation of the 
conjecture would have to be in the form of an impossibility proof as opposed to a 
constructive proof. We will reformulate this conjecture in a more tractable form, namely 
as a special instance of the separability problem, for which a large number of tools are 
present. As a steppingstone we first show how to translate the distillability problem into 
the problem of detecting Schmidt number 3. In the next subsection we then reformulate 
it as a separability problem. We will do so for one, two and n copies of the Werner 
states. 

Let us start with one copy of p, for which the answer is known. As p belongs to 
the set of C/?7-invariant states, p Ts will belong to the set of ?7i7*-invariant states (the 
isotropic states) and hence will be invariant under the C7C/*-twirl: Tjju*(p Tb ) = p Ts . 
From this follows that 

W\p TB \t/}) = Tr(P^ T *) = Tt(T tw .(i^)/ B ), 

here Tjju*{P^) is the operator P^ = after application of the C/C/*-twirl. Thus 

we do not need to check over the whole set of Schmidt rank two vectors, but instead 
over the restri cted set of isotropic states with Schmidt number 2. Recall (see (12. ip and 



Ref. 



2571 . 



2861 ] ) that if we parametrise the isotropic states as p a = 11 + otP then p a has 



Schmidt number k when 

d{kd- 1) 



a < 



d-k 

From this follows that p = 1 + (3F is 1-distillable if and only if Tr(p TB p a ) < 0, with 
a = d ^d-2^ ■ Going through the algebra we recover that all Werner states pp with 
13 < -1/2 are 1-distillable. 

Let us now look at two copies of the Werner states p® 2 ] as pointed out before, these 
states belong to the larger class of the C/C/yi^P-invariant states. Thus the relevant 
dual set is the set of the UU*VV* F- invariant states. For convenience we write the 
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operators in the order of the indices 1,2,3,4 and omit these indices. With this in mind 
we parametrise the UU*VV* F-inv&ri&nt states as 

P = Q®Q + x(Q ®P + P<g>Q) + yP<g>P, 

with x,y > 0. These states are separable for d 2 — 2d{x — 1) + (1 — 2x + y) < and 
d 2 — (1 — 2x + y) < as depicted in Figure 1431 The set of Schmidt number 2 states 
contains at least the convex hull of the points (x,y) 

= / (3d-4)(d+l) 2(d + l) 2 (d-l) \ 
V 2d -4 ' d-2 ) 

obtained by twirling the Schmidt rank two vectors from equation (j4. 10p . All the states 
lying above the CA-line have a Schmidt number larger than two. This follows from 
the fact that the operator 1 ® 1 — 4j-P ® P is positive on Schmidt number 2 states. 
If the Werner states are 2-undistillable then also all the states lying to the right of 
the AB -line have Schmidt number at least three. This follows easily by evaluating the 
expectation value of p at (1 — |P)® 2 . Note that it is sufficient to show that all the 
states p on for instance the I?I}-segment (see Figure l4~5|) have Schmidt number 3. This 
can be seen as follows: let p approach E, and suppose we can show that each p has 
Schmidt number 3. This would imply the existence of a hyperplane W separating p 
from the set of Schmidt number 2 states. For p arbitrarily close to E this hyperplane 
would be parallel to the AB-segment, otherwise cutting it. Therefore this would show 
that all the states lying to the right of the AB-hyperplane have Schmidt number at 
least three. 

For the general case we need to consider the set of U\U* ■ ■ ■ U n U*F-'mvariant states. 
Here the subindex in Ui refers to the subsystem the unitary operator is acting and the 
F denotes any permutation of the subsystems. In what follows we will call these states 
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C/jC/^F-invariant states. We parametrise them as 

p = Q®« + Xl (P &> Q^- 1 + Q&P® Q® n ~ 2 + •••)+ 
x 2 {P m ® Q® n ~ 2 +P®Q®P® Q® n ~ 2 + •••)+ 

• • • + x n P® n , 

where we have found it convenient to now use the normalised Q = Q/(d 2 — 1). The 
relevant hyper plane is given by 

i=i v y 

We can also write this as 1 + X^ILi (^)( — = with ctij = (f[ — lYxi. From now 
on we will continue to work in these normalised variables. Next we will generalise the 
idea developed for two copies. First we need to show that the hyperplane is spanned 
by Schmidt number 2 states. Then in order to check distillability, it will be enough to 
find the boundary between two and three Schmidt number along a line from the origin 
to an interior point of the points spanning the hyperplane. 
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An independent set of Schmidt number 2 states spanning the hyperplane is easily 
obtained as follows (compare to the case for two copies). Let = -^lOl)®™ -1 <E> 
(|00) + 1 11)). Twirling this state will yield an L^L^i 7 - invariant state with coordinates 
(^,0, • • • ,0). Similarly, twirling = -^|00) <g> |01)®"- 2 <g> (|00) + |11)) will yield an 
C/jC/^F-invariant state with coordinates (x\ ^ 0, x 2 ^ 0, 0, • • • , 0). In general 

v2 

will yield a state with coordinates x% ^ for i < + 1 = for i > + 1. It is 
evident that all points will lie on the hyperplane and that they form an independent 
set, spanning the hyperplane. An interior point can for instance be obtained from the 
first point, as + e(l — -), e, • • • , e) for sufficiently small e. It can readily be verified 
that this point belongs to the hyperplane by using the identity Y17=2 (?)( — -0* = 

Thus n-undistillability of the Werner states beyond the 1-distillability boundary is 
equivalent to the statement that the [/jCT^F-invariant state with coordinates + e(l — 
i) + S, e, • • ■ , e) has Schmidt number 3 for e > small enough and all 5 > 0. 

4.3.4 As a separability problem 

In this section we cast the distillability problem as a special instance of the separability 
problem. Theorem 14.2.31 provides the key ingredient to this, which applied to the 
conjecture implies that 

W x = P2®X Tb , (4.12) 

is an entanglement witness for all n > 1 with X = p® n and p = (1 — F/2). We will 
again use the local symmetry, and present a dual positive formulation of the conjecture. 
The analysis will be analogous to the reformulation as a Schmidt number problem and 
a continuous comparison of this section with the previous one is very instructive. 

For one pair, we need to prove that P2 (1 — |-P) is an entanglement witness. 
As before, it will be sufficient to characterize the subset of the separable states of the 
general UU*VV* invariant states. Here U acts on a two-dimensional Hilbert space. 
We parametrize the UU*VV* -mvariant states as 

P = Q2 ® Q + xP 2 <g> Q + yQ 2 <8> P + zP 2 P, 
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Figure 4.6: The state space of the [/C/*yV r *-invariant states. The point C lies in the 
xz-plane such that the points A, B, C and E lie in the same plane. 



with x,y,z > 0. The separable states are a subset of the states with positive partial 
transpose, which satisfy the inequalities 

1 + 3x 

3z < — y ABCE-plane 

d — 1 

I — x 

z < h y ECD-plane 

d+1 

x — 1 

z > -; h y EAD-p\ane. 

d — 1 

Thus the states with positive partial transposition are contained in the polyhedron 
spanned by the points OABCDE as in Fig 14.61 By twirling the pure separable states 
|0101), |0100), 1 0001) and X^ij=olW7') one °bt ams the states represented respectively 
by the points O, A, C, D and E. Note that B is not in this list. Indeed, as we know 
that P2 ® (1 — f -P) is an entanglement witness, all states satisfying z > 2x/(d — 2) are 
entangled. In particular the states in the tetrahedron spanned by the points ABCO are 
PPT entangled. Conversely, knowing that the polyhedron ABCO is PPT entangled 
immediately proves that P 2 ® (1 - f P) is an entanglement witness. 
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For two copies, we need to consider the UU* (ViV*V 2 V£ F)-'mvaii&nt states which 
we parametrise as 

P = Q 2 ® {Q ® Q + yi(P ® + Q ® -P) + 2/2P ® P)+ 

P 2 ® (x Q ®Q + a:i(P0Q + QoP) + x 2 P <g> P). 

We will not attempt to completely classify the separable subset, instead it is enough to 
look at what happens in the neighbourhood of the hyperplane 



Tr 



/ / d\ \ ® 2 l /d — 2\' z 

P 2 ® [Q + (l - -JPJ pj = x - (d - 2)xi + (-g-j ^2 = 0. 



d-2\2 



We now show that this hyperplane is spanned by separable states. Consider the fol- 
lowing pure separable states and their coordinates (yi,y 2 , xo, x±, x 2 ) after action of the 
UU*(V 1 V 1 *V 2 V 2 *F) twirl: 

|01)|01)|01) : (0,0,0,0,0) 
|01)|00)|01) : ^-^L_, 0,0, 0,0^) 

|oi)|oo)|oo): (^^p ^rr, 0,0,0) 

^ , . , . ( 1 d - 2 3d - 4 2 \ 

i,j=0 \ / 

These coordinates can be verified after a tedious but straightforward calculation. 

A point in the interior of the convex hull of these points can be obtained by averaging 
these coordinates. In this case one obtains the state with coordinates xq = 2(d — 
2)/(15d), xx = (5d - 6)/(30d(d - 1)), x 2 = 2/(15d(d - 1)), yi = d(2d - 3)/(5(d - 
l)(d 2 — d — 1)) and y 2 = 1/ (5(d 2 — d— 1)). From this follows that the Werner states are 
2-undistillable if and only if the UU* (ViV*V 2 V2 ^-invariant states with coordinates 
(yi,y 2 ,xo,xi + e,x 2 ) are entangled for all e > 0. 
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Let us now move to n copies. The relevant set of states is the set of UU* (ViV* F)- 
invariant states 

P = Q2 ® Q® n + Q2 ® [yi(P ® Q® n ~ l + Q®P® Q®"- 2 + •••)+ 
y 2 (P m ® Q® n ~ 2 + •••) + ••• + VnP®^ 1 } +P2® [x Q® n + 
Xl (P ® Q® n ~ l + Q®P® Q® n ~ 2 + •••) + ••• + x n P® n ], 

with Xi,yi > 0. The relevant hyperplane is given by 

i=i v 7 

Renormalising Xj, this can be rewritten as Y17=i (1) (~ ^Y x i = with Xj = — l)*Xj. 
Next we will show that this hyperplane touches the set of separable states by construct- 
ing a set of 2n + 1 separable states spanning the hyperplane. The first n + 1 states are 
obtained by twirling 

\ij)) k = |oi)|oo)® fc |oi) 0n - fc , 

for k = 0, ■■■ , re. The twirled state will be ?7?7*(V^V^*i ? )-invariant and will satisfy 
Xj = for all i and / for j < /c and = for j > fc. The last n states are 
obtained by twirling 

1 

i,j=0 

for = 0, ■■■ ,n — 1. The twirled state will be LW*(V^F)-invariant and the x 
coordinates will satisfy Xj = for j > n — k and Xj / for j < n — k. Therefore 
the coordinates of the 2n + 1 states are linearly independent and the choice of \4>)k 
guarantees that the states will lie in the hyperplane. 

An interior point in the convex hull of these 2n + 1 points can be obtained by 
choosing xo = e and jjj = xj = e, for e sufficiently small. One verifies that this point 
belongs to the hyperplane using the identity Y17=o ©( — -0* = ~~ 1- From this follows 
that the conjecture is equivalent to the statement that the UU* (ViV* F)-mvaxiaxit states 
with coordinates y« = e , Xj = e and x\ = e + 5, where j = 1, ■ • ■ , n and i = 0, 2, ■ ■ • , re, 
are entangled for e > small enough and all S > 0. 
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Discussion 



1 10f | . However, 



The general separability problem has been proven to be NP-hard 
above we showed that the distillability problem can be reformulated as the question 
of entanglement of a particular set of one parameter states. The first and by far still 
the m ost elegant tool for detecting entanglement is the partial transposition criterion 



2221 ] . H owev er, it is not useful to solve the dual entanglement problem, as it was proven 



in Ref. 



1771 ] that W p m is a decomposable entanglement witness if and only if p has 



positive partial transposition. Hence W ®n will be either a non-decomposable witness 

Pw 

or no entanglement witness at all depending on whether p^ 1 is 1-undistillable or 1- 
distillable. Thus in order to solve the dual problem, we will either need some powerful 
tool for detecting PPT-entanglement or a tool for proving separability. However, in the 
latter case, the original formulation (see (|4. llj) ) seems easier, and for this purpose we 
presented an efficient algorithm for detecting distillability in Section [4.2.51 

Probably the most powerful method for detecting entan glem ent is the complete 
family of separability criteria introduced by Doherty et al. jQ, Q], as reviewed in 
Section 11.2.51 Recall that their method yields an analytical provable entanglement 
witness from the output of the algorithm. In principle therefore, the distillability 
problem can be solved for any number of copies, using our dual formulation together 
with the algorithm associated with the complete family of separability criteria. The 

can also be employed for this goal. 



numerical method introduced by Eisert et al. 
However, note that the reformulation of the distillability problem as a separability 
problem is unnecessary here, as their method allows for a simple way of checking 
whether or not an operator is an entanglement witness. Using analogous techniques, it 
can be extended to testing whether or not an operator is a Schmidt number 3 witness 



891 ] . It should be stressed that the family of criteria introduced by Doherty et al. has the 
major advantage that an analytical canonical entanglement witness can be extracted. 
Unfortunately, we were unable to test either algorithm, as both methods seem only 
practical for low dimensional systems. The interesting case, two copies of the Werner 
states for d = 3, could hence not be tested. It would be worth investigating whether 
these numerical methods could be simplified for the highly symmetrical states we are 
interested in. 
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Numerical solutions are one option, another possible approach would be one of a 
more indirect nature. A powerful method for proving that a certain state is entangled 
is to show that, when shared by two parties, the state can enhance typical quantum 
operations such as teleportation or d istill ation. We will come back to this topic in 
Section T4. 41 One of the main results is 



2(5 ] 



2031 ] that every entangled state can enhance the 
so-called conclusive teleportation fidelity of some other state. This characterisation of 
entangle d st ates is promising as a way of proving that a certain state is PPT entangled. 
In Ref. 1651 ] a class of PPT states was constructed which was shown to provide overall 
convertibility of pure entangled states. In particular it was shown to be able to increase 
the Schmidt number of a pure state. Now the PPT entangled states we obtained for 
one copy of the Werner states (tetrahedron ABCO in Figure B~6j) are of a similar form, 
and although not necessarily in an optimal way, they too can help increase the Schmidt 
number of a pure state. Similar activation effects can be expected from the conjectured 
PPT states derived from two copies of the Werner states. 

4.3.5 The rainbow states 

In this section, as a curiosity, we present a set of states which we conjecture to exhibit 
all the different forms of distillability. The set of states was inspired by the results of 
Ref. [1651]. 



The states are t7£/W^-invariant, but with subsystems of unequal dimension (we 

omit the indices {1,2} and {3,4}) 

de — 1 „ me — 1 _ 1 — (m + d)e + dm5 _ _ 

p = t m t d + — — l m ® F d + F m <g> t d + i —!■ F m <g) F d . 

d m dm 

In what follows we will assume that 3 < m < d. The set of density operators is 

restricted by the following inequalities 

1 , x m + d 
l + S + 2e + — - e + 1)^- > 
md md 

1 m+d r 1 1 
1- — + e— --5 + --->0 

md md m d 

1 . . m + d 

l + S-2e + — + (1- e) — > 

md md 

The partial transpose is given by 

P Tb = Qm ®Qd + meP m ®Q d + deQ m (g) P d + md5P m (g> P d , 
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and it is easy to see that the states are NPT iff e < or 5 < 0. Now we will show 
that these states also include some PPT entangled ones. It is easy to see from equation 
(I2.1.ip in Section \2.1\ that l d — dP d /m is positive on Schmidt number m states; it then 
follows from Lemma 12.1.21 that 

mP m <g> (t d - —P d ) 
m 

is an entanglement witness. The partial transpose of an entanglement witness is again 
an entanglement witness: 

F m ®(t d -—F d ). 

m 

Applying this entanglement witness on p we find that the states satisfying 

em 2 (d 2 - 1) + dm5(m - d) < 

are (PPT) entangled. 

Let us now look at the distillation properties. Similar arguments to the ones we 
used for the UUVV F-mvaiiant states apply here. For e or 5 sufficiently small, one 
can find n-undistillable states. The following vectors provide apparently the optimal 
boundaries for the 1-distillable states (we will only be interested in the states having 
5 > 0). 

1. |V) = |00)a|00) b + |10) a |10)b gives rise to 

2 + 2(de -l)/d + 4(me - l)/m + 4(1 - (m + d)e + dm6)/(md) < 

2. \ip) = |00) a |01}b + |10>a|H)b gives rise to 

e<--\. (4.13) 
m 2 

Now let us take two pairs 

„ _ -n 34 M fl 12 , de - 1 34 i2 me ~ 1 p 34 M -n 12 , 

pi - l m <g> TL d H - — lL m <g> H — — t m <g> ! d + 

1 — (m + <f)e + <im5 o A „ 19 



and 



„ _ -n 78 *->. -n 56 , " e 1 -n 78 *->. P 56 , me 1 P 78 *->. fl 56, 
a m 



1 — (m + <i)e + dm<5 
dm 



78 ^ 7756 
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Figure 4.7: The rainbow states: the top line of the triangle covers the whole spectrum 
of different types of distillability. 
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Taking both pairs together, and projecting upon Pj'^P^' 6 , we end up with a UUVVF- 
invariant state in m 2 !8> m 2 with 

, _ e(d 2 5 + d 2 -l) 
€ ~ d 2 e 2 + d 2 -l 
e 2 (d 2 -l)+d 2 5 2 
d 2 e 2 + d 2 -l ' 

Now we know when these states can be distilled. The states distillable with this pro- 
tocol are shown in Figure 14.71 One verifies that projecting upon P™ 5 <S> Pra performs 
worse. The set 1 — K + e^^f — 5 + - — 4 = contains states of all kinds: 1-distillable, 

ma ma ma ' 

n-undistillable but n + 1-distillable, NPT undistillable (conjectured), PPT bound en- 
tangled and separable. 



4.4 Activation of bound entanglement 

As bound entanglement cannot be distilled, it cannot be directly used for typical quan- 
tum information tasks, such as teleportation. Activation of bound entanglement is any 
use of bound entanglement that allows for operations which are not possible by means 
of LOCC only. 

The first of these effects was discovered by the Horodeccy 1491 ] . They showed that 



a specific state can be quasi-distilled with the aid of an unlimited supply of particular 
PPT bound entangled states, while this is not possible without them. Although this 
is a relatively weak result, it was the first indication that bound entanglement can be 
activated . Af ter this first result, a series of negative results were obtained. Linden and 
Popescu 18s| showed rigorously that some PPT bound entangled state cannot teleport 



states with better than classical fidelity. Vedral [292] subsequently demonstrated that 
the PPT distillable entanglement (that is the distillable entanglement of a state in the 
presence of an unlimited supply of bound entangled states) is still bounded by the 



relative entropy of entanglement (see also Ref. 731] ). 



2701 ]. and their results were 



Major progress on activation was made by Shor et al. 

□ 

later generalised by Eggeling et al. [87]. We discuss this last paper in some detail 
and sketch the main proofs. Their work is based on the following theorem by Rains 



2381 ] . which gives a characterisation of the maximal PPT fidelity of distillation F^. 
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This quantity is denned as the maximum overlap a state can have with the maximally 
entangled state P+ after application of a PPT and trace preserving operation. 

Theorem 4.4.1. The maximal PPT fidelity of distillation of a d-level bipartite state p 
is given by 

FAp) = maxTr(P + A(p)) = maxTr(/9j4), 
A A 

where the maximum is taken over all PPT and trace preserving maps A or all Hermitian 
operators A such that < A < 1 and -1 < dA TB < 1. 

Proof, (sketch) Since P+ commutes with all unitaries of the form U®U* we may assume 
without loss of generality that A(p) is an isotropic state 

A(p) = Tr(pB){t - P + ) + Tt(j>A)P+, 

where it is easy to see that the coefficients can always be written as linear functionals 
on p. By the definition of A, the operator A(p) has to be a proper quantum state, so 
we have to require that < A, B < 1 and (d 2 — l)B + A = 1. All that is left to do is 
account for the PPT property of A, that is if p > then 

A(p T BfB = Tr(pB TB )(t + Tr(pA TB )^ > 0. 

It is easy to see that this is equivalent to — 1 < dA Te < 1. □ 
We are now in the position of proving the main result of Eggeling et al. Ref. [s2]. 

Theorem 4.4.2. With the aid of a bound entangled state any bipartite NPT state is 
1-distillable. 

Proof. Denoting P n the projector onto the negative eigenspace of p TB we define 

A = ±(t-eP? B ), 

with < e < min(2, \\P^" B \\^~) and || • ||oo the operator norm (the maximum singular 
value). It is easy to check that the operator defined in this way meets the conditions 
of Rains' Theorem. So that 
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where M{p) > is the negativity. Therefore Rains' theorem implies that there exists a 
PPT and trace preserving operation A such that 

Tr(A(p)P + ) = Tr(p'P + ) > I. 

and such states are 1-distillable (see Theorem l4.1.3p . The proof ends by noting that any 
PPT preserving map can be implemented locally using a shared PPT bound entangled 
state (see Section I3.1.3p . □ 



17% It 



An alternative proof using entanglement witnesses can be found in Ref. 
was later also shown that any NPT state can be 1-distilled using only an infi nitesimal 
amount of bound entanglement, as this is possible for the Werner states 3101 ] . 



The implications of this result are non-trivial. As pointed out in Ref. [2701 ] . it 
probably implies the superadditivity of the distillable entanglement. Indeed, assuming 
the existence of NPT bound entanglement we have that there exist states p\ (NPT) 
and p 2 (PPT) such that E D (p x ) = E D (p 2 ) = 0, but 

E D {pi®p2) > 0, 

where the inequality follows from the previous theorem. Th is eff ect is called su peractiv a- 



271| ] (see also Ref. [80j, 273]). 



tion, and has been proven to exist in a multipartite setting 
Equally it follows from the existence of NPT bound entanglement that the distillable 
entanglement is nonconvex, as convexity would imply 

= E D (pi ® |00><00|) + E D (p 2 |11)(11|) 
> E D (pi ® |00>(00| + p 2 ® |11)(H|) = E D (p) > 0. 

The last inequality follows from the fact that the states p\ (g) |00)(00| and p 2 ® |11)(11| 
can be locally discriminated. So that when Alice and Bob are given a large number 
of states p they can recover a large supply of p\ and p 2 , which they can then distill 
by tensoring them together. In any case, suc h states p are examples of 2-distillable 



states which are not 1-distillable. In Ref. 



1301 ] Michal Horodecki brilliantly argues on 



why convexity of the distillable entanglement is not to be expected. The argument 
is based on the fact that the distillable entanglement is a measure c harac terising; the 



entanglement of source rather than that of individual states (see also 



300]). 
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Another activation effect was demonstrated by Ishizaka [1651 ] , who showed that with 
the aid of specific bound entanglement, pure entanglement is completely interconvert- 
ible: all pure entangled states can be converted into each other with some probability, 
regardless of the Schmidt rank (this was already implicit in Ref . [9] ) . Theorem 14.4.21 
implies that all NPT bound entanglement can be activated, as well as some particular 
PPT entangled states. The question remains if every PPT bound entangled state can 



2031 ] . He 



be activated. Recently, this question was answered affirmatively by Masanes 
showed that any entangled state can enhance the teleportation fidelity of some other 
state. A special case can be formulated in the context of distillation: 

Theorem 4.4.3. For any entangled state p there exists a 1-undistillable state a such 
that p® a is 1-distillable. 

Although this result is of tremendous importance, the proof works by contradiction, 
and reveals nothing about the state a. Not surprisingly, as if this was the case, it 
would solve the separability p roble m! A generalisation of this result _to mult ipartite 
systems was presented in Ref. 



2041 ] . Related results were obtained in 34j, 12971 ] , where 
it was shown that the generalised robustness (which is nonzero for entangled states) 
has an operational meaning as the maximum percentual increase an entangled state 
can provide in the fidelity of teleportation of another state. 



APPENDIX A 



Two Matlab programs 



A.l General permutation 

The m-file below is a simple implementation of an arbitrary permutation of a density 
matrix, as denned in (|1.12p . It is listed here, as some complicated implementations 
of for instance the partial transpose criterion are circulating. The procedure below is 
transparent and general. It can also be used to exchange subsystems within a density 
matrix. 

function rho2=permcrit (rho ,p) 

y o permcrit (rho , [1 2 4 3]) calculates the partial transposition 
%permcrit (rho , [2 4 13]) calculates the realigned matrix 
%permcrit (rho , [3 4 12]) exchanges the first with the second subsystem 
%permcrit (rho , [3 5 4 2 6 1]) random permutation for tripartite system 
[dl ,d2] =size(rho) ; [s ,t] =size (p) ; n=t/2; d=round(dl~ (1/n) ) ; q(p)=l:t; 
if s~=l II floor(t/2)~=ceil(t/2) || ~isequal(q(p) , 1 :t) || length(q)~=t 
error ('Not a valid permutation'); end; 

if d~n~=dl || dl~=d2, error (' Incompatible dimensions'); end; 
rho2=zeros(dl,dl) ; vec=zeros (t , 1) ; dnw=power (d,n-l : -1 : 0) ; 
for i=0:dl-l, for j=0:dl-l 
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r=l; c=l; i2=i; j2=j; 

for w=l:n, vec(w*2-l)=f loor (i2/dnw(w) ) ; vec(w*2)=f loor(j2/dnw(w)) ; 
i2=i2-vec(w*2-l) *dnw(w) ; j2=j2-vec (w*2) *dnw(w) ; end; 
vec=vec(p) ; 

for w=l:n, r=r+vec (w*2-l) *dnw(w) ; c=c+vec(w*2) *dnw(w) ; end; 
rho2(r,c)=rho(i+l, j+1) ; 
end ; end ; 

A. 2 Distillability test 

The m-file below checks whether a state is 1-distillable with the algorithm outlined 
in Section 14.2,51 Input parameters are the partial transposed state, the dimension 
of the Hilbert spaces of the parties and the number of tests. As the tests involve 
randomisation, it is recommended to set the number of tests at least to 10 4 , and to 
repeat the test a number of times. If the output is negative, however, the state is 
guaranteed to be distillable. 

function mineig=distest(rhotb, d, nrtests) 
idxs=zeros (1) ; mineig=l; precision=-10~ (-8) ; 
for yy=l:d 

idxs(2*yy-l)=(yy-l)*d+l; idxs(2*yy)=(yy-l)*d+2; 
end; 

if issparse (rhotb) , EE=sparse(d*d,d*d) ; 
else, EE=zeros(d*d,d*d) ; end; 
for aa=l: nrtests 

[UU,R] =qr(randn(d)+i*randn(d) ) ; 

UU=UU*diag(exp(2*pi*i*rand(d,l))) ; UU(3:d, :)=0; 

for k=l:d:d*d, EE(k:k+d-l ,k:k+d-l)=UU; end; 

B=EE*rhotb*EE' ; mreig=min(real (eig(f ull (B(idxs , idxs) ) ) ) ) ; 

if mreig<mineig, mineig=mreig; UUb=UU; end; 

if mineig<precision, return; end; 
end; 
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for loop=l :nrtests/d 
for counter=l:2 

tel=UUb(counter , : ) ; 

for mn=l : d 

UUb (counter , mn) =randn ( 1 , l)+i*randn(l , 1) ; 
UUb(counter , : )=UUb( counter, : )/norm(UUb(counter, :)) ; 
for k=l:d:d*d, EE(k:k+d-l ,k:k+d-l)=UUb; end; 
B=EE*rhotb*EE' ; mreig=min(real (eig(full (B(idxs , idxs) ) ) ) ) ; 
if mreig<mineig, mineig=mreig; tel=UUb (counter , :) ; 
else, UUb(counter , : )=tel ; end; 
end; 

if mineig<precision, return; end; 
end ; end ; 
mineig=l ; 
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